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PREFACE 


The  CURVES  computer  program  described  in  this  report  is  an  out- 
growth of  Rand's  military  cost  analysis  research  activity.  The  program 
provides  a user-oriented  tool  for  estimating  by  least-squares  procedures 
the  parameters  and  statistical  characteristics  of  several  equations  com- 
monly used  in  the  derivation  of  cost-estimating  relationships. 

A central  feature  of  cost  analysis  research  is  the  development  of 
predictive  relationships  by  which,  for  example,  the  costs  of  new  mili- 
tary equipment  or  activities  can  be  estimated  from  data  on  past  equip- 
ment and  activities.  Although  several  standard  computer  programs  are 
available  for  curve-fitting  and  statistical  analysis  of  data,  they  often 
are  not  well  suited  to  the  purposes  of  cost  analysis.  Most  standard 
statistical  programs  are  designed  to  accommodate  very  large  data  sets 
and  to  provide  a wide  range  of  appropriate  statistical  tests.  The  typi- 
cal cost  analysis  application,  however,  involves  few  data  points  (usually 
less  than  100)  and  requires  rather  selective  curve-fitting  routines  and 
statistical  tests,  circumstances  that  make  the  standard  programs  both 
cumbersome  and  fairly  expensive  to  operate.  Moreover,  no  single  standard 
statistical  program  is  likely  to  include  all  of  the  functional  forms 
most  useful  to  the  cost  analyst.  Hence,  a special-purpose  program  such 
as  CURVES  is  both  more  convenient  and  more  economical  to  operate  than  a 
standard  program. 

The  CURVES  program  described  here  is  a substantial  updating  and  ex- 
tension of  the  program  in  an  earlier  Rand  report:  RM-5762-PR,  CURVES: 

A Five-Function  Curve-Fitting  Computer  Program  (December  1968) , by  one 
of  the  present  authors,  H.  E.  Boren,  Jr.  The  new  CURVES  program  adds 
three  functions  (logarithmic-linear  and  two  semilogarithmic-linear  forms) 
and  several  new  statistical  and  operational  features. 

The  CURVES  program  was  developed  as  a by-product  of  research  on  the 
cost  of  advanced  military  aircraft  and  missiles.  Its  use  is  not  re- 
stricted to  advanced  hardware,  however,  or  even  to  cost  analysis  appli- 
cations. It  should  be  useful  to  analysts  throughout  the  Air  Force  and 
elsewhere  in  the  Defense  Department  who  are  concerned  with  describing 
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causal  relationships  in  functional  form.  This  report  was  undertaken 
as  part  of  the  Project  RAND  research  task  entitled  "Cost  Analysis 
Methods  for  Air  Force  Systems." 

Every  effort  has  been  made  to  remove  errors  from  the  CURVES  pro- 
gram described  in  this  report.  However,  no  guarantee,  expressed  or 
implied,  is  made  as  to  either  the  numerical  or  the  logical  accuracy 
of  the  program.  Information  concerning  any  errors  or  difficulties 
found  in  the  use  of  the  program  or  documentation  will  be  greatly 
appreciated  by  the  authors. 

During  the  period  when  this  report  was  prepared.  Captain  Gerald 
W.  Corwin  was  on  duty  at  The  Rand  Corporation  in  the  Management  Sciences 
Department.  He  is  at  present  with  the  Cost  Analysis  Division,  Direc- 
torate of  Management  Analysis,  Office  of  the  Comptroller  of  the  Air 
Force,  Headquarters  United  States  Air  Force. 

This  report,  supersedes  R-1753-PR,  which  was  first  distributed  in 
December  1975.  It  is  being  reissued  in  its  present  form  to  acquaint 


users  with  several  recent  major  modifications  to  the  CURVES  program. 
The  modifications  are  discussed  in  Appendix  C,  and  the  updated  program 
listing  is  provided  in  Appendix  D. 
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SUMMARY 

This  report  describes  an  extension  of  a FORTRAN-IV  curve-fitting 
(regression  analysis)  computer  program  (CURVES)  developed  in  1968  to 
facilitate  the  derivation  of  cost-estimating  relationships  for  advanced 
military  equipments.  The  program  makes  least-squares  determinations  of 
the  parameters  of  any  of  eight  types  of  equations  selected  by  the  user, 
given  a set  of  observations  on  the  dependent  and  independent  variables 
of  interest.  The  types  of  equations  that  can  be  fitted  using  CURVES 
are:  linear,  quadratic,  power,  asymptotic-power,  exponential,  log- 

arithmic-linear, and  two  types  of  semilogarithmic-linear . Except  for 
the  quadratic  and  asymptotic-power  equations,  up  to  seven  independent 
variables  may  be  used. 

Because  of  the  importance  of  learning  curve  theory  in  cost  analysis 
work,  CURVES  has  been  expanded  to  include  logarithmic-linear  as  well  as 
semilogarithmic-linear  functions.  Although  the  power  form  counterpart 
of  the  logarithmic-linear  equation  was  included  in  the  previous  version, 
estimates  of  the  parameters  obtained  from  the  regressed  logarithmic- 
linear  equation  are  often  preferred  to  those  obtained  from  the  regressed 
power  equation. 

Other  new  features  have  been  added.  A correlation  matrix  of  the 
input  data  is  provided  for  all  fitted  equations  using  more  than  one  in- 
dependent variable,  as  well  as  a variance-covariance  matrix  of  the  esti- 
mated coefficients.  Also  included  are  standard  errors  and  Student's 
t-ratios  of  the  parameters,  significance  levels,  beta  coefficients,  and 
the  Durb in-Watson  statistic.  A plot  routine  is  incorporated  for  provid- 
ing various  plots  of  the  data.  The  plots  available  to  the  user  are: 

(a)  residual  Y versus  fitted  Y,  (b)  observed  Y versus  fitted  Y,  and 
(c)  observed  Y versus  any  one  of  the  independent  variables.  In  the 
case  of  a single  independent  variable,  plot  (c)  can  be  used  to  provide 
a plot  of  the  regression  equation.  The  program  is  fairly  small  (about 
92,000  bytes  of  core),  fast  in  execution  time,  and  hence  cheap  to  operate. 

Section  I describes  the  overall  features  of  CURVES  and  the  options 
available  to  the  user.  Section  II  discusses  the  equations  available  for 
regression  in  the  program,  including  an  examination  of  the  nonlinear  ones. 
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which  require  special  methods  for  solution.  The  section  also  contains 
a brief  summary  of  the  statistics  used  in  the  program.  Section  III 
presents  specific  details  of  the  operation  of  the  program,  including 
the  options  available  to  the  user,  Section  IV  describes  the  program 
outputs.  A sample  output  for  four  runs  is  shown,  together  with  a card 
listing  of  the  deck  setup  required  for  the  runs. 

Appendices  A and  B treat  the  mathematical  considerations  relating 
to  nonlinear-least-squares  solutions.  Appendix  C discusses  recent 
modifications  that  have  been  made  to  CURVES,  and  Appendix  D presents 
an  updated  listing  of  the  CURVES  computer  program. 
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I.  INTRODUCTION 


One  of  the  central  tasks  of  cost  analysis  is  the  development  of 
cost-estimating  relationships,  predictive  models  that  mathematically 
describe  the  cause  and  effect  connections  affecting  the  resources 
(costs)  required  to  produce  given  outputs.  In  military  cost  analysis, 
one  often  needs  to  estimate  approximate  relationships  that  describe 
the  cost  of  advanced  aircraft  in  terms  of  weight,  speed,  production 
rate,  and  similar  parameters.  One  of  the  several  means  by  which  esti- 
mating relationships  may  be  derived  is  through  application  of  curve- 
fitting and  statistical  analysis  techniques  to  empirical  (historical) 
data.  The  CURVES  computer  program  described  in  this  report  is  specifi- 
cally designed  to  provide  the  computational  facility,  mathematical 
equations,  and  descriptive  statistics  most  often  needed  in  cost  analysis. 

The  widespread  use  of  similar  techniques  in  many  fields  of  research 
has  led  to  development  of  numerous  standard  curve-fitting  and  statistical 
computer  programs.  Most  of  these  standard  programs,  which  provide  many 
desirable  features  and  options,  often  are  not  well  suited  to  the  needs 
of  cost  analysts.  The  typical  cost  analysis  problem  involves  small 
amounts  of  data,  often  fewer  than  100  data  points.  Most  standard  pro- 
grams are  intended  to  accommodate  very  large  data  sets.  One  consequence 
is  that  they  frequently  are  cumbersome  and  expensive  to  operate,  involv- 
ing large  computer  core  storage  needs  and  frequent  disk  operations,  when 
applied  to  small  data  sets.  Moreover,  such  programs  generally  do  not 
treat  the  nonlinear  equations  encountered  in  cost  analysis;  and  the 
CURVES  program  is  intended  to  provide  a more  economical,  more  convenient, 
and  better  tailored  tool  for  cost  analysis.  Although  it  was  developed 
in  the  context  of  research  on  the  cost  of  advanced  military  aircraft  and 
missiles,  it  is  applicable  to  a much  wider  range  of  situations  in  which 
mathematically  described  causal  relationships  are  needed. 

The  CURVES  program  was  written  originally  for  the  purpose  of  making 
ordinary  least-squares  determinations  for  the  parameters  of  five  types 
of  equations:  linear,  quadratic,  asymptotic-power,  and  exponential.^ 

^H.  E.  Boren,  Jr.,  CURVES:  A Five-Function  Curve-Fitting  Computer 

Program , The  Rand  Corporation,  RM-5762-PR,  December  1968. 
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This  report  presents  an  updated  version  of  CURVES  that  is  much  more 
powerful  and  flexible,  yet  cheaper  to  operate,  than  the  previous  ver- 
sion. In  addition  to  the  original  five  equations,  three  logarithmic 
equations  have  been  added.  They  are:  Ln  Y vs  Ln  X (logarithmic- 

linear)  , Ln  Y vs  X (semilogarithmic-linear) , and  Y vs  Ln  X (semilog- 
arithmic-linear ) . "*■  Except  for  the  quadratic  and  asymptotic-power 
equations,  all  equations  may  now  be  fitted  using  up  to  seven  independ- 
ent variables.  In  addition,  values  of  the  Y-intercept  may  be  prespeci- 
fied for  all  equations  except  the  power  and  exponential.  The  program 

has  been  rewritten  to  be  as  user-oriented  as  possible  in  terms  of  input 

2 

procedures.  CURVES  is  a fairly  small  program  (about  92,000  bytes  of 

core),  is  very  fast  in  execution  time,  and  is  designed  to  minimize  disk 

input/output  operations.  It  is  currently  in  use  on  the  Rand  IBM  370/158 

computer  but  is  adaptable  to  any  computer  system  that  accepts  FORTRAN 

and  has  enough  core  capacity  to  handle  the  program. 

A plot  routine  has  been  added  to  provide  the  following  plots: 

(a)  residual  Y versus  fitted  Y,  (b)  observed  Y versus  fitted  Y,  and 

(c)  observed  Y versus  any  one  of  the  independent  variables.  The  routine 

may  also  be  used  to  plot  the  regression  equation  for  a one-independent- 

variable  case.  Plots  use  letter  and  numeral  symbols  so  that  each  data 

point  may  be  easily  identified. 

The  statistics  calculated  in  CURVES  include  those  relating  to 

"goodness-of-f it  measures,"  such  as  sum  of  squares  of  Y residuals,  total 

2 

sum  of  squares,  coefficient  of  determination  (R  ),  standard  error  of 
estimate  of  Y (SEY),  coefficient  of  variation  (ratio  of  SEY  to  sample 
mean  of  Y),  mean  of  absolute  relative  deviations  of  Y,  and  the  F- 
statistic.  Also  included  are  standard  errors  of  the  parameter  estimates, 

3 

Student's  t-ratios,  significance  levels,  beta  coefficients,  and  the 
Durbin-Watson  statistic.  The  printbut  of  significance  levels  is  very 

"'"In  this  report  the  abbreviation  "Ln"  is  used  to  denote  a natural 
logarithm  (to  base  e = 2.71828...). 

2 

This  includes  about  25,000  bytes  for  two  variable-dimensioned 
arrays,  whose  sizes  can  be  changed  to  suit  the  user's  needs. 

3 

Formulas  for  calculating  significance  levels  were  obtained  from 
a study  at  Rand  by  D.  Tihansky  and  F.  Timson  in  April  1972. 
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useful  because  it  obviates  the  need  to  obtain  the  values  from  a Student's 
t-table.  Means  and  standard  deviations  of  the  dependent  and  independent 
variables  are  printed  as  well  as  a correlation  matrix  in  the  multivariate 
case.  The  variance-covariance  matrix  of  the  estimated  coefficients  is 
now  printed. 

CURVES  can  treat  up  to  several  hundred  data  points  depending  on  the 
type  of  equation  being  fitted,  the  number  of  independent  variables  being 
used,  and  whether  plotting  is  done.  A set  of  data  needs  to  be  entered 
only  once  even  if  several  regressions  are  to  be  run  on  it.  A variable- 
format  procedure  is  provided  the  user  so  that  data  may  be  entered  in  any 
order  on  the  input  cards.  Data  may  also  be  entered  from  tape  or  disk 
provided  that  the  data  are  in  the  appropriate  format.  CURVES  also  pro- 
vides for  variable  transformations  as  discussed  in  Appendix  C. 

The  CURVES  program  is  written  in  FORTRAN-IV  (G/H  level)  except  for 
one  assembler  subroutine  that  permits  reading  the  data  from  memory;  how- 
ever, the  subroutine  is  not  required  for  normal  program  operation.  If 
the  assembler  subroutine  cannot  be  made  available,  the  user  can  use  a 
"scratch"  disk  for  data  storage  and  retrieval. 

One  major  change  in  this  edition  of  CURVES  is  that  the  Gauss-Jordan 
method  of  solving  simultaneous  equations  is  used  instead  of  Cramer's 
Rule.  All  solutions  are  made  in  double-precision  arithmetic  either 
through  standard  algebraic  methods  for  the  linear,  quadratic,  and  log- 
arithmic equations  or  through  iterative  methods  for  the  other  equations. 
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II.  PROGRAM  CONSIDERATIONS 


EQUATION  TYPES 

The  equations  available  in  CURVES  were  chosen  principally  on  the 
basis  of  their  application  to  the  derivation  of  cost  analysis  estimat- 
ing relationships.  The  Y- intercept  value  A (or  Ln  A for  the  logarithmic- 
linear  equation)  may  be  specified  for  all  equations  except  the  power  and 
exponential.^  The  equations  are: 


1.  Linear 

Y=A+B  • XI  +C  • X2+  ...  + H • X7, 

2 . Quadratic 

Y=A+B*X1+C*  XI2, 

3.  Power 

Y = A • X1B  * X2C  • ...  • X7H, 

4.  Asymptotic-Power 

Y = A + B • X1C, 

5.  Exponential 

Y=e(A+B  * XI  + C • X2+  ...  + H • X7), 

6 . Logarithmic-Linear — Ln  Dependent  vs  Ln  Independent 

Ln  Y = Ln  A + B • Ln  XI  + C • Ln  X2  + ...  + H • Ln  X7, 

7.  Semilogarithmic-Linear — Ln  Dependent  vs  Independent 

Ln  Y = A + B • XI  + C • X2  + ...  + H • X7 , 

8.  Semilogarithmic-Linear — Dependent  vs  Ln  Independent 
Y = A + B • Ln  XI  + C • Ln  X2  + ...  + H • Ln  X7 , 

*ln  this  report  the  dependent  variable  is  represented  by  Y and 
the  independent  variables  by  XI,  X2,  ...,  X7.  When  only  one  independ- 
ent variable  is  considered,  XI  is  used. 
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where 

XI,  X2,  . 
A , B , C , 


Y = dependent  variable, 

. , X7  = independent  variables, 

H = parameters  to  be  estimated  by  least-squares  methods, 
Ln  = natural  logarithm,  base  e. 
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EQUATION  CHARACTERISTICS 

Examples  of  types  of  curves  that  can  be  fitted  in  the  program  are 
shown  in  Fig.  1. 

Linear  (1) 

The  linear  form  is  the  simplest  treated  here.  Its  characteristics 
are  well  known  and,  in  our  opinion,  need  no  further  elaboration.  The 
user  has  the  option  of  specifying  the  Y-intercept  A. 

Quadratic  (2) 

Sometimes  the  quadratic  equation  is  used  to  represent  points  that 
lie  along  a parabola.  However,  one  must  be  aware  that  a quadratic  equa- 
tion always  has  a maximum  or  minimum  point  (vertex) . This  means  that 
the  effects  of  changes  in  the  independent  variable  (XI)  on  the  dependent 
variable  (Y)  are  different  in  sign  on  either  side  of  the  vertex.  The 
coordinates  of  the  vertex  are  printed  in  the  output.  Again,  the  user 
has  the  option  of  specifying  the  Y-intercept  A. 

Power  (3) 

The  power  equation  is  one  of  the  more  common  equations  used  in  cost 
analysis  work.  A plot  of  its  logarithmic  counterpart,  the  logarithmic- 
linear  form,  is  known  as  the  "learning  curve"  or  "improvement  cost  curve." 
Conceptually,  a regression  of  the  power  form  does  not  result  in  the  same 
estimates  of  the  parameters  as  a regression  of  its  logarithmic-linear 
counterpart.  Appendix  A discusses  the  differences  between  the  power 
and  logarithmic-linear  regressions.  For  positive  exponent  B,  the  graph 
of  the  power  equation  always  passes  through  the  origin.  Therefore,  it 
should  never  be  used  if  a positive  Y-intercept  is  desired  or  logically 
required.  For  negative  B,  the  equation  is  undefined  at  XI  = 0,  has  a 
negative  slope,  and  approaches  zero  asymptotically  as  XI  becomes  infinite. 
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Asymptotic-Power  (4) 

An  examination  of  the  graph  of  the  asymptotic-power  equation  shows 
that  the  curve  has  a horizontal  asymptote  of  Y = A for  negative  C.  That 

C 

is,  as  XI  becomes  large,  the  second  term  (B  • XI  ) approaches  zero,  and 
hence  the  value  of  Y tends  to  A.  Consequently,  there  is  a leveling-off 
effect  for  negative  C.  This  equation  may  thus  be  used  to  represent 
points  that  lie  along  a curve  either  increasing  or  decreasing  to  a hori- 
zontal asymptote.  Like  thv  power  equation,  this  equation  is  undefined 
at  XI  = 0 for  negative  C.  For  positive  C,  the  Y-intercept  is  equal  to 

Q 

A.  As  XI  becomes  large,  the  second  term  (B  • XI  ) ultimately  becomes 
large  compared  with  A,  and  therefore  the  equation  behaves  like  the  power 

C 

function  (B  • XI  ) in  this  region  of  XI. 

Exponential  (5) 

The  exponential  form  is  used  to  represent  points  that  lie  along  a 

A 

curve  having  a positive  Y-intercept  (e  ) . As  XI  increases,  the  graph 
of  the  equation  rises  for  B > 0 and  falls  for  B < 0.  In  the  latter 
case,  the  Xl-axis  is  the  asymptote  of  the  curve.  The  logarithmic  coun- 
terpart of  the  exponential  equation  is  the  semilogarithmic  equation  in 
which  Ln  Y is  a linear  function  of  XI.  However,  for  the  same  reasons 
as  stated  for  the  power  and  logarithmic-linear  cases,  a least-squares 
regression  of  the  exponential  form  produces  different  estimates  of  the 
parameters  than  a regression  of  its  semilogarithmic-linear  counterpart. 

Logarithmic-Linear  (6) 

The  logarithmic-linear  equation,  also  known  as  the  learning  curve, 
is  an  equation  in  which  the  logarithm  of  the  dependent  variable  is  a 
linear  function  of  the  logarithms  of  the  independent  variables.  In  the 
program,  it  is  fitted  using  the  same  technique  as  for  the  linear  equa- 
tion. The  constant  term  may  be  specified  in  the  form  of  Ln  A. 

The  "learning"  process  is  a phenomenon  that  prevails  in  many  in- 
dustries, and  its  existence  has  been  verified  by  empirical  data  and 
controlled  tests.  The  basis  of  learning  curve  theory  is  that  each  time 
the  total  quantity  of  items  produced  is  increased  by  a constant  percent- 
age, the  cost  per  item,  or  average  cost  of  all  items  produced,  is  reduced 
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by  some  constant  percentage.  If  the  number  of  items  produced  is 
doubled,  then  the  percentage  to  which  the  cost  is  reduced  is  known 
as  the  learning  curve  slope.  For  example,  if  the  number  of  items  is 
increased  from  120  to  240,  and  the  cost  reduces  from  $100  to  $80,  then 
the  learning  curve  slope  is  80  percent.  Equations  of  this  type  may  be 
applied  either  to  the  cost  per  Nth  item  produced  (unit  cost  curve)  or 
to  the  average  cost  for  the  first  N items  produced  (cumulative  average 
cost  curve). 

The  learning  curve  slope  (S)  can  be  expressed  as  a function  of  the 
exponent  (B)  as  follows.  From  the  above  definition 


S 


S 


A • (2N)B 
A • NB 


or 


S = 2 


LojB_S 
Log  2 ’ 


where  S 
B 
N 
Y 

Log 


learning  curve  slope  (decimal) , 
exponent  of  quantity, 
quantity, 

dependent  variable  (cost,  manhours,  people,  etc.), 
logarithm  to  any  base. 


With  the  use  of  the  power  and  logarithmic  functions  and  the  plot 
routine  provided  in  CURVES,  the  user  now  has  the  capability  to  examine 
learning  curve  regressions  and  select  the  equation  most  appropriate  for 
the  set  of  data  under  study. 
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Semilogar ithmic  (7,  8) 

The  semilogarithmic  equations  exist  in  two  forms.  In  one  form, 
the  logarithm  of  the  dependent  variable  is  a linear  function  of  the 
independent  variables.  In  the  other,  the  dependent  variable  is  a linear 
function  of  the  logarithms  of  the  independent  variables.  The  graphs 
of  the  two  equations  produce  straight  lines  on  rectangular  coordinate 
paper  or,  in  terms  of  X and  Y,  straight  lines  on  logarithmic  paper, 
provided  that  the  proper  axis  is  scaled  in  logarithms.  The  constant  A 
may  be  prespecified  in  fitting  equations  using  these  functional  forms. 

Note  that  for  the  logarithmic-linear  equation  (Eq.  (6)),  the  constant 
term  is  specified  as  Ln  A,  whereas  for  the  semilogarithmic  cases,  it 
is  specified  as  A. 

NONLINEAR-LEAST-SQUARES  ESTIMATES 

Least-squares  estimates  of  the  parameters  of  an  equation  are  always 
unique  with  a closed,  algebraic  solution  provided  the  equation  is  linear 
with  respect  to  all  of  its  parameters.  Therefore,  for  this  program, 
regressions  of  the  linear,  the  quadratic,  and  the  three  logarithmic 
equations  produce  such  estimates  of  the  parameters.'*'  However,  the  power, 
asymptotic-power,  and  exponential  equations  are  not  linear  in  terms  of 
all  of  their  parameters.  Thus,  least-squares  estimates  of  their  param- 
eters cannot  usually  be  obtained  by  simple,  algebraic  methods  and,  as 
shown  later,  may  not  represent  absolute  minimums  of  the  sum  of  squares 
of  Y residuals.  They  must  be  obtained  in  some  other  way,  usually  through 
some  type  of  iterative  procedure.  (The  general  principles  of  such  pro- 
cedures and  other  mathematical  considerations  relating  to  nonlinear- 
least-squares  regressions  are  presented  in  Appendix  A.) 

For  the  power  and  exponential  equations,  a modified  Gauss-Newton 
method  is  used,  in  which  initial  estimates  of  the  parameters  are  obtained 
from  the  logarithmic-linear  regressions;  and  then  corrections,  guaranteed 
to  produce  convergence  to  a solution,  are  applied  to  those  initial  esti- 
mates. For  all  iterative  procedures,  a solution  is  reached  when  the 

^Under  this  definition,  the  quadratic  function  is  considered  to  be 
a linear  function  because  it  is  linear  with  respect  to  the  parameters 
that  are  to  be  estimated. 
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ratio  of  the  value  of  each  parameter  to  its  value  corresponding  to 
the  previous  iteration  differs  from  unity  by  some  predetermined  value 
(^10  ^ or  as  otherwise  specif ied) . ^ 

Least-squares  estimates  of  the  parameters  of  the  asymptotic-power 
equation  are  based  on  another  type  of  iterative  procedure  because  there 
appears  to  be  no  easy  way  to  obtain  the  initial  guesses  of  the  parameters 
required  for  the  modified  Gauss-Newton  method.  This  procedure  is  treated 
in  Appendix  B.  Because  the  modified  Gauss-Newton  method  cannot  be  used 
in  this  case,  the  equation  is  restricted  to  one  independent  variable. 

STATISTICAL  CONSIDERATIONS 

Program  statistics  contain  the  standard  measures  relating  to  good- 
ness of  fit,  such  as  sum  of  squares  of  residuals,  sum  of  squares  total, 

2 

coefficient  of  determination  (R  ),  standard  error  of  estimate  of  Y (SEY), 
coefficient  of  variation  (ratio  of  SEY  to  sample  mean  of  Y),  mean  of 
absolute  relative  deviations  of  Y,  F statistic,  and  the  Durbin-Watson 
statistic.  The  relative  deviation  of  Y at  the  ith  point  is  the  ratio 
of  the  Y residual  at  that  point  to  the  observed  value  of  Y.  Also  in- 
cluded are  standard  errors  of  the  parameters,  t-ratios,  significance 
levels,  and  beta  coefficients.  Means,  standard  deviations,  and  corre- 
lation coefficients  are  printed  for  the  input  data. 

The  standard  errors  of  the  parameters  for  the  nonlinear  cases  (power, 
asymptotic-power,  and  exponential)  are  calculated  as  follows.  In  obtain- 
ing the  least-squares  estimate  of  the  parameters  for  the  power  and  ex- 
ponential equations,  a matrix  of  partial  derivatives  with  respect  to  the 
parameters  is  calculated  and  inverted  during  each  iteration  in  order  to 
obtain  correlations  to  the  parameter  values.  As  the  iterations  converge, 
the  square  roots  of  the  diagonal  terms  of  the  inverted  matrix  multiplied 
by  the  standard  error  of  estimate  of  Y converge  to  the  standard  errors 
of  the  parameters  obtained  by  the  least-squares  estimates.  For  the 

^This  procedure  is  described  in  detail  in  C.  A.  Graver  and  H.  E. 
Boren,  Jr.,  Multivariate  Logarithmic  and  Exponential  Regression  Models, 

The  Rand  Corporation,  RM-4879-PR,  July  1967.  The  term  "exponential" 
there  is  equivalent  to  the  term  "power"  in  this  report. 
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I 

asymptotic-power  equation  for  which  the  Gauss-Newton  method  is  not 
used,  the  inverted  matrix  of  partial  derivatives  is  obtained  only  after 
the  estimates  of  the  parameters  are  computed  using  another  type  of 
iterative  procedure  discussed  in  Appendix  B. 

The  Durbin-Watson  statistic  is  used  to  test  for  serial  correla- 
tion.''' It  is  based  on  successive  differences  in  the  Y residuals. 
Therefore,  the  statistic  may  not  be  useful  unless  the  data  are  ordered 
in  some  meaningful  way. 

The  statistics  that  are  printed  out  by  CURVES  are  defined  in 
Table  1. 

^J.  Durbin  and  G.  S.  Watson,  "Testing  for  Serial  Correlation  in 
Least  Squares  Regression  II,"  Biometrika,  Vol.  38,  1951. 


Table  1 

STATISTICAL  EQUATIONS  USED  IN  PROGRAM 


Statistic 

Equation 

Degrees  of  freedom  for  error 

DF1  = N - M 

where  N = number  of  data  points 

M = number  of  parameters  estimated 

Degrees  of  freedom  due  to 
regression 

DF2  = number  of  independent  variables 

Total  degrees  of  freedom 

DFT  = DF1  + DF2 

Sum  of  squares  total 

(unspecified  Y-intercept) 

N 2 

SST  = ^2  (Y±-  Y)  , 

i=l 

where 

Y.  = Y value  for  ith  observation 

l 

Y = mean  of  observed  Y values 

Sum  of  squares  total 

(specified  Y-intercept  A) 

N 2 
SST  = ^2  <V  A)  » 

i=l 

Sum  of  squares  of  residuals 

SSE  - £ (Yrvcl)2, 
i=l 

where 

Yci  = fitted  value  for  ith 
observation 

, 

Standard  error  of  estimate 


S 


SSE 

T5FT 


Statistic 


Equation 


Coefficient  of  variation 


Coefficient  of  determination 


2 _ SSE 
R 1 SST 


Relative  deviation  for 
ith  residual 


Di  ■ -tt1  • q * 0 


Mean  of  absolute  relative  Y 
deviationsb 


Standard  deviation  of 
input  variables 


Student's  t-ratio  of 
parameter  estimates 


f ^ 

E (vrv); 


SDEV  =\ 


where 


= value  of  input  variable 
under  consideration 


V = mean  value  of  input  variable 


SE  ’ 


where 


B = parameter  estimate 


SE  = standard  error  of  estimated 
parameter 


I 


I 

I 

' 


I 

I 


* 
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Table  1 — (cont.) 


Statistic 


Level 

ing 


of  significance 
to  t-ratio 


relat- 


SIGLEV  = 


Equation 


(DF1+1)  / 2 
) dX  , 


where  T denotes  the  gamma  function. 

(Above  formula  approximated  by  series 
expansion  as  used  in  a study  by  D. 
Tihansky  and  F.  Timson  at  Rand  in 
April  1972.) 


3-coefficient 

coefficient 


of  regression 
on  variable  V 


SDEVV 

~s~devy 


Durbin-Watson  statistic 


DW 


E KW'Wi-rLi-i)!2 


SSE 


| 


, 

{ 


For  the  logarithmic-linear  case,  the  Y-intercept  is  Ln  A. 

If  any  Yj  is  zero,  the  corresponding  Dj.  cannot  be  calculated. 

In  such  a case  the  summation  is  reduced  to  fewer  than  N data  points. 
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III.  INPUT  PROCEDURES 


The  flow  of  operations  within  the  program  is  depicted  in  Fig.  2. 

The  program  is  structured  so  that  many  sets  of  data  may  be  entered  in 
which  each  set  constitutes  a run.’'-  Table  2 lists  the  maximum  number 
of  data  points  that  can  be  used  for  each  equation.  As  was  stated  pre- 
viously, the  maximum  number  depends  on  the  type  of  equation  being  re- 
gressed, the  number  of  independent  variables  being  used,  and  whether 
plots  are  to  be  made.  As  soon  as  each  set  is  read  in,  the  program 

operates  on  that  set  before  proceeding  to  the  next  set  of  input  data. 

2 

Each  set  of  data  may  be  entered  on  a separate  deck  of  cards.  However, 
several  or  all  of  the  sets  of  data,  if  space  on  the  cards  permits,  may 
be  entered  on  one  deck  of  cards,  thus  effecting  considerable  savings 
in  the  use  of  cards  and  in  the  effort  of  duplicating  a deck  of  cards 
containing  data  for  several  runs.  A variable  format  procedure  is  used, 
allowing  much  flexibility  in  the  format  of  the  input  data.  Data  may 
also  be  entered  from  disk  or  tape  provided  all  format  requirements  are 
met. 

Table  3 lists  the  types  of  cards  used  for  a job  consisting  of  one 
or  more  runs  (steps  6-10  may  be  repeated  as  often  as  desired). 

TITLE  CARD 

The  Title  card  must  be  entered  for  each  run.  It  contains  the  title 
for  the  current  run  (which  is  listed  at  the  top  of  each  output  page)  and 
may  consist  of  any  valid  characters;  all  80  columns  may  be  used.  An 
example  of  a Title  card  is  shown  in  Fig.  3 as  it  might  appear  in  a card 
arrangement  for  a CURVES  run. 

CONTROL  CARD 

The  Control  card  is  the  second  input  card  for  each  regression  and 

^A  series  of  one  or  more  runs  is  defined  here  as  a job  constituting 
one  session  on  the  computer. 

2 

The  use  of  the  word  card  in  this  report,  while  usually  meaning  the 
normal  punch  card,  can  in  the  case  of  data  "cards"  or  blank  "card"  mean 
physical  record  (card,  disk,  or  tape). 


Read 


of  operations 
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Table  2 

MAXIMUM  NUMBER  OF  DATA  POINTS  BY  EQUATION,  PLOT  OPTION, 

AND  NUMBER  OF  INDEPENDENT  VARIABLES 


(: 


Equation 

Number  of 
Independent 
Variables 

Without 

Plotting 

With 

Plotting 

Without  ID3 

With  ID 

Without 

ID  With  ID 

Linear 

1 

551 

441 

381 

306 

3 

367 

314 

254 

218 

7 

220 

200 

152 

138 

Quadratic15  ] 

1 

441 

367 

306 

254 

Exponential  ! 

3 

314 

215 

218 

190 

Semi logarithmic-linear 

7 

200 

183 

138 

126 

(Ln  dependent  variable) 

Power 

1 

367 

314 

254 

218 

Asymptotic  power'5 

3 

220 

200 

152 

138 

Ln-linear 

7 

121 

115 

84 

79 

Semilogarithmic- linear 

(Ln  independent  variable) 

ID  indicates  data-point  identifier. 
^For  one  independent  variable  only. 


NOTE:  This  table  is  based  on 


where 

N = maximum  number  of  data  points 
max  r 

C = maximum  number  of  cells  of  data  storage  = 2211 

0^  = number  of  cells  reserved  for  plotting  = 676 

p = plotting  requirement  (p  = 0 for  no  plotting;  p = 1 for  plotting) 

k = number  of  columns  required  for  data;  columns  are  required  for: 

(a)  dependent  variable  Y (always) 

(b)  each  independent  variable  (always) 

(c)  identifier  (only  if  used;  see  Control  Card, 
columns  2 to  10) 

(d)  computed  Y (always) 

(e)  residual  Y (always) 

(f)  Ln  Y for  all  equations  except  1 and  2 

(g)  Ln  X for  each  independent  variable  for  equations 
3,  4,  6,  8 


J 


Table  3 

CARD  TYPES  FOR  CURVES  PROGRAM 


First 

Run 

1. 

Title 

Required. 

2. 

Control  card 

Required . 

3. 

Specification 

cards 

Required;  used  to  supply  special  control 
information  such  as  format,  labels,  data 
location,  and  initial  guesses  for  iter- 
ative solutions.  A "Read"  or  "Read8" 
specification  card  is  always  required. 

4. 

Data  cards 

Required . 

5. 

Blank  card 

Required;  used  only  at  end  of  data  cards. 

Succeeding  Runs 

6. 

Title  card 

Required. 

7. 

Control  card 

Required. 

8. 

Specification 

card  (s) 

"Read"  or  "Read8"  card  is  always  required. 

9. 

Data  cards 

Required  only  if  input  was  not  saved  from 
previous  run. 

10. 

Blank  card 

Required  only  at  end  of  input  data  cards. 

After 

Final  Run 

11. 

Done  card 

Word  DONE  entered  in  Cols.  1-4  after  last 
run.  Terminates  job. 

is  required  (even  if  it  is  blank).  It  contains  the  control  data  to 
perforin  the  desired  regression.  Table  4 summarizes  all  the  information 
on  the  Control  card.  An  example  of  a control  card  is  shown  in  Fig.  3. 


Equation  Designator 

Column  1 is  used  for  the  equation  designator.  An  integer  from  1 
through  8 is  entered  to  designate  which  equation  is  to  be  used  for  the 
regression  in  the  run.  (A  zero  indicates  a plot-only  option.)  The 
equation  designators  are: 


Fig.  3 — Examples  of  Title,  Control,  Specification,  and  Data  cards 
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Table  4 

SUMMARY  OF  INFORMATION  ON  CONTROL  CARD 


Card 

Co  1 umn ( s ) 

Use 

Value 

Control  Description 

1 

equation  index 

Blank 

Use  previous  equation  index. a 

0 

No  regression,15  only  plot  Y vs  XM,  where  M = 1,  2,  ....  7 

1 

Linear  regression 

2 

Quadratic  regression 

3 

Power  regression 

4 

Asymptotic-power  regression 

5 

Exponential  regression 

6 

Ln-linear  regression 

7 

Semi  log-1  inear  (Ln  dependent  vs  independent)  regression 

8 

Semi  log-1  inear  (dependent  vs  Ln  independent)  regression 

2-10 

Variable  order 

I 

Identifier  (optional) 

Y 

Dependent  variable  (required) 

1 

First  independent  variable  (required) 

2-7 

Second  through  seventh  independent  variables  (optional) 

(Above 

values  may  be  in  any  order  and  indicate  the  order  of  the  data 

fields 

specified  on  the  format  card.  Values  must  be  packed,  left- 

justified;  nine  blanks  mean  "use  previous  variable  order".) 

n 

Plot  Y residual  vs  fitted  Y 

Blank 

Use  previous  plot  option  ,a 

0 

Do  not  plot. 

1 

Plot. 

12 

Plot  Y observed  vs  fitted  Y 

Blank 

Use  previous  plot  option. 

0 

Do  not  plot  . 

1 

Plot 

13 

Plot  Y observed  vs  XM 

Blank 

Use  previous  plot  option . a 

0 

Do  not  plot . 

1-7 

Plot  . 

8-9 

For  one -independent  variable  regressions,  plot  Y vs  XI  and 

regression  equation;  otherwise  plot  Y vs  Xl  only. 

14 

Zero/zero  plot  option 

Blank 

Use  previous  index  . 

0 

Plot  Y vs  XM  data  with  minimum  rectangular  grid. 

1 

Plot  Y vs  XM  data,  including  the  additional  point  (0,0). 

15 

Zero/negati ve  data 

acceptance  option 

Blank 

Use  previous  index. 

0 

Reject  all  nonpositive  data. 

1 

Accept  zero-value  data  points,  reject  negative  values. 

2 

Accept  all  data . 

16 

Format  for  Table  of  Residuals 

Blank 

Use  previous  index. a 

0 

Print  actual  values  rather  than  transformed  logarithms. 

1 

Print  Ln-transformed  variables  and  calculated  values 

for  Ln-linear  functions.  (Equation  index  > 6.) 

17 

Order  input  data  in  ascending 

dependent-variable  order 

Blank 

Use  previous  order  index. 

0 

Do  not  order  data  . 

1 

Order  data  . 

18 

Number  of  input  cards 

per  data-point  record 

Blank 

Use  previous  value.  Initial  setting  is  1 card/record 

(i.e.,  each  input  card  has  one  data  point). 

1-9 

Value  indicates  the  number  of  card  images  for  each  input 

record  (as  described  in  the  input  format  statement). 

19-20 

Iteration  limit 

Blank 

Use  previous  value.  Initial  setting  is  20. 

1-99 

Maximum  number  of  iterations  before  aborting  an  iterative 

solution  of  power  and  exponential  equations. 

21-30 

Specified  intercept 

Blank 

Normal,  unconstrained  regression 

Real 

Specified  intercept  (invalid  for  power  and  exponential 

equations) 

31-40 

Iteration  tolerance 

Blank 

Use  previous  value.  Initial  setting  is  10  . 

Real 

Solution-difference  tolerance  for  iterative  solutions 

41-80 

(d) 

a In i tl a 1 setting  is  zero. 
bSee  description  of  Col.  13. 

cZero  and  negative  values  cannot  be  used  safely  for  any  equation  except  linear  and  quadratic  because 
regressions  of  all  other  equations  require  logarithms. 

dSee  Appendix  C,  p.  58. 
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Equation 

Designator  Equation 


0 (Plot  only) 

1 Linear 

2 Quadratic 

3 Power 

4 Asymptotic-power 

5 Exponential 

6 Logarithmic-linear 

7 Semilogarithmic-linear  (Ln  Y vs  XI,  X2,  ....  X7) 

8 Semilogarithmic-linear  (Y  vs  Ln  XI,  Ln  X2,  Ln  X7) 


One  of  the  above  integers  should  be  entered  in  Col.  1 for  the  first  run 
(default  is  zero).  If  Col.  1 is  blank  after  the  first  run,  then  the 
value  for  the  previous  run  is  used.  Thus,  if  the  same  type  of  equation 
is  to  be  used  for  a series  of  runs,  its  designator  needs  to  be  entered 
only  for  the  first  run/  If  a zero  is  entered  as  an  equation  designator, 
no  regression  is  run;  instead,  the  input  data  are  plotted  as  Y versus  XM, 
where  M is  defined  in  Col.  13  as  one  of  the  independent  variables  being 
used  (if  Col.  13  is  left  blank,  a plot  against  XI  is  assumed).  It  is 
important  to  note  that  for  this  version  of  CURVES  a distinction  is  made 
on  the  Control  card  between  a blank  and  a zero  (0) . A blank  field  on 
the  Control  card  always  causes  the  program  to  retain  the  previous  value 
of  the  associated  indicator,  except  for  the  specified  intercept  (which 
reverts  to  the  condition  of  "no  specified  intercept") . 


Order  of  Variables 

Columns  2 through  10  indicate  the  order  of  the  variables  on  each 
data  card  as  described  in  the  variable-format  specification  statement. 
Depending  on  the  number  of  independent  variables  being  used  and  on 
whether  data-point  identifiers  are  being  used.  Cols.  4 to  10  may  be 


This  procedure  has  great  advantages  in  simplifying  the  input  oper- 
ations for  a job  consisting  of  more  than  one  run,  but  it  also  has  its 
drawbacks.  If,  for  example,  data  for  a run  are  removed  from  a block  of 
data  representing  a series  of  regressions  and  then  rerun,  control  in- 
formation that  should  be  passed  on  to  future  runs  may  be  lost  and  sub- 
sequent blanks  will  default  back  to  the  last  run  with  a nonblank  entry. 


Symbol  Type  of  Data 

I Identifier  (alphanumeric)  (optional) 

Y Dependent  variable  (required) 

1 First  independent  variable  (required) 

2, 3, 4, 5, 6, 7 Second  through  seventh  independent 
variables  (optional) 

The  independent  variables  are  identified  ordinarily  and  thus  require  that 
no  numerical  symbol  be  skipped;  i.e.,  use  of  3 implies  that  1 and  2 exist. 

Suppose  that  a set  of  data  is  to  be  entered  in  which  values  for 
three  independent  variables  are  located  in  Cols.  1-12,  13-24,  and  23-36, 
and  values  for  the  independent  variable  are  located  in  Cols.  37-48. 

Suppose  also  that  an  identifier  (a  six-digit  integer)  is  in  Cols.  55-60. 

Then  1,  2,  3,  Y,  and  I (123YI)  would  be  entered  in  Cols.  2-6  to  show  the 
above  order  for  a format  of  "(4F12.0,  6X,  A6).u 

Note  that  if  the  format  statement  is  written  with  tab  formats  (i.e., 
"T") , the  order  applies  to  the  order  of  fields  as  defined  in  the  format 
statement.  Thus,  for  the  above  example,  a format  of  "(T55,  A6,  Tl,  4F12.0)" 
would  require  an  order  of  1123Y. 

The  data-point  identifiers  may  be  up  to  eight  characters  long,  with 
a corresponding  format  specification  as  large  as  A8. 


) May  be  in  any  order 
( from  left  to  right. 


Plot  of  Residuals  vs  Fitted  Values 

Column  11  is  used  to  indicate  whether  a plot  of  the  residuals  versus 
the  fitted  values  is  desired.  If  a 1 (or  any  positive  digit)  is  entered 
in  Col.  11,  the  plot  is  generated;  if  a 0 (zero)  is  entered,  the  plot  is 
not  generated.  As  with  all  indices  on  the  control  card,  a blank  causes 
the  previous  value  to  be  retained. 

Plot  of  Observed  Y vs  Fitted  Values 

Column  12  is  used  similarly  to  Col.  11  to  generate  a plot  of  ob- 
served values  of  Y versus  their  fitted  values. 
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Plot  of  Observed  Y vs  Independent  Variable  or  Plot 
of  Regression  Equation 

Column  13  is  used  to  indicate  the  independent  variable  to  be 
plotted  with  the  observed  Y.  In  this  case,  one  of  the  digits  1 through 
M is  entered,  where  M is  the  index  of  the  independent  variable  to  be 
plotted.  In  addition,  for  the  case  of  a single  independent  variable, 
if  Col.  13  contains  an  8 or  9 instead  of  a 1,  the  resulting  regression 
equation  is  plotted  on  the  same  graph  using  dots  as  the  plot  symbols. 

In  all  cases  where  the  plotting  index  value  is  determined  to  be  invalid, 
the  plot  is  not  generated. 

For  all  plots  except  that  of  the  regression  equation  (dots),  26 
letters  and  nine  numerals,  or  a total  of  35,  are  used.  The  numeral 
zero  is  excluded  so  as  not  to  conflict  with  the  letter  "0".  For  plots 
of  data  with  more  than  35  points,  the  symbols  are  repeated  in  blocks  to 
account  for  all  data  points.  As  an  example,  for  37  data  points,  the 
points  in  order  from  1 through  37  would  be  represented  by  A,  A,  B,  B, 

C,  D,  E,  ...,  Z,  1,  2,  3,  ...,  9.  Each  symbol  used  in  a plot  is  also 
listed  beside  the  corresponding  data  point  in  the  Table  of  Residuals 
so  that  the  data  point  can  be  readily  identified. 

Plot-size  Option 

Column  14  is  used  as  a designator  for  a plot-size  option  for  the  Y 
versus  X plot.  Normal  operation  of  all  plotting  is  to  determine  the 
minimum  and  maximum  ordinate  and  abscissa,  and  then  produce  the  minimum 
size  (largest  scale)  rectangular  plot  containing  all  the  data.  In  the 
case  of  the  Y vs  X plot,  especially  for  a positive  Y-intercept,  the 
point  (0,0)  can  be  added  to  the  plot-points  (to  depict  the  first  quad- 
rant) by  entering  a 1 in  Col.  14. 

Data  Value  Acceptance 

Column  15  is  used  as  an  indicator  for  data  value  acceptance.  Nor- 
mally, a negative  value  for  a data  point  terminates  the  job  because 
negative  values  cannot  be  used  in  some  cases  for  the  logarithmic  equa- 
tions and  in  all  cases  involving  the  power,  asymptotic  power,  and  ex- 
ponential equations.  A zero  value  causes  the  program  to  reject  the  data 
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point  (except  the  origin,  which  has  a special  meaning  as  described 
below)  but  to  continue  the  processing  (reading).  The  latter  condition 
allows  a data  set  to  be  processed  even  though  there  are  some  missing 
data.  (A  blank  is  read  as  zero  in  a floating-point  format.) 

In  some  cases,  it  is  desirable  to  accept  zero  as  a valid  data 
value  (e.g.,  a linear  regression  with  a zero/one  dummy  variable);  to 
do  so,  a 1 is  entered  in  Col.  15.  However,  note  that  a regression  of 
any  equation  except  the  linear  and  quadratic  requires  logarithms  of 
some  or  all  of  the  input  variables.  If  the  program  attempts  to  take 
the  logarithm  of  a zero  or  negative  value,  the  job  terminates.  If  Col. 
15  contains  a 2,  all  data  values  (positive,  zero,  and  negative)  are 
accepted  except  the  origin,  which  is  interpreted  to  signify  that  the 
reading  of  data  is  complete. 


Output  Status  for  Logarithmic  Eauations 


Column  16  is  used  to  indicate  the  output  status  for  the  logarithmic- 
equations — equation  designators  6,  7,  and  8.  For  each  of  these  equa- 
tions, the  regression  is  performed  on  the  logarithms  of  some  or  all  of 
the  variables.  If  the  column  is  left  blank,  all  logarithmic  results 
for  the  Table  of  Residuals  are  printed  as  nonlogarithmic  (antilog)  data 
for  better  readability  and  understandability . That  is,  the  original 
input  data  (before  logarithmic  transformation)  are  printed,  the  fitted  Y 
values  are  exponentiated  if  in  Ln  Y form  (equation  designator  value  of 
6 or  7),  and  the  residuals  and  relative  Y deviations  are  calculated 
based  on  nonlogarithmic  data.  To  use  a logarithmic  model  with  all  ap- 
propriate output  in  the  Table  of  Residuals  in  logarithmic  data,  a 1 is 
entered  in  Col.  16.  Regardless  of  the  value  entered  in  Col.  16,  the 
Summary  Table  always  produces  statistics  for  the  regressed  equation 
whether  or  not  the  equation  is  logarithmic. 


Ordering  of  Input  Data 


Column  17  is  used  to  designate  whether  the  input  data  are  to  be 
ordered  from  low  to  high  values  of  Y.  A value  of  1 signifies  that  the 
data  are  to  be  ordered.  For  the  first  run  a blank  (or  zero)  signifies 
that  the  data  are  not  to  be  ordered;  for  subsequent  runs  a blank  signi- 
fies that  the  value  of  the  order  designator  for  the  preceding  run  is  to 
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be  used.  Again,  this  is  done  so  that  if  all  runs  in  a series  are  to 
be  either  ordered  or  unordered,  the  order  designator  need  only  to  be 
entered  for  the  first  run.  To  reset  the  order  indicator  to  zero  (no 
ordering),  simply  enter  a 0 (zero). 

Figure  3 also  shows  an  example  of  a control  card  in  which  a linear 
regression  is  to  be  made  on  the  input  data  (1  in  Col.  1).  The  data  are 
to  be  ordered  with  respect  to  Y (1  in  Col.  17). 

Number  of  Input  Cards  per  Data-Point  Record 

Column  18  is  used  to  indicate  the  number  of  input  cards  per  data- 
point  record.  If  the  input  data  are  entered  from  the  card  reader,  the 
entry  in  Col.  18  is  the  number  of  physical  (card)  records  per  logical 
record  (data  point);  for  disk  or  tape  input  data,  see  Format  Card  in 
the  next  subsection  under  SPECIFICATION  CARDS.  The  initial  setting  is 
one  card  per  record;  if  a value  of  zero  is  entered,  a value  of  1 is 
substituted.  The  maximum  number  of  cards  per  record  is  nine.  This 
designator  is  used  primarily  to  block  the  input  data  into  individual 
records  when  the  READ  MEMORY  option  is  used;  a slash  is  not  permitted 
in  the  variable  format  statement  when  using  memory  because  it  repeats 
the  same  record  instead  of  advancing.  This  value  has  no  use  for  the 
READ  ONCE  or  READ8  ONCE  options.  The  initial  setting  is  one  card  per 
record.  A blank  value  repeats  the  previous  value. 

Maximum  Number  of  Iterations  Allowed 

Columns  19  and  20  are  used  to  specify  the  maximum  number  of  iter- 
ations to  be  allowed  before  aborting  the  iterative  solution  of  the 
power  and  exponential  equations.  The  default  setting  is  100. 

Y-Intercept  Specification 

Columns  21  through  30  are  used  to  specify  the  equation  intercept 
term  (regression  constant).  This  is  the  A value  for  all  equations  ex- 
cept the  Ln-linear,  for  which  it  is  Ln  A.  The  intercept  for  the  power 
and  exponential  equations  may  not  be  specified.  The  desired  intercept 
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is  entered  in  floating-point  format  anywhere  within  the  field.  The 
implied  decimal  point  location  is  at  the  right  end  of  the  field,  after 
Col.  30.  If  the  field  is  blank,  the  unspecified  regression  equation 
is  assumed. 

Iteration  Tolerance 

Columns  31  through  40  are  used  to  specify  the  iteration  tolerance 
(DELTA).  For  regressions  of  the  nonlinear  equations,  an  interative 
solution  is  achieved  when  the  ratio  of  the  value  of  each  parameter  to 
its  value  for  the  previous  iteration  differs  from  unity  by  an  amount 
equal  to  or  less  than  DELTA.  DELTA  is  initialized  at  10  . The  user 

may  change  this  by  specifying  another  constant  in  Cols.  31-40,  prefer- 
ably in  decimal  format,  for  example,  0.00000001.  The  implied  decimal 

-12 

point  is  prior  to  Col.  31.  Any  value  of  DELTA  less  than  10  or 
greater  than  10  ^ is  reset  to  10  ^ . The  authors  recommend  that  DELTA 
not  be  changed,  eccept  for  very  special  reasons;  a DELTA  of  10  ^ should 
be  sufficient  to  obtain  all  parameters  within  the  accuracy  printed. 

SPECIFICATION  CARDS 

There  are  basically  four  specification  cards  that  provide  control 
information  for  the  run  similar  to  the  control  card  but  that  requires 
more  space  to  express.  They  are:  Format,  Label,  Guess,  and  Read.  They 

follow  the  control  card  and  may  be  in  any  order  except  that  the  Read 
card  must  be  the  last  one.  For  any  such  card,  the  specification  name 
is  entered  beginning  in  Col.  1;  the  information  following  the  name 
should  begin  in  Col.  9. 

Format  Card 

The  Format  card  indicates  where  the  data  are  located  on  the  data 
cards.  This  card  must  begin  with  the  word  FORMAT  in  Cols.  1-6  followed 
by  a left  parenthesis.  A matching  right  parenthesis  closes  the  format 
specification.  Information  within  the  parentheses  must  conform  to  the 
rules  for  FORTRAN  formats.  In  addition,  except  for  the  alphanumeric 
identifiers,  all  input  data  must  be  in  real-number  (floating-point) 
formats.  One  or  two  continuation  cards  may  be  used,  each  of  which  must 
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likewise  contain  the  word  FORMAT  in  Cols.  1-6.  After  the  final  right 
parenthesis  of  the  format,  the  user  may  make  any  comments  desired 
(e.g.,  identifying  the  variables  the  format  refers  to),  because  data 
beyond  the  last  parenthesis  are  ignored  by  the  program.  If  a Format 
card  is  omitted  after  the  first  run,  format  specifications  are  carried 
over  from  the  previous  run.  The  Format  card  shown  in  Fig.  3 could  be 
used  for  the  previous  example,  as  discussed  before. 

Label  Card 

The  Label  card  is  used  to  identify  the  variable  with  an  eight- 
character  label.  The  labels  are  placed  next  to  the  regression  parameter 
or  variable  to  which  they  refer  in  the  output  Summary  Table;  they  also 
appear  as  table  headings  in  the  Table  of  Residuals.  If  an  alternate 
form  (LABELS)  is  used,  the  Summary  Table  also  contains  a printout  of 
the  names  of  the  regression  variables  immediately  below  the  title. 

(This  alternative  may  allow  the  user  to  use  the  same  title  with  several 
runs,  but  to  distinguish  the  runs  by  the  variable  list;  it  may  also 
reduce  the  requirements  of  Title  card  preparation.) 

If  no  Label  (or  Labels)  card  is  used,  the  default  (for  each  run) 
is  listed  below.  Labels  do  not  carry  over  from  one  run  to  the  next; 
they  must  be  entered  each  time.  Table  5 shows  the  format  of  the  Label 
card . 

Guess  Card 

The  Guess  card  is  used  to  provide  initial  guesses  of  the  parameter 
values  for  the  iterative  solutions  in  fitting  the  power  and  exponential 
equations.  The  normal  solution  procedure  in  the  program  is  first  to 
make  a least-squares  regression  of  the  corresponding  logarithmic-linear 
(or  semilogarithmic-linear)  equation  to  obtain  preliminary  estimates  of 
the  parameters  (the  values  are  printed  out  under  the  heading  INITIAL 
GUESS  in  the  Summary  Table)  and  then  iterate.  The  user  may  make  the 
initial  guesses  by  entering  a Guess  specification  card  in  the  format 
shown  in  Table  6.  The  implied  decimal  point  is  at  the  right  end  of 
each  field.  The  authors  recommend  that}  except  for  special  circum- 
stances, the  user  should  not  attempt  initial  guesses  for  the  iterative 
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Card  Columns 


Table  5 


FORMAT  OF  LABEL  SPECIFICATION  CARD 


Value 


LABELS 


Specification  Description 
Indicates  the  type  of  specification  card 


LABELS 


17-24 


25-32 


33-40 


Not  used 


Alphanumeric  Heading  for  identifier,  if  used  (default 
is  LABELS!*);  if  no  identifier,  Xl*Kl*Ml*K. 

Alphanumeric  Heading  for  dependent  variable  (default 
is  WWWW). 

Alphanumeric  Heading  for  first  independent  variable 
(default  is  XXXXlXXl*) . 

Alphanumeric  Heading  for  second  independent  variable, 

if  used  (default  is  1*1*!*X2!*1*1*  for  all  func- 
tions except  quadratic  and  asymptotic; 
quadratic  default  is  XXX1**21*;  asymptotic 
default  is  EXPONENT). 


73-80  Alphanumeric  Heading  for  seventh  independent  variable, 

if  used  (default  is  Xl*l*X71*l*X)  . 

3A  blank  column  is  indicated  by  ]*.  (See  pp.  21-22  for  definition 
of  variables.) 


Table  6 


FORMAT  OF  GUESS  SPECIFICATION  CARD 


Card  Columns 


17-24 

25-32 


Value 


GUESS!* 


Specification  Description 

Indicates  the  type  of  specification  card 
Not  used 

Initial  guess  for  A 
Initial  guess  for  B 
Initial  guess  for  C,  if  used 


65-72 


Initial  guess  for  H,  if  used 


See  p.  5 for  description  of  A,  B,  C,  ...  H use. 
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oases.  It  is  quite  likely  that  by  using  the  Guess  card,  it  will  take 
the  program  much  longer  to  converge  to  the  solution,  if  indeed  it  con- 
verges at  all. 


Read  Card 

The  final  specification  card  is  the  Read  card.  This  card  (which 
is  a required  input)  signifies  that  the  data  are  ready  to  be  processed. 
(An  alternate  reading  file  for  the  data  cards  is  discussed  below  under 
the  READ8  option.)  The  format  of  the  Read  card  is  shown  in  Table  7. 

Anything  other  than  ONCE,  MEMO,  or  DISK  in  Cols.  9-12  is  equiva- 
lent to  and  will  result  in  the  data  being  reread  from  the  save-data 

file  last  used  (disk  or  memory).  If  the  last  option  was  ONCE,  a new 
set  of  cards  will  be  read  (because  no  data  were  stored).  The  initial 
setting  is  ONCE.  Since  rereading  of  saved  data  requires  no  more  data 
input,  data  cards  are  not  expected  nor  may  they  be  present  in  the  input 


Table  7 

FORMAT  OF  READ  SPECIFICATION  CARD 


Card  Columns 

Value 

Specification  Description 

1-6 

READM 

or 

READ  8)5 

Indicates  the  type  of  specification  card;  data 
are  read  in  from  card  unit  5 (i.e.,  sequentially 
from  the  normal  card  reader)  or  from  user-defined 
unit  8. 

00 

1 

Not  used 

9-12 

ONCE 

The  save-data  option  parameter: 

Do  not  save  the  data,  just  process  it  for  this 
one  run. 

MEMO 

Save  the  card  images  in  memory  (^  100  cards)  for 
rapid  access  on  later  runs. 

DISK 

Save  the  card  images  on  scratch  disk  unit  4 (up 
to  the  maximum  number  of  regression  data  points) 
for  access  on  subsequent  runs. 

nu 

Use  saved  data 

13-80 

Not  used 
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stream  when  this  option  is  in  effect.  Figure  3 shows  an  example  of 
the  READ  ONCE  option. 

In  the  event  that  a data  input  stream  other  than  the  normal  se- 
quential card  input  (i.e.,  unit  5)  is  desired,  the  user  may  substitute 
READS#  for  READ##  in  Cols.  1-6.  This  alternative  sets  up  unit  8 for 
the  input  of  data  cards  only  (including  the  BLANK  card).  The  Title, 
Control,  and  Specifications  cards  always  come  from  unit  5. 

Any  data  saved  must  be  in  card  images.  When  using  memory,  data 
are  stored  in  a record  length  (in  bytes)  of  80  times  the  number  of 
cards  per  record  (maximum  of  720  bytes — 9 cards — in  the  program) . 
Therefore,  in  reading  back  such  stored  data  from  memory,  the  same 
record  length  must  be  used  in  the  format  statement;  e.g.,  FORMAT  (T200, 
3F8.0).^  A slash  mark  is  not  permitted  in  this  case,  because  it  repeats 
the  same  80-byte  record  in"tead  of  advancing  to  a new  record.  For  a 
scratch  disk,  data  are  stored  in  80-byte  records  and  hence  must  be  read 
back  the  same  way.  Therefore,  slash  marks  must  be  used  in  the  format 
statement  for  multiple  cards  per  record  when  using  disk. 

DATA  CARDS 

Each  data  point  must  contain  at  least  a pair  of  values,  one  for 
the  dependent  variable  (Y)  and  one  for  the  independent  variable  (XI). 

Each  set  of  data  constituting  a run  must  contain  at  least  as  many  data 
points  as  the  number  of  parameters  being  estimated  and  may  contain  up 
to  the  maximum  allowable  as  indicated  in  Table  1.  If  the  maximum  is 
exceeded,  an  error  message  explains  this  fact.  The  location  of  the 
data  on  the  card  must  be  in  exact  agreement  with  the  information  entered 
on  the  Control  and  Format  cards,  or  else  the  data  will  not  be  read 
properly.  The  numerical  data  (dependent  and  independent  variable  values) 
are  read  as  real  (floating-point)  numbers  and  data-point  identifiers  (if 
used)  as  alphanumeric  data. 

^The  largest  repetitive  skip  that  can  be  used  in  a FORTRAN  format 
statement  is  255  (i.e.,  255X  or  T255).  If  a larger  skip  is  required, 
simply  break  it  up  into  two  or  more  skips.  For  example,  a skip  of  300 
columns  may  be  entered  in  the  format  statement  as  150X,  150X;  or  100X, 
200X;  etc. 


h- 
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For  any  data  card,  if  either  the  Y field  or  any  of  the  X fields, 
but  not  all,  is  blank  or  contains  the  value  zero,  that  card  is  usually 
skipped.^  However,  if  all  X and  Y fields  are  blank  (zeros),  and  the 
identifier  field,  if  used,  is  blank  or  contains  the  word  BLANK,  the 
reading  of  input  data  for  the  run  is  terminated  at  that  point  (see 
BLANK  CARD  below). 

If  the  identifier  field  is  neither  blank  nor  contains  the  word 
BLANK  and  all  the  data  are  zero,  the  data  point  is  merely  rejected. 
Thus,  the  user  may  identify  places  where  data  are  needed  in  the  data 
set,  but  for  which  data  are  not  yet  available. 

Figure  3 also  shows  a data  card  containing  data  in  the  specific 
format . 

BLANK  CARD 

Each  set  of  data  cards  constituting  a run  must  always  end  with  a 
blank  card.  This  card  is  used  to  terminate  the  reading  of  the  input 
data  for  a given  run.  There  must  be  as  many  blank  cards  as  there  are 
cards  per  record  (see  description  under  CONTROL  CARD,  Column  18) . As 
an  option,  the  user  may,  if  saving  data  in  memory  or  on  disk,  type  the 
word  BLANK  in  Cols.  1-5  of  a single  blank  card. 

DONE  CARD 

After  each  set  of  data  is  read  and  processed,  the  machine  cycles 
back  to  read  another  set  of  data.  To  terminate  a run  or  a series  of 
runs,  the  word  DONE  is  entered  in  Cols.  1-4  on  what  would  be  the  Title 
card  for  the  next  run.  This  causes  the  program  to  print  a large  DONE 
on  a separate  page  and  stop. 

SUMMARY 

Figure  4 briefly  summarizes  the  information  on  the  Title,  Control, 
Specification,  Data,  Blank,  and  Done  cards.  Note  that  the  Read  card 

^As  was  stated  previously  under  Data  Value  Acceptance,  an  option 
is  provided  to  allow  zero  or  negative  values  to  be  accepted  by  the  pro- 
gram. Nevertheless,  the  default  condition  is  to  skip  zero  fields  and 
to  terminate  the  run  on  reading  negative  values, 
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must  always  be  the  last  Specification  card- for  any  run.  Figure  5 shows 
the  order  of  two  data  decks  for  a series  of  three  runs.  In  the  figure, 
the  second  deck  also  contains  data  for  the  third  run.  Hence,  it  is  to 
be  saved  for  rereading  during  the  third  run. 


Third  run 

( Data  for  ibis  run  are  to 
be  read  back  from  memory 
or  off  disk.) 


Second  run 
(Data  for  this  run 
are  to  be  read  in 
the  same  format 
as  the  first  run.  ' 
In  addition,  the 
cards  are  to  be 
written  into  memory 
or  onto  disk  since 
they  contain  data 
for  the  third  run.) 


f Done  card 
Read  card 


/ Format  cord 
/ Control  card 
Title  card 


Blank  card 


( Data  deck 
^ Read  card 
/ Control  card 
f Title  card 
Blank  card 


( Data  deck 
f Read  card 
Format  dord 
Control  card 
Title  card 


First  run 


Fig.  5 — Arrangement  of  two  data  card  packs  for  three  runs 
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IV.  PROGRAM  OUTPUT 


Figure  6 is  an  example  of  program  output  for  fitting  linear, 
power,  logarithmic-linear,  and  quadratic  equations.  The  output  also 
includes  examples  of  the  plots  available  in  the  program.  For  compar- 
ability and  continuity,  the  input  data  for  the  first  two  runs  are  the 
same  as  in  the  example  in  the  previous  CURVES  report  (RM-5762-PR) . 

The  particular  data  shown  have  no  real  meaning  and  are  presented  only 
for  the  purpose  of  displaying  the  output.  Because  the  headings  are 
self-explanatory  and  the  statistics  have  already  been  discussed,  no 
explanation  of  them  is  given  here. 

Figure  7 is  a listing  of  the  input  data  for  the  four  runs  whose 
output  is  shown  in  Fig.  6.  The  data  are  in  the  same  format  on  the 
cards  as  shown  in  Fig.  3.  As  Fig.  7 shows,  the  input  data  for  the 
first  run  (linear  regression)  are  read  in  a P-format  in  order  to  scale 
the  data.  (The  scale  option  in  the  previous  version  of  CURVES  has 
been  deleted  because  scaling  can  be  accomplished  with  the  use  of  the 
P-format.)  Note  that  for  the  third  run,  the  word  BLANK  in  Cols.  1-5 
is  used  rather  than  a blank  card  to  terminate  reading  of  the  data 
cards.  Either  method  is  acceptable  to  CURVES  for  saved  data. 

Table  8 lists  various  central  processing  unit  (CPU)  execution 

times  for  typical  runs  made  on  the  Rand  IBM  370/158  computer.  (The 

sample  output  shown  in  Fig.  6 required  2.5  seconds  of  execution  time 

* 

and  92,000  bytes  of  core  storage,  based  on  an  H-compiled  program.) 


If  data  are  stored  on  a scratch  disk  or  tape,  instead  of  in 
memory,  additional  core  storage  may  be  required. 
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rg  rg  rg  rg  rg 


•AINIA'JP  RELATIVE  DEVIATION  - -0.18814,  *AFAN  ABSOLUTE  RELATIVE  DEVIATION  « 0.03995,  VAX  I HUH  RELATIVE  DEVIATION  » 0.08T55 
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TEST  RUN  1 -- 

LINEAR  WITH 

THREE  INDEPENDENT  VARIABLES 

Title  Card 

1 123YI 

1 

Control  Card 

FORMAT  (1PF12 

.0,  0P3F12. 

0,  6X , AS) 

XI  SCALED 

BY  10**  (- 

1) 

Format  Card 

R FA D ONCE 

Read  Card 

20. 19283 

10. 34567 

4.27653 

27.94618 

100404 

> 

50.37289 

10.  0 

20.0 

80.09382 

10131  3 

105.  35467 

55. 27618 

12.16547 

100.47238 

102121 

62. 17862 

18. 88886 

29.17654 

89.92674 

101717 

35.25671 

42. 16543 

27.17864 

33.60992 

1 0050  5 

80. 24561 

35. 2541  1 

15.25672 

83.40204 

10141  4 

12.15437 

3.26751 

31.26884 

40.71304 

100707 

65. 33821 

16. 27865 

1.27553 

84.84CC2 

101515 

89.26718 

109.26547 

44.  27861 

54.53670 

101010 

100.0 

100.  C 

100. C 

107.75168 

102323 

121.  17625 

175. 26876 

52.17625 

40.97917 

100808 

2. 37658 

15.28765 

30.98716 

19.44080 

100202 

106.  26789 

81. 27543 

107.26784 

137.0 

102727 

56. 27182 

19.2671 3 

41.26517 

88.26627 

101616 

15.26718 

40. 15672 

63.17772 

24.06789 

100303 

\ Data  Cards 

79. 26182 

12.25418 

90.2671 3 

147.55076 

103030 

51.22268 

86.12357 

24.16273 

14.39608 

1C0101 

97.26173 

84.23456 

93.25671 

111.35247 

102424 

62. 24518 

10.27625 

17.24561 

92.51133 

101919 

32. 27615 

8.  18761 

16.27615 

51.71629 

100909 

38. 27b1 5 

27.37677 

3.28716 

40.40568 

100606 

49. 28817 

14. 11167 

54.28817 

94.38172 

102020 

94.24689 

48. 23418 

35.28861 

104.27262 

102222 

40.  19824 

15.0 

39.22218 

71.09687 

101212 

117. 33922 

95. 18293 

100.0 

131.24108 

102626 

75.12345 

76.11111 

88.99112 

90.28862 

101818 

10.0 

20.0 

101.0 

59.  C 

10111  1 

74. 28192 

18. 22184 

97.22215 

138.46147 

102828 

89.22186 

12. 16524 

8.01 187 

122.69347 

102525 

j 

64.  16253 

102.33728 

114.23477 

143.72068 

102929 

/ 

Blank  Card 

TEST  RUN  2 — 

POWER  WITH 

THREE  INDEPENDENT  VARIABLES 

Title  Card 

3123YI 

Control  Card 

FORMAT  (4F12. 

0,  6 X , A 6 ) 

Format  Card 

R EAD 

Read  Card 

35.  18762 

10.26781 

16.25673 

351.20949 

200606 

\ 

102. 15672 

54.28716 

40.25617 

1510.21359 

202929 

78.28719 

66.26718 

100.0 

661.12500 

20151  5 

102. 17628 

112. 27162 

130.18719 

880.20145 

202020 

70.27168 

108. 26152 

132. 17817 

535.55565 

201313 

43. 13425 

59.  17236 

18.21926 

800.  C 

201818 

10.21822 

65.27819 

95.11118 

60.25169 

200101 

67.22218 

48.18273 

92.23116 

530.19078 

201 11  1 

100. C 

100.0 

100.0 

980.76950 

2C2222 

> Data  Cards 

120.11234 

77. 19283 

105.25172 

1138. 31287 

202525 

34. 18273 

26.18273 

5.26174 

972.57176 

202121 

29.  12678 

8. 16279 

75.25619 

115.19201 

200202 

20.23457 

14.16289 

4.25617 

46C. 16664 

200909 

57.  19287 

102.0 

17.  26153 

1345.54035 

202727 

120. 17865 

90.  24133 

105.25671 

1158. 45163 

202626 

25.14567 

34.25671 

90.0 

1 43.CU087 

200303 

17. 26155 

76. 25144 

10.15782 

370.88634 

2CC7C7 

65.24138 

110.27816 

21.25517 

1472.22902 

202828 

28. 12816 

9. 18826 

3.27715 

779.6280C 

201717 

41.  15287 

118.  26132 

77.23518 

425.47155 

200808 

i 

35.1231 1 

19.23518 

37.26153 

342.15482 

2C0505 

Fig.  7 — Deck  setup  for  sample  runs  shown  in  Fig.  6 
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121.^7157 

55.44361 

1 10. 23145 

10CC. 58691 

2C2323 

23.15o78 

104.28715 

19.49175 

461.97847 

201C10 

111.25411 

23.18726 

76.1623  3 

81 4. 8790C 

2C1  919 

1d1 .28761 

93.27615 

178.49977 

1125.67637 

202424 

44. 2765  1 

85.46716 

27.19279 

755.  794d  7 

201 61 6 

71.16253 

ICO.  16273 

145. 27168 

531.15436 

2C1212 

112.18862 

147.23557 

25.18892 

-65C.01C76 

2C3C30 

38.29418 

27.16827 

75.47164 

247,96037 

2 0C  409 

89.26615 

5.27715 

45. 2351 9 

585. 28353 

201414 

Data  Cards 
(Continued) 


Blank  Card 

TEST  RUN 

3 

— LN-LINEaR  kITH  SEVEN  IFEEFEEDENT 

VARIABLES  (MAXIMUM) 

Title  Card 

6IY 1234567 

c 

Control  Card 

FOR  HAT 

(A7 

$ 3 X , 

6F5. 

0) 

Format  Card 

READ 

MEMORY 

T 

X • 

E.  , 

READ  CARDS, 

SAVE 

in  n 

EMCRY 

Read  Card 

AIFHa 

1.0 

2.1 

1.1 

6.  6 

. 94 

77 

76 

14  \ 

BETA 

2.0 

4.1 

2.9 

8*9 

. 81 

86 

12 

92  ' 

SA8.MA 

3.1 

9.0 

3.7 

3.2 

. 61 

98 

St 

64 

DELTA 

4.0 

8.0 

3.5 

3.7 

. 63 

45 

66 

75 

E FS ILCN 

7.0 

13.5 

4.1 

3 < 9 

.22 

71 

6 4 

91 

2 ETA 

5.0 

9.9 

3.8 

7.4 

. 99 

81 

42 

43 

El  A 

9.9 

1 9.  0 

3.7 

7.6 

. 36 

84 

93 

76 

> Data  Cards 

THETA 

10. C 

21.0 

6.9 

2.6 

.40 

59 

29 

80 

IOTA 

11.0 

24.0 

8.5 

1 .8 

.12 

6 5 

69 

33 

KAPPA 

12.0 

25.0 

7.3 

2.3 

. 05 

4 1 

34 

37 

LA68DA 

15.0 

19.0 

3.4 

0.6 

. 35 

14 

29 

09 

MU 

17.0 

36.  5 

8.4 

5.6 

.18 

02 

73 

83 

OMEGA 

20.0 

40.0 

7.6 

8.7 

. 4 4 

46 

25 

21  / 

BLANK 

Blank  Card 

TEST  RUN 

4 

— QUADRAT 

IC  k 

ITH  AIL 

PLOTS,  SAME  Y 

AND  X2  (LION  XI)  DATA  AS  RUN  3 

2IY1 

1 191 

Control  Card 

FOE.  NAT 

(Ab 

, 2X, 

F5.0 

. SX 

, F5.0) 

Format  Card 

LABEL 

HEADING 

T CTLCC JTPCPULA  T NFOP** 

2 

Label  Card 

READ 

Read  Card 

DONE 

Done  Card 

\ ' 

Fig.  7 (cont.) 
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Table  8 

APPROXIMATE  CPU  TIMES  (SECONDS)  FOR  VARIOUS 


CURVES  RUNS  ON 

IBM  370/158 

Number 

of  Data 

Points 

Number 

Equation  Index  of  Variables  20 

60 

100 

( 1 

.3 

.6 

1.0 

1,2,6, 7, 8a  { 3 

.4 

.8 

1.3 

1 7 

.6 

1.2 

1.6 

i 1 

.6 

1.7 

2.3 

3b  ) 3 

1.6 

5.5 

7.0 

\ 7 

4.6 

12.7 

17.3 

4C  1 

1.9 

4.8 

7.7 

( 1 

.5 

1.0 

1.3 

5b  \ 3 

.7 

1.4 

2.2 

\ ^ 

1.3 

2.5 

3.4 

Xinear,  quadratic,  and  logarithmic  equations, 
each  having  an  algebraic  solution  for  parameter 
estimates. 

kpower  and  exponential  equations,  each  using  modi- 
fied Gauss-Newton  iterative  procedure  to  obtain 
parameter  estimates. 

CAsymptotic-power  equations  using  iterative, 
incremental-stepping  procedures  to  obtain  parameter 
estimates . 
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Appendix  A 

NONLINEAR-LEAST-SQUARES  CONSIDERATIONS 


LOGARITHMIC  AND  NONLOGARITHMIC  EQUATIONS 

The  usual  procedure  for  deriving  least-squares  estimates  of  the 
parameters  of  the  power  or  exponential  equation  is  first  to  convert 
the  equation  into  a logarithmic-linear  (or  semilogarithmic-linear) 
equation.  One  then  has  an  equation  for  which  least-squares  estimates 
of  the  parameters  can  be  obtained  by  simple  algebraic  means.  However, 
note  that  these  least-squares  estimates  are  not  the  same  as  the  least- 
squares  estimates  of  the  parameters  that  specify  the  original  equations. 
This  may  be  seen  by  considering,  for  example,  the  power  equation  and 
its  logarithmic  form. 

Let 


Y = A • X1B  • X2C  • X3°  • ...  • X7H 


and 

Ln  Y = Ln  A + B • Ln  XI  + C • Ln  X2  + D • Ln  X3  + . . . + H • Ln  X7  . 

For  a least-squares  solution,  one  is  interested  in  minimizing  the 
sum  of  squares  of  the  Y residuals  (denoted  by  Q) . ^ Therefore,  for  the 
power  equation, 


o ■ l «i  - Ycl>' 

i=l  1 1 


and  for  the  logarithmic  equation 


Q'  = l (Ln  Y - Ln  Y )‘ 
i*=l  1 C1 


Throughout  this  discussion,  Q is  used  to  represent  the  sum  of 
squares  of  the  Y residuals. 


i 

t; 
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or 


where  N = number  of  data  points, 

= observed  value  of  dependent  variable  for  ith  data  point, 

Y = fitted  value  of  dependent  variable  for  ith  data  point. 

In  the  logarithmic  case,  the  sum  of  squares  of  the  actual  differ- 
ences (residuals)  between  the  observed  and  fitted  Y values  is  not 
being  minimized,  rather  the  sum  of  squares  of  the  logarithms  of  the 
ratios  of  those  values  is  being  minimized.  Depending  on  the  observa- 
tions, the  two  procedures  may  produce  substantially  different  estimates 
of  the  parameters  A,  B,  C,  ...»  H. 

It  may  also  be  seen  that  any  statistic  based  on  the  sum  of  squares 
of  Y residuals,  such  as  the  coefficient  of  determination,  may  be  mis- 
leading if  used  to  compare  the  logarithmic  form  with  its  nonlogarithmic 
counterpart.  For  the  logarithmic  form,  such  statistics  are  based  on 
logarithms  and  hence  have  different  meanings. 

Regression  theory  states  that  if  the  error  term  on  the  dependent 
variable  Y is  an  additive,  normally  distributed  random  variable  with 
a mean  of  zero,  then  a least-squares  fit  will  lead  to  maximum  likeli- 
hood estimates  of  the  regression  coefficients.  Therefore,  for  the 
power  form,  the  following  is  assumed: 

Y = A • X1B  • X2C  • X3D  • ...  • X7H  + £ 

where  the  errors  are  independent,  normally  distributed  random  variables 
with  mean  zero  and  a common  variance.  For  the  logarithmic  form,  one 
has 

Ln  Y = Ln  A + B Ln  XI  + C Ln  X2  + D Ln  X3  + ...  + H Ln  X7  + Ln  6 


or 
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Y = A • X1B  • X2C  • X3°  • . . . , X7H  • 6 , 

where  the  error  terms  Ln  6 satisfy  the  conditions  specified  for  the 
errors  in  the  previous  case.  In  this  case,  the  error  term  is  multi- 
plicative. 

The  question  as  to  whether  the  regressed  power  equation  or  its 
regressed  logarithmic  form  is  more  appropriate  for  a set  of  data  depends 
on  many  factors  including  the  error  term  associated  with  the  data  and 
what  criterion  is  used  for  a "good  fit."'*'  However,  one  of  the  best 
tests  for  comparison  is  to  examine  the  plot  of  Y residuals  versus  fitted 
values  (or  residuals  of  Log  Y versus  fitted  Log  Y for  logarithmic  case). 
The  "better"  model  (power  versus  log)  will  show  a more  random  normal 
distribution  of  the  Y residuals  around  the  zero  line.  This  plot  is 
available  In  the  program  for  such  an  examination. 

NONLINEAR  SOLUTIONS 

It  is  a necessary  condition  that  the  first  partial  derivatives  of 
Q with  respect  to  the  parameters  must  be  zero  in  order  that  Q be  min- 
imized. This  is  not,  unfortunately,  a sufficient  condition  for  a func- 
tion that  is  not  linear  with  respect  to  all  of  its  parameters.  The 
reason  for  this  is  that  if  Q could  be  graphed  (in  multi-dimensional 
space)  for  a nonlinear  function,  there  might  be  other  critical  points — 
such  as  saddle  points  or  relative  maxima  or  minima  points — where  the 
first  partial  derivatives  would  also  be  zero.  A test  that  checks  for 
this  possibility  involves  examining  the  matrix  of  second  partial  deriva- 
tives of  Q,  which  is  a generalization  of  the  second-derivative  test  for 
a one-parameter  case.  If  this  matrix  is  positive-definite  for  all  param- 
eters in  a region  containing  a solution,  it  can  be  shown  that  the  solu- 
tion represents  an  absolute  minimum  for  Q in  that  region  and  is  the  only 

2 

solution  in  that  region.  However,  if  the  matrix  is  not  positive-definite 

*Tor  the  interested  reader,  this  question  is  treated  in  Graver  and 
Boren,  RM-4879-PR. 

2 

H.  0.  Hartley,  "The  Modified  Gauss-Newton  Method  for  the  Fitting 
of  Non-Linear  Regression  Functions  by  Least-Squares,"  Technometrics , 

Vol.  3,  No.  2,  May  1961,  pp.  273-274. 
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at  all  points  in  that  region,  then  there  may  be  other  "solutions"  for 
the  same  set  of  data. 

For  regressions  of  the  power  and  exponential  equation  involving 
very  large  (s  10^)  or  very  small  (s:  10  values  of  X or  Y input  data, 

the  matrix  of  partial  derivatives  may  not  be  inverted  accurately  enough 
to  give  reasonable  corrections  to  the  parameters.  Consequently,  in 
such  cases  there  may  be  no  convergence  to  a solution  of  the  parameters. 
This  problem  can  usually  be  remedied  by  rescaling  the  input  data  (using 
the  P-format) . If  input  data  are  rescaled,  only  the  estimate  of  the 
parameter  A is  changed  in  the  regression. 

In  summary,  one  should  be  aware  that  for  a nonlinear  equation  as 
defined  in  this  report,  the  "solution"  obtained  may  not  represent  an 
absolute  minimum  for  Q.  The  only  sure  way  to  know  is  to  try  all  com- 
binations of  the  parameters  for  each  data  sample  to  determine  all 
"solutions"  and  to  then  determine  which  solution  gives  the  lowest  sum 
of  squares  of  Y residuals.  For  practical  reasons  this  is  very  diffi- 
cult to  do.  However,  one  must  remember  that  an  attempt  is  being  made 
to  find  a solution  to  a function  that  adequately  represents  the  data. 
Whether  or  not  there  are  solutions  in  other  unknown  regions  may  be 
rather  unimportant  if  the  solution  that  is  found  is  satisfactory  to 
the  analyst — that  is,  if  it  satisfies  the  analyst's  criterion  for  a 
good  fit.'*’ 


xFor  further  information  on  nonlinear  least-squares  solutions, 
see  N.  R.  Draper  and  H.  Smith,  Applied  Regression  Analysis , John  Wiley 
& Sons,  Inc.,  New  York,  London,  Sydney,  Chap.  10,  1966,  pp.  263-304. 


I 
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Appendix  B 

LEAST-SQUARES  ESTIMATION  FOR  ASYMPTOTIC-POWER  EQUATION 
Y = A + B • X1C 


REGRESSION  EQUATIONS 

To  obtain  least-squares  estimates  of  the  parameters  A,  B,  and  C 
of  the  asymptotic-power  equation,  the  following  procedure  is  used. 

First,  let  the  residual  corresponding  to  Y be  defined  by 

e.  = Y.  - Yci  = Y.  - (A  + B • XI J)  , (1) 

where  A,  B,  and  C are  least-squares  estimates  of  the  parameters. 

The  requirement  for  a least-squares  fit  for  N data  points  is  that 
the  sum  of  squares  of  the  Y residuals  (denoted  by  Q)  shall  be  a minimum; 
here. 


Q = 


N 


l 

i=l 


<Yi 


(2) 


If  Q is  to  be  a minimum,  the  partial  derivatives  of  Q with  respect  to 
the  parameters  A,  B,  and  C must  be  zero: 


3A  = 0A  " ^ ' l (Yi  “ A - B • XlJ) 

0B  = -2  • l (Y±  - A - B • XlJ) 

Qc  = ~2  ‘ I <Yi  " A ‘ B • XlJ) 

Simplifying  and  rearranging  terms  gives: 


|: 


• B • Xl^  • Ln  Xli  = 0 . 


1 
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L 


l = A * N + B • l XlJ  , (3) 

l Y±  * XlJ  = A • l XlJ  + B • l X1?C  , (4) 

T Y.  • X1C  • Ln  XI.  = A • y X1C  • Ln  XI.  + B • y X1?C  • Ln  XI . . (5) 

u i l l L l l L i l 

The  problem  then  becomes  one  of  solving  Eqs . (3),  (4),  and  (5)  for 
the  parameters  A,  B,  and  C,  given  a set  of  independent  observations  of 
Y and  XI.  Except  in  very  special  cases,  the  equations  cannot  be  solved 
by  ordinary  algebraic  methods  but  must  be  solved  by  iterative  tech- 
niques. First,  A can  be  eliminated  from  Eqs.  (3)  and  (4)  by  multiplying 
r C 

Eq.  (3)  by  \ Xl^  and  Eq.  (4)  by  N and  then  subtracting  the  two  equations. 
Having  done  this,  one  can  solve  for  B in  terms  of  C.  That  is 


B 


l Xic  - N • y Y • 

_ 1 L 1 


- N • l X12C 

L l 


(6) 


Therefore,  for  a given  set  of  observations  of  Y and  XI,  if  C is  known, 
B can  be  solved  from  Eq.  (6),  and  A can  then  be  solved  from  Eq.  (3). 


A = 


l Yi  - B * l X1C. 


(7) 


The  solution  of  A,  B,  and  C must  also  satisfy  Eq.  (5).  Let  G 
represent  the  difference  of  the  members  in  Eq.  (5)  as  follows: 


G = y Y • Xl^  • Ln  XI.  - A • y Xl^  • Ln  XI  - B • y X12C  • Ln  XI,  . 
‘•ii  i ^-i  i ri  i 


(8) 


G will  be  zero  only  when  A,  B,  and  C are  a solution. 


PROGRAM  SEQUENCE  OF  OPERATIONS 

The  sequence  of  operations  in  the  computer  program  is  as  follows. 

First,  the  various  summations  involved  in  Eqs.  (6),  (7),  and  (8)  are 

obtained  using  C = -8.001  (initially).  Then  B and  A are  determined 

2 

from  Eqs.  (6)  and  (7).  After  these  calculations  are  made,  the  value 
of  G is  obtained  from  Eq.  (8),  and  its  algebraic  sign  is  noted.  Unless 
A,  B,  and  C are  a solution,  G will  not  be  zero.  The  machine  then  steps 
the  value  of  C by  +0.1,  repeats  all  of  the  summations  and  calculations, 
and  checks  the  algebraic  sign  of  G again.  This  procedure  is  continued 
until  the  algebraic  sign  of  G is  reversed,  signifying  that  a solution 
lies  somewhere  between  the  previous  value  of  C and  the  value  of  C at 
this  cross-over  point. 

At  this  point,  the  program  begins  an  iterative  operation  in  which 
at  each  cross-over  point  the  incremental  step  is  halved  and  the  direc- 
tion of  advance  is  reversed.  This  iterative  procedure  is  done  as  many 
times  as  desired  to  give  any  degree  of  accuracy  required  for  C.  In 
the  program,  this  procedure  is  repeated  until  the  ratio  of  the  value 
of  each  of  the  parameters  A,  B,  and  C from  one  iteration  to  the  next 
differs  from  unity  by  an  amount  equal  to,  or  less  than,  10  ^ (or  as 
otherwise  specified). 

3 

The  search  for  roots  continues  to  C = -0.001.  After  this  point  , 
is  reached,  the  program  begins  another  search  starting  at  C = +0.001 
and  proceeding  by  increments  of  +0.1  out  to  +8.001.  If  no  solution 
at  all  is  found  within  these  limits,  a statement  to  this  effect  is 
printed,  and  the  program  continues  on  to  the  next  run.  Any  time  a 
solution  is  found  for  A,  B,  and  C,  the  sum  of  squares  of  Y residuals 
(Q)  is  determined  and  compared  with  the  corresponding  value  for  the 

^Acknowledgment  is  made  to  James  Johnston  (formerly  at  Rand)  for 
his  suggestions  in  the  initial  programming  aspects  of  this  problem. 

2 

If  A is  specified,  then  that  value  is  used  instead  of  calculating 
A from  Eq.  (7).  Also,  the  equations  for  B and  G are  changed.  However, 
the  procedure  of  solving  for  B and  C is  similar  to  the  case  in  which  A 
is  not  specified. 

3 

Because  a zero  value  for  C results  in  a degenerate  case,  the 
search  purposely  avoids  a region  very  close  to  zero  for  C. 


1 


1 
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previous  solution  (if  there  was  one) . The  solution  that  gives  the 
lowest  sum  of  squares  of  Y residuals  is  stored  temporarily  for  compar- 
ison with  any  future  solution  so  obtained.  In  this  way,  when  the 
search  is  completed  and  if  there  is  a solution,  that  solution  will 
generally  represent  the  lowest  sum  of  squares  of  Y residuals  in  the 
region  searched. 

Any  "solution"  found  in  the  specified  range  for  C represents  a 
solution  for  which  the  partial  derivatives  of  Q with  respect  to  the 
parameters  are  zero.  The  Q value  for  that  solution  is  also  compared 
with  the  Q values  for  the  end  points  of  C to  make  sure  that  Q is  not 
decreasing  to  some  other  minimum  outside  the  range  of  C.  As  of  now, 
we  have  not  been  able  to  determine  any  requirements  for  Q to  have  a 
unique  minimum  but  have  observed  that  for  various  sets  of  data,  Q 
seems  to  have  a unique  minimum  in  the  region  searched.  Even  if  it 
does  not,  the  minimum  of  the  relative  minima  will  usually  be  found. 

As  stated  before,  there  is  apparently  no  proof  that  other  minima  can- 
not exist  outside  the  range  searched,  which  cannot  be  determined  by 
the  above  method,  even  when  a solution  has  been  found  in  the  prescribed 
range.  However,  this  may  be  unimportant  if  the  "solution"  found  satis- 
fies the  analyst's  criterion  for  a good  fit. 

The  above  limits  on  C and  the  increments  of  0.1  were  chosen  on 
the  basis  of  what  is  believed  to  be  a reasonable  search  range  for  C, 
of  economic  computer  operating  time,  and  of  the  extent  to  which  the 
search  range  should  be  covered  in  order  to  lessen  the  chances  of  miss- 
ing a root.  Although  two  roots  could  conceivably  be  missed  in  the  in- 
crement of  0.1,  indicating  that  the  G function  goes  from,  say,  a posi- 
tive to  a negative  to  a positive  value  within  an  interval  of  C equal 
to  0.1,  this  seems  rather  unlikely.  Such  a function  would  have  to 
behave  very  erratically,  and  test  results  indicate  that  this  function 
does  not  behave  in  this  manner. 

Perhaps  it  should  be  noted  that  a degenerate,  or  trivial,  case 
results  if  C = 0 or  if  all  Y values  are  constant  or  if  all  XI  values 
are  constant.  Any  of  these  conditions  results  in: 


* 


I; 


Y = constant. 
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Appendix  C 

MODIFICATIONS  TO  CURVES 

INTRODUCTION 

The  purpose  of  this  appendix  is  to  report  on  several  modifications 

that  have  been  made  to  the  CURVES  Cost  Analysis  Curve-Fitting  Program, 

* 

reported  in  an  earlier  Rand  Report  R-1753-PR.  A new  listing  of  the 
modified  CURVES  computer  program  is  presented  in  Appendix  D. 


VARIABLE  TRANSFORMATIONS 

A major  modification  has  been  incorporated  in  CURVES  to  allow  for 


variable  transformations  of  the  following  kinds:  (a)  power,  (b)  log- 

arithmic, and  (c)  binary.  Each  of  these  is  discussed  below. 

Power 

A power  transformation  may  be  made  on  any  variable  as  follows: 

E 

V - (V  ) P , 

P P 

where  V = variable  p, 

P 

p = 0 through  7 (p  = 0:  Y variable,  p = 1:  XI  variable,  p = 2: 

X2  variable,  ...,  p = 7:  X7  variable), 

E = real-number  exponent  for  variable  p,  with  range  -10.0  < E < 10.0, 
P P 

-*■  = "transformed  to." 

To  fit,  for  example, 

1 / Y ■ A + B • XI  , 

-1.0  is  entered  for  Eg.  To  fit 

*H.  E.  Boren,  Jr.,  and  Capt.  G.  W.  Corwin,  CURVES:  A Cost  Analysis 

Curve-Fitting  Program , The  Rand  Corporation,  R-1753-PR,  December  1975. 


Y = A + B • XI  + C • , 


0.5  is  entered  for  To  fit 


Y = A + B / Nl  , 


-.333  is  entered  for  E^.  To  fit  an  equation  of  the  form 
Y=A+B  • X1+C/X2+D  / /X3  , 

data  for  Y and  the  three  X variables  are  entered  as  if  a linear  regres- 
sion is  to  be  run.  Then  a value  of  -1.0  is  entered  for  and  -0.5 
for  Ey  The  values  for  each  variable  are  then  raised  to  the  appropriate 
exponent  and  stored  in  the  same  cells  as  the  original  values,  thereby 
replacing  the  original  values.  For  CURVES,  any  time  a variable  is 
transformed,  the  original  values  of  the  variable  are  lost  for  that  run. 
However,  if  the  original  values  are  initially  stored  in  memory  or  on 
disk,  they  can  be  recovered  for  subsequent  runs.  Any  variable  so  trans- 
formed is  indicated  as  such  in  the  Table  of  Residuals  of  the  output  by 
the  word  MODIFIED  over  the  name  of  the  variable.  All  transformation 
factors  are  also  printed  at  the  end  of  the  Table  of  Residuals. 

Values  for  E are  entered  on  the  Control  card  (see  Table  4 — p.  20) 

P 

in  5-column  fields  beginning  in  Col.  41.  The  fields  are  listed  in  Table 
C.l  below. 


Table  C.l 

TR/  NSFORMATION  FIELDS  ON  CONTROL  CARD 


Columns 

Transformation 

Factor 

Variable  to  be 
Transformed 

41-45 

E0 

Y 

46-50 

E1 

XI 

51-55 

E2 

X2 

56-60 

E3 

X3 

61-65 

E4 

X4 

66-70 

E5 

X5 

71-75 

E6 

X6 

76-80 

E7 

X7 
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If  a 0 (zero)  or  1.0  is  entered,  or  the  field  is  left  blank,  no 
transformation  takes  place  for  the  variable  corresponding  to  that  field. 
The  implied  decimal-point  location  is  at  the  right  end  of  each  field; 
a punched  decimal  point  overrides  the  implied  location.  Similar  to  most 
of  the  other  information  on  the  Control  card,  any  transformation  factor 
that  has  been  entered  is  retained  from  run  to  run  unless  superseded  by 
a new  value.  Thus,  if  the  same  transformations  are  to  be  made  for  a 
series  of  regressions,  the  transformation  factors  need  only  be  entered 
for  the  first  run. 

Logarithmic 

Any  variable  may  now  be  transformed  to  its  logarithm  as  follows: 

V -*  Ln  (V  ) , 


where  V,  p = same  as  before, 

Ln  = natural  logarithm. 

Although  the  latest  published  version  of  CURVES  treats  logarithmic 
equations,  the  program  cannot  treat  logarithms  of  individual  variables. 
For  example,  it  is  now  possible  to  fit  an  equation  of  the  form 

Y=A+B*Xl+C*Ln  (X2)  + D • X3  , 

in  which  logarithms  are  taken  only  of  the  X2  variable.  To  designate  a 
logarithmic  transformation,  a value  of  88.0  is  entered  in  the  transforma- 
tion field  corresponding  to  the  variable  to  be  transformed  by  logarithms. 
Thus,  for  the  above  equation,  an  88.0  would  be  entered  for  in  Cols. 
51-55  of  the  Control  card  (Table  C.l).  In  that  case,  natural  logarithms 
would  be  taken  of  the  values  entered  for  the  X2  variable . 

Binary 

The  following  kinds  of  binary  operations  may  be  made  in  CURVES: 


1.  V + V 
P P 


’"'-■W'-w 


PI 


■ 
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t j 
' 

4 I 


J 


I 

i 

I 


2. 


V ->  V 
P P 


3.  V -*■  V + V , 
p p q 


4.  V -*  V - V , 

p p q 


where  V,  p 

q 


same  as  before, 

subscript  of  second  variable  involved  in  operation 
(0  through  7). 


To  effect  any  of  the  above  binary  operations,  a two-digit  number 
ranging  from  10.0  through  47.0  is  entered  in  the  field  corresponding 
to  the  variable  to  be  transformed.  The  left  digit  of  the  number  indi- 
cates which  of  the  above  binary  operations  is  to  take  place.  A 1 
signifies  multiplication,  a 2 division,  a 3 addition,  and  a 4 subtrac- 
tion. The  second  digit  is  the  q subscript  representing  the  second 
variable  involved  in  the  operation.  For  example,  suppose  that  the 
following  equation  is  to  be  fitted: 


(Y  / XI)  = A + B • XI  . 


A 21.0  would  be  entered  for  Eq  on  the  Control  card  in  Cols.  41-45  cor- 
responding to  the  Y field.  The  digit  2 in  21  signifies  that  divisions 
are  to  be  made  on  the  Y values,  and  the  digit  1 in  21  signifies  that 
the  divisors  are  to  be  the  XI  values.  In  this  case,  after  the  Y and  XI 
values  are  entered,  each  Y value  is  replaced  by  the  division  of  it  by 
the  corresponding  XI  value.  The  division  thus  becomes  the  new  Y value. 

Another  example  of  a binary  operation  is 

Y = A + B • (XI  / X3)  + C • (X2  + X4)  . 

Here,  Y and  the  "four"  independent  variables  XI  through  X4  must  be 
entered  through  the  usual  input  process  for  CURVES.  Thus,  on  the  Con- 
trol card  a linear  regression  with  four  independent  variables  would  be 
indicated.  However,  after  the  operations  are  performed  in  accordance 
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with  the  equation,  only  two  independent  variables  are  to  be  used  in 
the  regression,  namely,  (XI  / X3)  and  (X2  + X4)  . This  is  accomplished 
as  follows: 

In  Cols.  46-50  of  the  Control  card,  corresponding  to  the  XI  vari- 
able, the  two-digit  transformation  factor  23  is  entered  for  E^.  This 
signifies  to  the  program  that  each  value  of  the  original  XI  variable 
is  to  be  replaced  by  the  value  of  XI  / X3.  In  Cols.  51-55,  corre- 
sponding to  the  X2  variable,  the  two-digit  number  34  is  entered  for  E2. 
This  signifies  that  each  value  of  the  original  X2  variable  is  to  be 
replaced  by  adding  to  that  value  the  corresponding  value  of  X4.  To 
prevent  X3  and  X4  from  being  used  as  independent  variables  in  the  re- 
gression, a 99.0  is  entered  for  E^  in  the  next  field  on  the  Control 
card,  corresponding  to  the  X3  variable,  namely.  Cols.  56-60.  Any  time 
a 99.0  is  used  in  a transformation  field,  the  variables  corresponding 
to  that  field  and  to  any  remaining  transformation  fields  are  not  used 
in  the  regression.  Therefore,  if  a 99.0  is  used,  it  must  be  the  last 
transformation  factor  entered  on  the  Control  card.  Also,  whenever  a 
99.0  is  entered,  the  transformation  factors  for  that  field  and  for  the 
remaining  fields  are  set  to  zero.  This  prevents  unwanted  carryovers 
of  transformation  factors  from  previous  runs.  Obviously,  a 99.0  can 
never  be  entered  on  the  Control  card  in  the  fields  corresponding  to 
the  Y or  XI  variables  because  those  two  variables  are  always  required 
for  a regression.  If  a 99.0  is  to  be  entered  on  the  Control  card  in 

the  transformation  fields,  it  must  be  entered  for  an  E after  Col.  50 

9 

and  only  after  a transformation  factor  E has  been  entered  with  a value 

P 

in  the  range  of  10.0  fi  E S 47.0,  where  p < q. 

P 

In  the  above  example  in  which  Y is  regressed  against  (XI  / X3)  and 


(X2  + X4) , the  values  of  X3  and  X4  are  shown  in  the  Table  of  Residuals 
of  the  output  with  the  words  NOT  USED  over  the  (X3)  and  (X4)  headings, 
because  they  are  not  used  as  independent  variables  in  the  regression. 


SAMPLE  OUTPUTS 

Examples  of  outputs  resulting  from  regressions  involving  the  trans- 
formations discussed  here  are  given  in  Figs.  C.l  through  C.6.  The  first 
regression  involves  an  equation  of  the  form: 

Y = A + B/Xl  + C*Ln  (X2)  + D • + E • #K4  . 
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Fig.  C . 5 — Regression  of  Y = A + B • (XI  + X3)  + C • (X2  - X4)  + D 
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Data  for  the  examples  were  taken  from  precalculated  values  and 
therefore  represent  near-perfect  fits.  This  was  done  in  order  to  check 
the  results.  Because  XI,  X2,  X3,  and  X4  were  all  transformed,  the  word 
MODIFIED  appears  over  each  of  their  headings  in  the  Table  of  Residuals 
shown  in  Fig.  C.2.  The  transformation  factors  are  shown  at  the  bottom 
of  the  figure.  To  obtain  the  original  X values,  one  simply  reverses 
the  transformation  process.  That  is,  to  obtain  the  original  X values 
for  the  first  run,  the  XI  values  are  raised  to  the  -1  power,  the  X2 
values  are  exponentiated  (using  base  e),  the  X3  values  are  squared,  and 
the  X4  values  are  cubed.  However,  as  was  mentioned  before,  if  the 
original  values  are  stored  in  memory  or  on  disk,  they  can  be  reread 
(recovered)  for  another  run. 

The  third  page  of  output  (Fig.  C.3)  shows  a regression  involving 
binary  operations  of  the  form: 


Y = A + B • (XI  / X3)  + C • (X2  • X4)  . 

In  the  Table  of  Residuals,  the  word  MODIFIED  appears  over  the  XI  and  X2 
headings,  because  those  values  were  transformed.  The  values  of  X3  and 
X4  are  listed  but  with  the  heading  NOT  USED  over  (X3)  and  (X4)  to  indi- 
cate that  they  were  not  used  as  independent  variables  in  the  regression. 
Because  each  original  value  of  XI  was  divided  by  each  value  of  X3,  the 
original  values  of  XI  can  be  calculated  as  follows: 

XI  (new)  = XI  (old)  / X3  , 
or  XI  (old)  = XI  (new)  • X3  . 


Therefore,  if  each  XI  value  shown  in  the  Table  of  Residuals  is  multi- 
plied by  the  corresponding  value  of  the  X3  variable,  the  original  value 
of  XI  is  obtained  for  each  data  point.  A similar  process  can  be  used 
to  obtain  the  original  values  of  X2.  Figure  C.4  shows  another  output 
involving  a regression  of  an  equation  of  the  form: 


Y = A + B • (XI  + X3)  + C • (X2  - X4)  . 
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Lastly,  Figs.  C.5  and  C.6  show  an  extension  of  the  results  of  Fig. 
C.4  in  which  X3  and  X4  are  allowed  to  remain  in  the  regression.  The 
equation  is  then: 

Y = A + B • (XI  + X3)  + C • (X2  - X4)  + D • X3  + E • X4  . 

Because  X3  and  X4  are  allowed  to  remain  in  the  regression,  the  heading 
NOT  USED  does  not  appear  in  the  X3  and  X4  headings  in  the  Table  of 
Residuals.  A listing  of  the  input  data  for  the  four  runs  shown  in 
Figs.  C.l  through  C.6  is  shown  in  Fig.  C.7.  For  convenience,  card 
column  numbers  are  shown  at  the  top  of  the  figure. 

A summary  of  the  transformation  factors  is  given  in  Table  C.2.  An 
updated  listing  of  the  CURVES  program  including  a new  subroutine  TRANS 
is  included  in  Appendix  D. 

VARIANCE-COVARIANCE  MATRIX 

Except  for  a one-parameter  case,  the  CURVES  program  now  prints  (in 
scientific  notation)  the  variance-covariance  matrix  of  the  estimated 
coefficients  on  the  first  page  of  the  output.  For  a one-parameter  case, 


Y = B • XI  , 

the  variance  of  B is  simply  the  square  of  the  standard  error  of  B,  al- 
ready printed  at  the  top  of  the  page  to  the  right  of  the  value  of  B. 

MISCELLANEOUS 

For  the  power  and  exponential  cases,  the  default  value  of  the  it- 
eration limit  has  been  changed  from  20  to  100.  It  was  found  that  a 
limit  of  20  is  not  always  sufficient  for  such  regressions. 
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10  20  30  40  50  60 

12345678  90 12 34567  890 1 23 45 o7 890 1 2 34567 8 90 1 2 3 45678 90 1 2 34567890 1 2 3 4 5 < COLO  HNS 


TEST  BUN  5 

--  Y = 

A 

♦ B/X  1 

♦ C * LN  ( X 2) 

♦ D * 

SORT ( X 3 ) ♦ 

E * CUBE  ROOT ( X4 ) 

1Y1234 

1 

-1 

8ft 

.5. 3333 

FORMAT  (5F10.0) 
BEAD  ONCE 

11.2691599 

2.567 

8.725 

1. 276 

23.  1 768 

1 1.  1324268 

12.  738 

1.2876 

.8926 

10.287 

19.9600498 

23.987 

35.276 

17.827 

9.287 

13.2325164 

15. 

10. 

1 . 

3. 

13.0287067 

25.762 

10. 

. 538 

. 92c 

14. 2181985 

30. 127 

.3872 

13.287 

18.298 

13.9622199 

1. 

10. 

20. 

30. 

19. 3950887 

26.387 

11.251 

19.555 

7.295 

12.3779224 

5.332 

21.637 

. 3724 

10.278 

TEST  BUN  6 

— Y = 

A 

♦ B » 

( X 1 / X 3 ) ♦ C 

* 

( X 2 *X  4) 

23 

14 

99 

READ 

7.21034603 

2.367 

.8726 

8. 271 

1 . 345 

14. 6467017 

10.991 

5.2881 

.87165 

1 2. 382 

73.0106287 

. 8276 

10. 276 

5.287 

8.297 

190. 143445 

15.765 

20.  764 

9.111 

11.701 

7.79735533 

. 987 

. 3654 

10. 

3. 

126. 33447  . 

9436 

10. 567 

8.351 

14.54  t 

TEST  RUN  7 

--  Y = 

A 

♦ B * 

( X 1 ♦ X 3 ) + C 

* 

( X 2 - X 4 ) 

33 

44 

99 

READ  MEMORY 

4. 21554396 

.6743 

20. 

2. 5b7 

9. 554 

2.5090353  1 

.785 

25.784 

4. 678 

3. 224 

1.8302536  3 

.653 

15.  922 

1 . 287 

1 . 

15.9004728 

.332 

20. 

1. 

3. 

16.294881  . 

785 

18.431 

.555 

. 888 

7.8104  1 

17. 

2. 

3. 

4.76103915 

2.111 

14.045 

.882 

4. 

TEST  RUN  8 

--  Y = 

A 

♦ B * 

(XltXI)  ♦ C 

♦ 

( X 2 -X  4) 

♦ D 

* X 3 

♦ E * X4 

33 

44 

BEAD 

DONE 


Fig.  C.7 — Input  card  arrangement  to  generate  outputs  shown  in 

Figs.  C.l  through  C.6 
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Table  C.2 

SUMMARY  OF  TRANSFORMATION  FACTORS 


Type  of 

Factor  (E  ) Transformation 

p 

E 

-10.0  < E < 10.0  (V  ) p 

P P 

E = 88.0  Ln  (V  ) 

P P 

10.0  £ E £17.0  V • V 

p p q 

20.0  £ E £ 27.0  V /V 

p p q 

30.0  £ E £ 37.0  V +V 

p p q 

40.0  £ E £ 47.0  V - V 

p p q 

E = 99.0  Do  not  use  V or  any  remain- 

P P 

ing  variables  in  regression. 

Ep  and  remaining  transforma- 
tion factors  (E  ,,,  ...  E.,) 
p+i  / 

set  to  zero. 

E = other  values  Error 

P 

NOTES:  V = variable  p. 

P 

Vq  = variable  q;  q is  designated  by  second  digit 
of  E^;  first  digit  of  E^  (1,  2,  3,  or  4) 
designates  type  of  combination — multiplica- 
tion, division,  addition,  or  subtraction, 
respectively,  when  10.0  £ E^  £ 47.0. 

E = transformation  factor  for  variable  V . 

P P 

p,  q = 0 through  7 (0  = Y variable,  1 = XI  variable, 

2 = X2  variable,  ...,  7 = X7  variable). 

Ln  = natural  logarithm. 


u u u u 
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Appendix  D 


UPDATED  LISTING  OF  CURVES  PROGRAM 


CURVES:  A COST  ANALYSIS  CUBVE-FITTI NG  COHPUTEF  PBOGBAB,  BAND  HAIN0010 

REPORT  B-1753-1-PB,  EY  H. E.  BOREN,  JR.  AND  G.B.  CORWIN,  HAIN0020 

SEPTEHBER  1976  HAIN0030 

HAIN0040 

IBPLICIT  REAL* 8 (A*H, C-Z)  HAINC050 

TBE  FOLLOWING  IS  THE  CONPLBTE  COEBON  HAIN0060 

COHHON  /Cl/  W,  N RAX,  KOLUBN,  NIV,  NIVP1,  NVAR , NP,  LABEL,  ID  HAIN0070 

COHBON  /C2/  IV  AL  (10)  , IND,  IFV,  IN,  18,  JX  (9)  , NPAGE  HAIN0080 

COHHON  /C3/  ABORT,  LINEAR,  IA,  IGUESS,  NXSIT,  NYCOL,  NYC,  NY  DEV  HAIN0090 

COHHON  /C4/  S (8) , SYX(8,8),  COV(0,8),  FHT(30),  HEAD  (9) , LNH (8)  HAIN0100 

COHHON  /C5/  AN,  P(8),  SE(8),  TB(8),  SIGLEV(8),  BETA  (8) , VHSQ(8)  HAIN0110 

COHHON  /C6/  VHEAN  (8)  , T (8)  , H(8,8),  PI  (8)  , SDEV  (8)  , RH{8,8)  HAIN0120 

COHHON  /C7/  YCEPT , DELTA,  RDEVH,  EA,  DW,  XV,  YV,  SEYSQ  HAIN0130 

COHHON  /C8/  DPT,  DPI,  DP2 , CD,  SEY,  CV,  YDEVSQ,  FVALDE,  SST  HAIN'im 

COHHON  /C 9/  7(2211)  HAINCISO 

COHHON  /CO/  DATA  (10,101)  HAIN0160 

COHHON  /Cl/  EX V ( 8) , NEXV,  NVT  HAIN01 70 

LOGIC AL*4  ABORT,  LINEAR,  IA,  IGUESS  HAINC180 

EQUIVALENCE  (IEQ,  IVAL(1)),  (NCARDS,  IVAL(9)),  (LIH,  IVAL(IO))  HAINC190 

DO  10  I = 1,  8 HAINP20O 

EXV(I)  * 0.D0  HAIN0210 

B H (I,  I)  = 1 . DO  HAIN02  20 

10  IVAL(I)  = 0 HAIN0230 

NCARDS  - 1 HA INC 240 

LIH  = 100  HAIN0250 

NPAGE  * 1 HAIN0260 

DELTA  = 1.D-7  HAIN0270 

WRITE  (6,  20)  HAIN0280 

20  PORHAT  ( 1 H 1 , 42X,  'CORVES  REGRESSION  ANALYSIS  COHPUTER  PROGRAH'/  HAIN0290 
1 59X,  '(JULY  1976)'  ) HAIN0300 

SET  ABOBT-POR-ERROR  INDICATOR  TO  FALSE  HAIN0310 

30  ABORT  = .FALSE.  HAIN0320 

CALL  READ  HAIN0330 

CALL  INPUT  (V)  HAIR0340 

CHECK  ERROR  INDICATOR  HAIN0350 

IF  (ABORT)  GO  TO  30  HAIN0360 

IP  (IEQ  .EQ.  0)  GO  TO  80  HAIN0370 

CALL  PRINT  HAIN0380 

CALL  SUHS(V)  HAIN0390 

GO  TO  (40,  40,  50,  60,  50,  40,  40,  40),  IEQ  HAIN0400 

40  CALL  LINE  ( V)  HAIN04IO 

GO  TO  70  HAIN0420 

50  CALL  EXPO(V)  HAIN04 30 

GO  TO  70  HAIN0440 

60  CALL  ASYH(V)  HAIN0450 

70  IP  (ABORT)  GO  TO  30  HAIN0460 

CALL  STAT(V)  HAIN0470 

CALL  TVAL(V)  HAIN0480 

CALL  OUT1  HAIN0490 

80  CALL  OOT2 ( V)  HAIN0500 

GO  TO  30  HAIN05 10 

END  HAIN0520 
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LABEL,  ID 


1YC,  NY  DEV 


/C7  / 
/Cl/ 


SUBROUTINE  READ 
IMPLICIT  REAL*8 (A-H.O-2) 

COHHON  /Cl/  HI,  NHAX,  KOLUHN,  NIV,  NIVP1,  NVAfi,  NP, 

/C2  / IVAL(IO),  (12(2)  , IN,  18,  JX  (»)  , NPAGE 
/C3  / ABORT,  LINEAR,  IA,  IGUESS,  NXSST,  NYCOI., 

/C4/  7,1(136),  PHT ( 30 ) , H EA D ( 8 ) 

/C6 / 22  (80)  , PI  (8) 

YCEPT,  DELTA 
E XV (8) 

ABORT,  LINEAR,  IA,  IGUESS 

EX  V 1 (8)  , EXV2  (8)  , SVHEAD  (9)  , KX(9),  K V (10) 

ALPH  A ( 4)  , IVALI(IO),  IPHA(IO),  A NA  HES (6 ) , 

(LIST1  , ALIST  (2)  ) , (IEO,  IVAL  (1)  ) 

KV/'  Y ' , ' 1 • , *2' , • 3' , '4 ', • 5' , ' 6’ , ■ O' , ' • / 
LONCE/'ONCE' /,  HEHO/' HEHO' /,  LDISK/'DISK'/ 
PARLB/'  X 1 **2 ' /,  EXPGN/' EXPONENT'/ 

IbLANK  /•  '/,  BLANK/'  '/ 

SVHEAD/' LABEL', ' Y*,'  XI', • X2','  X3», 

• X 4 * , • XS'  , • X6 ' , ' 11' / 

ANAHE S/'READ',' FORHAT', 'LA BEL',' LABELS' , 'GUESS' , 
'READS'/,  JXSAVE/'  0 ',' K ',' D ',' E ',' R ',' ?• 


COHNJN 
COHHON 
COHNON 
COHHON 
COHHON 
COHHON 
LOGICAL*  4 
DIHEN SION 
DIHENSION 
EOUIV  ALENCE 
DATA 
DATA 
DATA 
DATA 
DATA 


, JXSAVE  (9) 
ALIST(IO) 


DATA 


SUBROUTINE  FOR  READING  TITLE,  CONTROL,  AND  SPECIFICATION  CARDS 


IVAL: 


(1)  (2) 

IEQ  N P 1 


(3) 
NP  2 


(8) 

NP  3 


(5) 

L22YES 


(6) 

KEEP2 


(7) 

LNO'JT 


(8) 

IORD 


(9)  (10) 

NCARDS  LIH 


C 

C 

c 


CALL  TITLER 

READ  (5,  10)  KX,  L PH  A , ALPHA,  I V A 1. 1 , YCEPT,  DELTA1,  EX  V 1 , EX  V2 
10  FORHAT  (1X,9A1,T1,A1,9X,8A1,A2,4A5,T1,I1,9X,8I1,I2.F10.0,F10.10, 

1 8P5.0,  T41,  8AS) 

CHECK  EACH  FIELD  OF  CONTROL  CARD  FROH  COLUHN  1 THROUGH  COLUSN  40. 
IF  NOT  BLANK,  SET  APPROPRIATE  DESIGNATOR  TO  VALUE  SO  ENTERED. 
(VARIABLE  TRANS FOR  HAT  ION  FACTORS  ARE  ENTFRED  IN  COLUHNS  41-80.) 


20 


I BL  ANK)  IVAL  (J)  = IVAL  1 (J  ) 


.OR.  ALPHA ( 1 ) . 

.or.  alpha  (3)  . 
DELTA  .LT.  1.D- 


30 


40 

50 


60 

70 

90 


DO  20  J = 1,  10 
IF  (LPHA(J)  .NE. 

CONTINUE 
IA  = .FALSE. 

IP  ( A LPH  A ( 2)  .NE.  BLANK 
IF  (A  LPHA  (4)  .NE.  BLANK 
IF  (DELTA  .GT.  0.1  .OR. 

I FHT  = -10 
LABEL  = 0 
IN  = 0 
18  = 5 

IGUESS  = . FALSE. 

DO  30  I = 1,  9 
HEAD(I)  = SVHEAD  (I) 

IF  (I  EC  .EQ.  2)  HEAD  (4) 

READ  (5,  50)  ANAHE,  ALIST 
FORHAT  ( A6 , A2,  -\8) 

DO  60  J =1 , 6 
IF  (ANAHE  .EQ.  A NAHES (J) ) 

CONTINUE 
WRITE  (6,  80)  ANAHE,  ALIST 
FORHAT  ( '0  ILLEGAL  NAHE  OH  SPECIFICATION 


NE. 

NE. 

12) 


BLANK) 
BLANK) 
DELTA  = 


IA  = . 
DELTA 
1.  D-l 


TRUE  . 

= DELTA  1 


= PARLB 


GO  TO  (160,  90,  110,  110,  130,  150),  J 


CARD.  THIS  JOB  HAS  ', 


1 'BEEN  TERHI NATED. '/  1H  , A6,  A2,  9A8) 
STOP 

90  I FHT  = I FHT  ' 

IF  (IFHT  . tj  . 


10 

30) 


GO  TO  70 


PEADOO 10 
RE AD0020 
RE  ADOO  30 
READO04C 
RE  ADOO  50 
RE ADOO 60 
RE AD0070 
P E ADOO  80 
H E ADOO  90 
READ01 00 
READ91 10 
READ01 20 
READ01 30 
READ01 40 
READ01 50 
RE  ADO  1 60 
R E AD0 1 70 
PEAD01 80 
READ0190 
READ0200 
READ02 10 
READ0220 
READ0230 
READ0240 
READ0250 
READ0260 
READ0270 
READ0280 
READ0290 
R E ADO  3 00 
READ 0310 
READ0320 
READ0330 
R E ADO  3 40 
READ0350 
RE  ADO  360 
RE  ADO  3 70 
READ03  80 
RE  ADO  3 90 
READ0400 
READ04 10 
READ0420 
READ04  30 
READ0440 
READ04  50 
READ0460 
READ0470 
READ0480 
READ0490 
READ0500 
RSADC5 10 
RE  ADOS  20 
R E ADOS  30 
READ0540 
READ0550 
READ0560 
READ0570 
READ0580 
RE  ADOS  90 
READ0600 


BE 


I 
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DO  100  I = 1,  10 

REAE06  1 0 

100 

FMT(I*IFMT)  = AL1S1(I) 

READ0620 

GO  TO  40 

READ0630 

110 

LABEL  = J - 2 

REAC0640 

DO  120  1-1,9 

READ0650 

120 

riEAD(I)  = ALIST  (I*  1) 

READ0660 

GO  TO  40 

READ0670 

130 

IGUESS  = .TRUE. 

READ0680 

C 

CONVERT  ALPHANUMERIC  TO  NUMERIC 

READ0690 

CALL  MEMRE  (ALIST  (2)  , 64) 

PEAD0700 

READ  (1,  14C)  PI 

READ07  10 

140 

FORMAT  (6F8.0) 

READG720 

GO  TO  40 

READ0730 

150 

18  = 8 

READ0740 

160 

IP  ( L 1ST  1 .EQ.  LONCE)  IN  = 1 

READ0750 

IF  (LIST  1 .EQ.  MEMO  ) IN  = 2 

READ0760 

1 

IP  ( L 1ST  1 .EQ.  LDISK)  IN  = 3 

READ0770 

IF  (I  EQ  .LQ.  4)  HEAD  (4)  = EXPON 

READ0780 

! 

IP  (NPAGE  .GT.  2 .OR.  IFMT  .GE.  0)  GC  TO  180 

READ0790 

WRITE  (6,  170) 

READ0800 

170 

FORMAT  (//'9FI  ST  RUN  MUSI  HAVE  FORMAT  CARC.  THIS  JOB  HAS 

REAE0R1 0 

1 'BEEN  TERMINATED. ') 

READ0820 

STOP 

READ0830 

c 

USE  PREVIOUS  ORDER  IF  V A RI A3 1 E -CRDE R INDEX  AIL  BLANK 

REAE0840 

180 

DO  190  I = 1,  9 

READ0850 

IF  ( K X ( I ) .NE.  IBLANK)  GO  TO  20C 

READ0860 

190 

CONTINUE 

READ0870 

i 

c 

ERROR  IF  FIRST  RUN  DOES  NCT  HAVE  V A B I A E LE-CR D EP  INDEXES 

READ0880 

; 

IF  (NFAGE  .LE.  2)  GO  TO  3 30 

REA  DOR  90 

1 

GO  TO  290 

R EAD0900 

c 

DETERMINE  VARIABLE  ORDER 

READ09 10 

t 

j 

200 

JMAX  = 2 

READ0920 

NBLNKS  = 0 

READ0930 

1 

DO  250  K = 1 , 9 

READ0940 

DO  210  J = 1,  10 

READ0950 

IF  ( K X ( K ) .EQ.  K V (0)  ) GO  TO  240 

READ0960 

210 

CONTINUE 

READ0970 

GO  TO  340 

READ0980 

J 

240 

JXSAVE(K)  = J 

READ0990 

! 

IF  (J  .EQ.  10)  NBLNKS  = 1 

READ1000 

c 

EUR  OR  FOR  IMEEDDED  ELANKS  IN  VARIABLE-ORDER  INDEX 

READ  10  10 

IF  (NBLNKS  .GT.  0 -AND.  J .LT.  10)  GC  TC  340 

READ1020 

IF  (J  .LE.  8 .AND.  JMAX  .IT.  J)  JMAX  = J 

READ1030 

250 

CONTINUE 

READ1040 

c 

ERBOR  FOR  A REPEATED  VARIABLE  INDEX 

READ1050 

DO  260  1=1,9 

READ1060 

IP1  = I ♦ 1 

READ1070 

DO  260  J = IP1,  9 

READ  1080 

IP  (JXSAVE(I)  .NE.  JXSAVE(J))  GO  TO  260 

READ1090 

IF  (JXSAVE(I)  .NE.  10)  GO  TO  340 

READ  1 1 00 

260 

CONTINUE 

READ1110 

c 

ERROR  IP  ALL  APPROPRIATE  VARIABLES  NCT  SPECIFIED 

REAC1120 

DO  230  J = 1 , JMAX 

READ1130 

DO  270  I = 1 , 9 

READ  1 1 40 

• 

IF  (JXSAVE(I)  .EQ.  J)  GO  TO  280 

READ1150 

270 

CONTINUE 

READ1 160 

GO  TO  340 

READ1170 

280 

CONTINUE 

READ  1 1 80 

c 

NUMBER  OF  INDEPENDENT  VARIABLES 

BEAD  1 190 

NIV  = JMAX  - 1 

READ1200 

IF  (NIV  .LE.  0)  GO  TC  340 

READ1210 

Mtl 
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NIVP1  = JM  AX 

READ1220 

NIVP2  = JMAX  ♦ 1 

REAC12  30 

c 

TOTAL  NUMBER  OF  PARAMETERS  (NP) 

READ1240 

290 

NP  = NIVF1 

READ1250 

IF  (NIV  .GT.  1 .AND.  (IFQ  .EC.  2 .OP.  IEC  .EC.  4))  GO  TO  340 

REAB1260 

IF  (IEC  -EC.  2 .CR.  lEv  .EC.  4)  NP  = 3 

READ127C 

D)  292  I = 1 , NIVP1 

READ1280 

IF  ( EXV2  (I ) .NE.  BLANK)  EXV(I)  = E X V 1 ( I ) 

READ1290 

292 

CONTINUE 

READ1 300 

ID  = 0 

READ1310 

DO  300  1=1,  NIVP2 

READ  13  20 

JX(I)  = JXSAVE(l) 

READ1330 

IF  (ax  (I)  .NE.  9)  GO  TO  300 

READ  1 340 

JX(I)  = NIVF2 

READ1 350 

ID  = 1 

READ1360 

300 

CONTINUE 

READ  1 370 

IF  (LABEL  .EC.  2)  HRITE  (6,  310)  (HEAD(J+1),  J = 1,  NIVP1) 

READ1380 

310 

FORMAT  (1H0,  21X,  A 8 , 1 WITH  7(A8,2X)) 

READ1390 

| 

IF  (ID  . EC.  0)  HEAD  ( 1)  = ELANK 

READ1400 

NVAR  = N1VP1  ♦ ID 

READ1410 

NYC  = NVAR  ♦ 1 

READ1420 

IF  (IEC  .EC.  2)  NYC  = NYC  ♦ 1 

READ1430 

NYDEV  = NYC  + 1 

READ1440 

KOLUMN  = NYDEV 

READ1450 

IF  (IEC  .GE.  3)  KOLUMN  = KOLUMN  ♦ NIVP1 

READ  1460 

IF  (IEC  .EC.  5 .OR.  IEC  -EC.  7)  KOLUMN  = KOLUMN  - NIV 

READ1470 

NMAX  = 2211 

READ1480 

C 

REDUCE  V-SIZE  BY  676  IF  PLOTTING. 

READ1490 

IF  ( I V A L (2)  ♦ IV  A L (3)  ♦ IVAL<4)  .GT.  0)  NMAX  = 1535 

READ1500 

NMAX  = NMAX  / KOLUMN 

READ1510 

NXSET  = 0 

READ1520 

j 

IF  (IEC  .Eg.  3 .OR.  IEC  . EC.  6 .OR.  IEC  .EC.  0)  NXSET  = NYDEV 

READ  1530 

! 

I 

NYCOL  = 1 

READ1540 

IF  (IEc  .EC.  J .OR.  (IE0.GE.5  .AND.  IEC.LE.7))  NYCOL  = NYDEV  ♦ 1 

READ155C 

IF  ( . NOT . I A .OR.  (IEC  -NE.  3 .AND.  IEC  .NE.  5))  RETURN 

READ1560 

WRITE  (b,  320) 

READ1570 

32  0 

FORMAT  ('OINTERCEPT  MAY  NOT  EE  SPECIFIID  FOR  THIS  FUNCTION.', 

READ1580 

1 

/ 1 H , 'THIS  JOB  HAS  BEEN  TERMINATED.') 

READ1590 

STOP 

READ1600 

1 330 

WRITE  (6,  350)  JXSAVE 

READ1610 

stop  READ1620 

340  WRITE  (6 , 3^0)  KX  READ1630 

350  FORtlAT  (//'  OTHER  E IS  AN  ERROR  IN  THE  V A F I A BLE -ORDER  INDEX  (',  READ1640 

1 'CONTROL  CARD)  ***', 9A 1 ,'*♦*'/ • THIS  JOE  HAS  BEEN  TER  HI N ATED . • ) REAC1650 
STOP  READ1660 

END  READI670 


it**-. 
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: ' 

SJBROUTI Nt  TITLER 

TITL0010 

COMMON  /C2/  Ml (23),  NPAGE 

7TTL0020 

REAL*4  TITLE  (20),  DATE(2) 

TITLOO  30 

DATA  IB/*  •/.  IH 1/* 1 '/,  DONE/1 

LONE*/ 

TITL0040 

READ  (5,  10)  TITLE 

TITL0050 

10  FORMAT  ( 20 A 4 ) 

TITL0060 

IF  (TITLE(I)  . EG  • LONE)  GC  TO 

30 

TITL0070 

IF  (NPAGE  .EC-  1)  CALL  DATER  (DATE) 

TITL0080 

< 

ENTRY  TITLE2 

TITL0090 

WRITE  (6,  20)  In,  TITLE,  DATE, 

NPAGE 

TITLO  100 

20  FORMAT  (A1/1H  , 2 0 A4 , 10X , ■ EAT E : 

',Z5,4X,' 

TIME:  * , 24 

, 4 X,  'PAGE: 

,13/) TITL01  10 

NPAGE  = NPAGE  + 1 

TITL0120 

IH  = IH 1 

TITLO 1 30 

; 

RETURN 

TITLG140 

C PRINT  TERMINATION  STATEMENT  IF 

ALL  DATA 

HAVE  EEEN 

PROCESSED. 

TITL0 1 BO 

f 

30  WRITE  (b,  40) 

TITL0160 

■ 

40  FORMAT  ( 1 H 1 / 10(1  HO/)  , 

TITLO  17C 

1 44  X , *******  **** 

* * 

*********  y 

TITL0180 

2 44X,  •*  * * * 

**  * 

* 

' / 

TITLO 1 90 

3 44X , '*  * * * 

* * ♦ 

* 

' / 

TITL0200 

m 

4 44X , •*  * * * 

♦ * * 

****** 

' / 

TITL02 1 0 

» 

5 44X,  '*  * * * 

* * * 

♦ 

' / 

TITL0220 

6 441,  '*  * * * 

* * * 

* 

' / 

TITL0230 

7 44X , •*  * * * 

* ♦♦ 

♦ 

’ / 

TITL0240 

844X,  *******  **•* 

♦ * 

*********  » 

TITL025C 

STOP 

TITL0260 

END 

TITL0270 

: 


1 


non 
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SUBROUTINE  INPUT(V) 
IMPLICIT  R E AL  *8 ( A - H, C-Z) 


DIMENSION 
COMMON  /Cl/ 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
LOGICAL*4 
DIMENSION 
DATA 


/C2/ 

/C3/ 

/C4/ 

/C5/ 

/C8/ 

/CX/ 

/CO/ 


V (NMAX.KCLUMN) 

N,  NMAX,  KOLUMN, 
M3(5) , KEEFZ,  M4, 
ABORT,  M 6 , IA 
Z 1 {1  36)  , FMT(30) 
AN 

DFT,  DEI,  DF2 
E XV (8)  , NEXV,  NVT 
DATA  (10,101) 
AEORT,  IL 
V DAT  A (10) 

E LAN  K/ • •/,  IM/5/ 


Ml,  N I V 1 1 , NVAR,  NP,  M2,  ID 
ICRD,  NCARES,  Ml 5 (3)  , IN,  18, 


, ABLANK/’BLANK’/,  1101/101/ 


JX(9) 


SUBROUTINE  FOR  READING  IN  AND  SAVING  DATA 


IF  (IN  .EQ.  O')  GO  TO  160 
IF  (IN  .EQ.  1)  GC  TO  150 
KOUNTB  = 0 

IF  (IN  .GE.  3)  GC  TO  80 
IBYTES  = 80  ♦ NCARDS 

C SET  INPUT  MEDIUM  NUHEEH  FOR  MEMORY. 

IM  = 1 

C SAVE  DATA,  LATER  IT  WILL  BE  READ  FROM  MEMORY  USING  ’MEMRE’. 

DC  40  I = 1,  1101 

READ  (18,  10)  (DATA(K,I),  K = 1,  10) 

10  FORMAT  (10A8) 

C CHECK  FOR  BLANK  CARD(S)  (MUST  BE  COMF1 ET ELY  ELANK)  OR  ’BLANK’. 

IF  (DATA  (1,1)  .EQ.  AELANK)  GC  TC  60 
DO  20  K = 1,  10 

IF  (DATA  (K, I)  .NE.  BLANK)  GO  TC  30 
20  CONTINUE 

KOUNTB  = KOUNTB  » 1 
IF  (KOUNTB  . EQ.  NCARDS)  GO  TO  160 
GO  TO  40 
30  KOUNTB  = 0 
40  CONTINUE 

50  NMAX  = 1101  / NCARES  - 1 
GO  TO  26C 

60  JCARDS  = I * NCARDS  - 1 

IF  (JCARDS  .GT.  1101)  GO  TO  50 
DO  70  J = I,  JCARDS 
DO  70  K = 1,  10 
70  DATA  (K, J)  = BLANK 
GO  TO  160 

C SET  INPUT  MEDIUM  FOR  DISK. 

80  IM  = 4 
REWIND  4 

NR  EC  = NCARDS  * NMAX 
DO  110  1=1,  NR  EC 

C READ  INPUT  DATA  AS  ALPHANUMERIC  DATA 

READ  (18,  10)  VDATA 

IF  (VDATA(I)  .EQ.  ABLANK)  GO  TO  120 
C WRITE  INPUT  DATA  ONTC  UTILITY  DISK. 

WRITE  (4,  10)  VDATA 
DO  90  K = 1,  10 

IF  (VDATA  (K)  .NE.  BLANK)  GO  TO  100 
90  CONTINUE 

KOUNTB  = KOUNTB  ♦ 1 


I NPTO0 10 
IN  PTOO  20 
INPT0030 
INPT0040 
INPT0050 
INPT0060 
I NPT0070 
INPT0080 
INPT0090 
INPT0100 
INPT01 10 
INPT0120 
INPTO 1 30 
INPT0140 
INPT01 50 
INPTO 160 
INPT0170 
INPT01 80 
INPT01 90 
IN  PT0200 
INPT02 10 
INPT0220 
INFT0230 
INPT0240 
I N PT  02  50 
INPT0260 
INPT0270 
INPT0280 
IN  PTC290 
INPT0300 
INPT0310 
INPT032C 
INPT0330 
INPTO 340 
INPT0350 
INPT0360 
INFT0370 
INPT0380 
INPT0390 
INFT0400 
INPT04 10 
INPT0420 
INPT0430 
INPT04U0 
INPT0450 
INPT0460 
INPT0470 
INPT048C 
INPT0490 
INPT0500 
INPT0510 
INPT0520 
INPT0530 
INPT0540 
INPT0S50 
IN  PT0560 
INPT0570 
IN  PT0580 
INPT0590 
INPTOfeOO 


n n n n 





- - 
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100 

110 

120 

130 

140 


150 

160 


170 


180 


190 

200 


210 


212 

214 

216 


IF  (KOUNTB  .EC.  NCARDS)  GC  TC  160 

GO  TO  110 

KOUNTB  = 0 

CONTINUE 

GO  TO  260 

DO  1 30  K = 1 , 10 

VDATA(K)  = BLANK 

DO  1 40  J = 1 , NCABCS 

WRITE  (4,  10)  VDATA 

GO  TO  160 

SET  INPUI  MEDIUM  FOE  CARDS. 

IM  = 5 

KRECRD  = 1 - NCARDS 
IF  (IM  .GE.  5)  IM  = IP 
IF  {IM  .EC.  4)  REWIND  4 

IF  SUBROUTINE  'MEMRE'  IS  NOT  ASSEMLLED,  00  NCT  USE  'READ  MEMORY' 
OR  'GUESS'  SPECIFICATION  OPTIONS.  FURTHER,  THE  DIMENSIONS  ON 
COMMON  /CO/  MAY  EE  REDUCED  TO  'DATA (1,1)'  IN  MAIN  AND  INPUT. 
CREATE  DUMMY  FROGRAMS  FOR  'MEMRE'  AND  'DATER'. 

NEG  = 0 

DO  200  I = 1 , NM  A X 
KRECRD  = KRECRD  ♦ NCARDS 

IF  (IM  .EC.  1)  CALI  MEMRE  (DATA  ( 1,  KRECRD)  ,IBYTES) 

READ  INPUT  DATA  FROM  CARDS,  FROM  MEMORY,  OP  FROM  UTILITY  DISK 
READ  (IM,  FMT)  ( V ( I , JX(J)),  J = 1,  N V A R ) 

NZ  = 0 

DO  180  K = 1,  NIVP1 

IF  ( V (I  , K)  .EC.  0 . DO ) NZ  = N Z + 1 

IP  (V  (I , K)  .LT.  0 . DO  .AND.  KEEPZ  .NE.  2)  NEG  = 1 
CONTINUE 

IF  (NZ  .LT.  N I VP  1 ) GC  TO  190 

IF  (.NOT. ID  .OR.  V (I.NVAR)  .Ey.  BLANK  .OB. 

1 V ( I , NV  A R)  .EC.  ABLANK)  GC  TC  2 10 
IF  (NZ  .GT.  0 .AND.  KEEPZ  .EQ.  C)  GO  TC  170 
CONTINUE 
GO  TO  260 

SET  N EQUAL  TC  NUMBER  OF  DATA  PCINTS. 

N = I - 1 

IF  (NEG  .EC.  1)  GO  TC  300 
AN  = N 
NEXV  = 0 
NVT  = 0 

DO  212  J = 1 , NIVP1 

IF  (EXV(J)  .NE.  0 . DC  .AND.  EXV(J)  .NE.  1.D0)  GO  TO  214 

CONTINUE 

GO  TO  216 

NVT  = 1 

CALL  THANS(V) 

IF  (ABORT)  RETURN 
TOTAL  DEGREES  OF  FREEDOM 
DFT  = AN  - 1 . CO  ♦ IA 
DEGREES  OF  FREEDOM  FOR  ERROR 
DPI  = N - NP  ♦ IA 

DEGREES  OF  FREEDOM  DUE  TO  REGRESSION 
DF2  = DFT  - DF 1 
IF  ( DF1  .LT.  0. DO ) GC  TO  280 
IF  (IORD  .EQ.  0)  RETURN 

ORDER  THE  DATA  FROM  LOW  TO  HIGH  VALUES  OF  Y. 

NK  = N - 1 
DO  250  I = 1,  NK 
IP1  = I ♦ 1 


INPT06 10 
INFT062C 
INPT0630 
INPT0640 
INPT0650 
INFT066C 
INPT0670 
INPT0680 
INPT0690 
INPT0700 
INFT07  1 0 
INPT07  20 
INPT0730 
INPT074C 
INPT0750 
INPT0760 
I NPT0770 
INPT0780 
INPT0790 
IN  PT0800 
INPT08 10 
INPT0820 
INPT0830 
INPT0840 
IN  PTO  850 
INPT0860 
INPT0870 
INPT0880 
INPT0890 
INPT0900 
INPT091C 
INPT0920 
INPT09  30 
INPT0940 
IN  PT0950 
INPT0960 
INPT0970 
INPT0980 
INPT099C 
I N PT 1000 
INPT10 10 
I N PT 1 020 
IN  FT  1 030 
INPT1040 
INPT1050 
IN  PT 1 06  0 
IN  PT 1 070 
INPT  1080 
I N PT 1 090 
IN  PT 1 1 00 
INPT11 10 
INPT1120 
IN  PT 1 1 30 
IN  PT 1 1 40 
INPT 1 1 50 
INPT1160 
IN  PT 1 1 7 0 
I N PT 1 1 80 
INPI1 190 
INPT 1200 
IN  PT 1 2 1 0 
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do  240  J = IP1,  N 

IP  ( V (1 , 1)  .LE.  V(J,1)|  GO  TO  240 

DO  230  K = 1,  NVAF 

TEMP  = V (J,K) 

L = J 

220  V(L,K)  = V(L-1,K) 

L * L - 1 

IP  (L  .Gl.  I)  GO  TO  220 
230  V(I,K)  = TEMP 
240  CONTINUE 
250  CONTINUE 
RETURN 

C ERROR  MESSAGES 

260  WRITE  (6,  270)  NMAX 

270  FORMAT  (//• ONUMB  EF  OF  INPUT  DATA  POINTS  HAS  EXCEEDED  ', 

1 'MAXIMUM  ALLOWABLE  (',13,').  THIS  JOE  HAS  bEEN  TERMINATED.') 
STOP 

2R0  WRITE  (6,  29C)  N 

290  FORMAT  (//'OTHE  NUMBER  OF  DATA  FCINTS  (’,11.  ')  IS  LESS  *, 

1 'THAN  THE  NUMBER  OF  PARAMETERS  TO  BE  SOLVED.  THIS  BUN  HAS  ', 

2 'BEEN  TERMINATED. • ) 

ABORT  = .TRUE. 

RETURN 

300  WRITE  (6,  310)  I 

310  FORMAT  (//'OA  NEGATIVE  VALUE  EXISTS  IN  THE  INFUT  DATA', 

1 ' (FOR  EXAMPLE,  DATA  CARD  NUMBER  ',  13,').  THIS  RUN  HAS  BEEN 

2 'TERMINATED.'  ) 

ABORT  = .TRUE. 

RETURN 

END 


IN  PT 1 220 
INFT1230 
INFT1240 
INPT1250 
I N PT 1 260 
INPT127C 
I N PT 1 280 
INPT 1290 
INFT1300 
IN  PT 1 3 1 0 
INPT1320 
I N PT 1 330 
INPT134C 
INPT1350 
INPT  1360 
INPT 1 370 
INPT1 380 
INPT 1 390 
INFT1400 
INPT 1 4 1 0 
INPT 1420 
INFT1430 
INFT1440 
I N FT  1 4 50 
INPT 1460 
INFT1470 
INPT1480 
INPT 1490 
I N PT 1 500 
IN  PT  1 5 1 0 
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. NE.  EXV  (J)  ( GC  TO  16C 


SUBROUTINE  TRANS(V) 

IMPLICIT  REAL*fl (A-H.C-Z) 

DIMENSION  V (NM  AX , KCLUMN ) 

COMMON  /Cl/  N,  NHAX,  KCLUMN,  NIV,  NIVF1,  Ml,  KF 
COMMON  /C2/  IEC 
COMMON  /C3/  ABORT 
COMMON  /CX/  EXV (8) , NEXV,  NVT 
LOGICAL*!*  AfaON  T 
) DO  140  J = 1 , NIVP1 

IF  (EXV  ( J)  . EC.  0.D0)  GO  TO  140 
IF  (EXV  (JJ  • EC',  08.  LO)  GO  TO  30 
IF  (EXV  (J)  .EC.  99. DC)  GO  TO  150 

IF  ( EXV  ( J)  .IE.  (-10. LO)  .OR.  EXV(J)  • GI . 47. DO)  GO  TO  160 
IF  (EXV  (J)  .GE.  10.10)  GO  TO  50 
EXPONENTIAL  TRAN SFCR M ATIO NS 
DO  20  I = 1,  N 
) V (I,  J)  = V (I  , J)  * * E XV  (J) 

GO  TO  140 

LOGARITHMIC  T RAN SFCR M ATIO NS 
) CO  40  I - 1,  N 
) V (I,  J)  = DLCG  (V  (I,  J)  ) 

GO  TO  140 
) KEXV  = EXV (J) 

IF  (FLOAT(KEXV)  .NE.  EXV(J))  GC  TO  16C 
NEXV  = 1 
KLV  = KEXV  / 10 
KMV  = KEXV  - 10  * KLV  ♦ 1 
IF  (KMV  .GT.  8)  GO  TC  160 
COMBINATIONS  CF  VARIABLES 
GO  TO  (60,  8C,  100,  120),  KLV 
) DO  70  I = 1,  N 
J V(I,J)  = V(I,J)  * V (I  , KM  V ) 

GO  TO  140 
) DO  90  I = 1,  N 
) V(I.J)  = V (L,J)  / V (I, KMV) 

GO  TO  140 
) DO  1 1 0 I = 1 , N 
) V(I,J)  = V(I,J)  ♦ V ( I , KMV) 

GO  TO  140 
) DO  130  I = 1,  N 
) V(I,J)  = V (I  , 0)  - V ( I , KM  V) 

) CONTINUE 
NEXV  = 0 
RETURN 

) IE  (J  .LE.  2 .OR.  NEXV  .EC.  0)  GC  TC  160 
NEXV  = NIVP1  ♦ 1 - J 
NIV  = NIV  - NEXV 
NIVP1  = NIV  ♦ 1 

IF  (I EC  .NE.  2 .AND.  IEQ  .NE.  4)  NP  = NP  - NEXV 
RETURN 

WRITE  ERROR  MESSAGE. 

) WRITE  (6,  170) 

) FORMAT  (//’ OA  TRANSFORMATION  FACTOR  HAS  NOT  BEEN  ENTERED  •, 
1 ’CORRECTLY.  THIS  RUN  HAS  BEEN  TERMINATED.’  ) 

ABORT  = .TRUE. 

RETURN 

END 


- V ( I , KM  V) 


TRANOO 1 0 
TPAN0020 
TRAN0030 
TBAN0040 
TRAN0050 
TRAN0060 
TRAN0070 
TRANOO0O 
TRAN0090 
TRAN0100 
TRAN 0 1 10 
TBAN0120 
TRAN013C 
TRAN01 40 
TR  AND  1 50 
TRAN0160 
TRAN017C 
TRAN0180 
TRAN0190 
TRAN020C 
TRAN0210 
TRAN0220 
TRAN0230 
TRAN0240 
TRAN0250 
TR  A NO  26  0 
TRAN0270 
TRAN028C 
TRAN0290 
TR  ANO  300 
TRAN0310 
TRAN0320 
TRAN0330 
TRAN0340 
TRAN0350 
TRAN0360 
TRAN0370 
TRAN0380 
TRAN0390 
TRAN0400 
TRAN04 10 
TRAN0420 
TBAN0430 
TRAN0440 
TR AN0450 
TRAN0460 
TRAN0470 
TRAN0480 
TRAN0490 
TRAN0500 
TRAN0510 
TRAN0520 
TRAN0530 
TRAN0540 
TRAN0550 
TRAN0560 
TRAN0570 


n n 
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SOBROUTINE  PRINT 
IMPLICIT  REAL*8 (A-H.C-Z) 

COMMON  /Cl/  HI (3) , NIV 
CONHON  /C2/  IEQ 

DATA  IE/'  •/,  IP/'**/.  IR/')'/ 

SUBROUTINE  FOB  PRINTING  SUBHEADINGS 


10 

20 


30 

40 

50 

60 


70 

80 

90 

100 


110 

120 


130 

140 


150 

160 


170 


GO  TO 
WRITE 
FORMAT 
C 

! ' ♦ G 
RETURN 
WRITE 
FORMAT 
RETURN 
WRITE 
FORMAT 

A1  , ■ * 
RETURN 
WRITE 
FORMAT 
RETURN 
WRITE 
FORMAT 

• c * 
: • g * 

BETURN 
WRITE 
FORMAT 
I A1, •♦ 
! A 1 , ' ♦ 
BETURN 
WRITE 
FORMAT 

•♦  C 
! •♦  G 
RETURN 
WRITE 
FORMAT 
I A1 , 

! A1 , • ♦ 
RETURN 
END 


(10,  3C , 50,  7C , 90,  1 10,  130,  150,  170),  IEQ 
(6,  20)  (IE,  1=1,  NIV) 

( ' OLIN  EAR  REGRESSION  — Y - A'.AI,'*  B * Xl',A1, 


* X2  ' , A 1 , ' ♦ D ♦ X3 ' , A1  , ' ♦ E 

* X6' , A 1, • ♦ H ♦ X7*  ) 


* X4 1 , A 1 , 1 ♦ F * X5 ' , A 1 , 


(6,  40) 

('OQUADRATIC  REGRESSION  --Y=A*B*X1*C*  X1»*2') 
(6,  60)  (IB,  1 = 1,  NIV) 

('OPOWER  REGRESSION  --  Y = A'.AI,'*  X1**B',A1,'*  X2**C', 
X3**D',A1,'*  X4**E ' , A 1 , ' • X5**F',A1,'*  16**G',A1,'»  X7**H 


(6,  80) 

{'0ASYMPTCTIC-POWER  REGRESSION  -- 


Y = A ♦ B * X1**C') 


(6,  100)  (IE,  IP,  I = 1,  NIV)  , IR 
( '0EXPCNENT1AI  REGRESSION  — Y = EXF(A',2A1,'  B * X1',2A1 
X 2 ' , 2 A 1 , ■ D * X3',2A1,'  E * X4',2A1,'  F * X5',2A1, 

X6»  ,2A1,'  H * X7  • , A 1) 

(6,  120)  (IB,  1=1,  NIV) 

(•CLOG-LINEAR  REGRESSION  — IN  Y = LN  A'.AI,'*  B * LN  XI* 
C * LN  X2  * , A1 , * ♦ C * LN  X3',A1,'*  E * LN  X4', 

F * LN  X5 ' , A 1 , ' ♦ G * LN  X6',A1,'*  H 


(6.  140)  (IE,  1=1,  NIV) 

(•OSEMILOG  REGRESSION  --  IN  Y = A'.AI,  B 

* X2',A1,**  C * X3 ’ , A 1 , ' ♦ E * X4,,A1,1 ♦ F * X5',Al, 

* X6  • , A 1 , ' * H » X7  • ) 


IN  X 7 • ) 

* XI  • , A 1, 


(6,  160)  (IE,  1=1,  NIV) 

(•OSEMILOG  REGRESSION  --  Y = A'.AI, •+  E * LN  XI', 
C * LN  X2 ' , A1 , • ♦ C * LN  X3',A1,'*  E * LN  X4', 

F * LN  X5',A1,'*  G * LN  X6',A1,'*  H * LN  XT) 


PRNT0010 
PHNT002C 
PRNT0030 
PRNT0040 
PRNT0050 
PRNT0060 
PRNT0070 
PRNT0080 
PRNT0090 
PRNT0100 
PRNT01 10 
PRNT0120 
PBNT01 30 
PRNT014C 
PRNT0150 
PPNT0160 
PRNT017C 
PRNT01 80 
PRNT0190 
1 ) PRNT0200 
PRNT02 1 0 
PRNT0220 
PRNT0230 
PRNT0240 
P8NT02S0 
, PRNT0260 
PRNT0270 
PRNT0280 
PRNT0290 
PRNT 0300 
, PRNT0310 
PRNT0320 
PRNT03  30 
PRNT034C 
PRNT0350 
PRNT0360 
PRNT0370 
PRNT0380 
PRNT0390 
PRNT04  00 
PRNT0410 
PRNT0420 
PRNT0430 
PRNT0440 
PRNT0450 
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SUEROUTINE  SUMS { V ) 

IMPLICIT  REAL«8  (A-h,C-Z) 

DIMENSION  V (NMAX.KCLUMN) 

COMMON  /Cl/  N,  NMAX,  KOIUF.N,  NIV,  NIVE1,  NVAR,  M1(2),  ID 
COMMON  /C2/  IEQ 

COMMON  /C3/  M2,  LINEAR,  IA,  K3,  NXSLT,  NYCOL,  M4  , NYUEV 
COMMON  /C  4/  S { 8)  , SYX(8,8) 

COMMON  /C5/  AN,  PB(4C),  VMSQ(8) 

COMMON  /C  6/  VMEA  N (8)  , T(fl),  H(8,8),  Z1(8),  SCfV(8) 

COMMON  /C7/  YCEPT 
COMMON  /C8/  EFT,  Z2(7),  SST 
LOGICAL  * 4 LINEAR,  IA 
C 

C SUBROUTINE  FCR  COMPUTING  VARIOUS  SUMS,  MEANS,  AND  STANDARD 

C DEVIATIONS  CF  THE  INFUT  EAT  I 

C 

LINEAR  = IEQ  .LE.  1 .OR.  IEQ  .GE.  8 
IF  (IEQ  .NE.  2)  GC  TC  20 
C MODIFICATIONS  FOR  QUADRATIC  CASE 

NIV  = 2 
NIVP1  = 3 
NVAR  = 3 ♦ IE 
DO  10  I = 1,  N 

IF  (IE  .GT.  0)  V (1,4)  = V (I,  3) 

JO  V (1,3)  = V (1,2)  * V (1,2) 

C INITIALIZE  SUMS. 

20  DO  30  J = 1 , NIV  P 1 
S(J)  = 0 . DO 
DO  30  N = J,  NIV  F 1 
30  SYX  (J,K)  = 0 . DO 

IF  (IEv  -LE.  2)  GO  TO  50 
J 1 = J 
J2  = NI VP  1 

IF  (IEQ  .EQ.  5 .CR.  IEQ  . EQ.  7)  J2  = 1 
IF  (IEQ  .EQ.  4 .CR.  IEQ  . EQ.  8)  J1  = 2 
DC  40  J = J1,  J2 
DO  40  I = 1,  N 

40  V ( I , J*  NY  DEV ) = DLCG (V  ( I , J ) ) 

50  L'O  80  I = 1,  N 
Y I = V (I  ,NYCOL) 

S(1)  = S(1)  ♦ YI 

SYX  (1  ,1)  = SYX  (1  , 1)  ♦ YI  * YI 

DO  70  J = 2,  NIVP1 

XIJ  = V(I,J+NXSET) 

S(J)  = S(J)  ♦ XIJ 

SYX  ( 1 ,J)  = SYX  (1  , J)  + YI  * XIJ 

DO  60  K = J,  NIVF1 

SYX  (J,K)  = SYX  (J  ,K ) ♦ XIJ  * V (I  ,K*  NXSET) 

60  CONTINUE 
70  CONTINUE 
80  CONTINUE 

C TOTAL  SUM  OF  SQUARES 

SST  = SY  X ( 1 , 1 ) - S ( 1 ) * S(1)  / AN 
C MEANS  AND  STANDARD  DEVIATIONS  OF  THE  INPUT  DATA 

DO  90  J = 1,  NIVP1 
VMEA  N (J)  = S < J)  / AN 

VMSQ(J)  = SYX  (J, J ) - S(J)  * S(J)  / AN 
90  SDEV(J)  = DSQRT  (VHSQ  (J)  / (AN  - 1 . DO)  ) 

IF  (IEQ  .EQ.  4)  RETURN 
IF  (IA)  GO  TC  120 


SUMS0010 
SUHS0020 
SUMS0030 
SUMS0040 
SUMS0050 
SUMS0060 
SUMS0070 
SUMS0080 
SUMS009C 
SUMSO  100 
SUMS01 10 
SUMS01 20 
SUMS01 30 
SUMSO 140 
SUMS0150 
SUHS016C 
SUMS017C 
SUMSO 180 
SUMSO 190 
SUMS0200 
S0MS02  10 
SUMS0220 
SUHS0230 
SUHS0240 
SUMS0250 
SUMS026C 
SUMS0270 
SUMS0280 
SUMS0290 
SUMS0300 
SUMS0310 
SO  H SO  320 
S0MS0330 
SUMS0340 
SUMSO  350 
SUMS0360 
SUMS037C 
SUMS0380 
SUMS0390 
SUMS0400 
SUMS04  10 
SUBS0420 
SUHS04  30 
SUHS0440 
SUNS0450 
SUHS0460 
SUBS0470 
SUBS0480 
SUBS0490 
SUBS0500 
SUBS0510 
SUBS0520 
SUBS05  30 
SDBS0540 
SUBS0550 
S0BS0560 
SUBS0570 
SUBS0580 
SUBS0590 
SUBS0600 
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DO  11C  J = 2,  NIVP1 

SUHS06 1 C 

1 (J~  1 ) = SYX(1,J)  - S(1)  » S(J)  / AN 

SUNS0620 

DO  100  K = J,  NIVP1 

SUMS06  30 

100 

H(J-1,K-1)  = SYX(J,K)  - S (J)  ♦ S(K)  / AN 

SUHS0640 

110 

CONTINUE 

SUMS065C 

RETURN 

SUMS0660 

120 

DO  140  J = 2,  NIVP1 

S0MS067C 

T(J-1)  = SYX(1,J)  - YCEPT  * S(J) 

SUKSC60C 

DO  130  K = J,  NIVP1 

SUMS0690 

130 

H (J-  1,K-  1)  - SYX  (J,K) 

SUHS0700 

140 

CONTINUE 

5UBS07  10 

RETURN 

SUMS072C 

FND 

SUBS0730 

SUBROUTINE  LINE { V) 

IMPLICIT  REAL*8 (A-H,C-Z) 

DIMENSION  V (NHAX ,KCLDMN) 

COMMON  /Cl/  N,  NHAX,  KOLUMN,  NIV,  N I V F 1 , Ml,  NP 
COMMON  /C2/  IEC,  H2(9),  IND 

COMMON  /C  3/  ABORT,  M3,  IA,  M4 , NXSET,  M5,  NYC 
COMMON  /C5/  Z1,  P (8) 

COMMON  /C6/  VHEAN(8),  T(8) 

COMMON  /C7/  YCEPT,  Z2{4),  XV,  YV 
LOGICAL*4  ABORT,  IA 
DATA  INDI/'LINE1/ 

SUBROUTINE  FCR  DETERMINING  L E AST -SyU AH ES  SOLUTIONS  OF  PARAMETERS 
FOR  LINEAR  FUNCTIONS  OF  FORM 

Y = A ♦ B*  X 1 ♦ C *X  2 + C*X3  ♦ E*X4  ♦ F*XS  + G*X6  ♦ H*X7 

NHERE  A MAY  EE  SPECIFIED 

IND  = IN  D 1 

SOLVE  FOR  PARAMETERS  (P)  . 

CALL  SOLVE (NIV) 

IF  (ABORT)  RETURN 
P (1)  = VMEAN  (1) 

DO  10  J = 2,  NP 
P(J)  = T(J-1) 

10  P < 1 ) = P(1)  - P(J)  * VMEAN(J) 

IP  (IA)  P ( 1 ) = YCEPT 

Y-COMPUTED  AND  Y-RESIDUAL  VALUES 
DO  30  I = 1,  N 
V (I  , NYC)  = P(1) 

DO  20  J = 2,  NIV F 1 

20  V (I , NYC)  = V (I,  NYC)  ♦ P(J)  * V(I,J*NXSET) 

30  CONTINUE 

IF  (IEQ  .NE.  2)  RETURN 
XV  = -P  (2)  / (2.  DO  * F (3)  ) 

YV  = P ( 1 ) ♦ P (2)  * XV  ♦ p (3)  * XV  * XV 

NIV  = 1 

NIVP1  = 2 

RETURN 

END 


LINE0010 
LINE0020 
LINE0030 
LINE004C 
LINE0050 
LINE0060 
LINE0070 
LINEC08C 
LINE0090 
LINE0100 
LINE01 10 
LINE0120 
LINEO  130 
LINE0140 
LINEO 1 50 
LINE0160 
LINE0170 
LI  N EC  1 8 0 
LI NEO 1 90 
LI N EO  200 
LINE0210 
LINE022C 
LINE0230 
LINE0240 
LINE0250 
LINE0260 
LINE0270 
LI N EO  28  0 
LINE0290 
LINE0300 
LINE0310 
LINE0320 
LINE0330 
LINE0340 
LINE0350 
LINE0360 
LINEC37C 
LINE0380 
LINE0390 
LINE0400 
LINE041C 


LP-P'i  ' 


■! 


; 


: 


[ 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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SUBROUTINE  EXPO(V) 

IMPLICIT  HEAL*8  (A-H.C-Z) 

DIMENSION  V (NMAX.KCLUMN) 

COMMON  /Cl/  Ml.  NMAX  , KOLUMN,  NIV,  M 2 (2)  , NP 

COMMON  /C2/  1EQ,  M3(«),  IND 

COMMON  /C3/  ABORT,  M4(2),  IGUESS 

COMMON  /Lb/  VMEAN(8),  T(8),  Z1(64),  FI  (8) 

LOGICAL*  4 AEORT,  IGUESS 

IATA  I ND2/ • EXFG •/ 


SUBROUTINE  FOR  DETERMINING  LE A ST -SQU A RS3  SOLUTIONS  OP  PARAMETERS 
POR  POWER  FUNCTIONS  CP  FORM 


Y = A » X 1 **E  * X2**C  ♦ X 3 **  C * X4»*E  * X5»*F  * X6**G  * X7**H 
OR  FOR  EXPONENTIAL  FUNCTIONS  OF  FORM 

Y = EX  F ( A ♦ B*  X 1 ♦ C*X2  ♦ C*X3  4 E»X4  ♦ F*XS  ♦ G*X6  ♦ H*X7) 


SET  SUBROUTINE  INDICATOR 
IND  = IN  C2 

IF  (IGUESS)  GC  TO  20 

FIRST,  DETERMINE  LE A ST -SQUAR E S SOLUTICNS  OF  PARAMTEBS  PI  (J) 

FIR  LOG-LINEAR  FORM  (FONEH) 

J.N(Y)  = LN(A1)  ♦ El  ♦ LN(X1)  4 Cl  * LN(X2)  4 ...  4 Hi  * LN(X7) 
OR  SEMILCG-LINEAB  FORM  (EXPONENTIAL) 

L N ( Y)  = A1  4 B1  * XI  4 Cl  * X2  4 . . . 4 HI  * X7 
CALL  SOLVE  (NIV) 


IF  (ABORT)  RETURN 
P1(1)  = VHEAN(I) 

DO  10  J = 2,  NP 
PI  (J)  = T(J-1) 

10  PI  (1)  = PI  (1)  - PI  (J) 


VHEAN ( J ) 

IP  ( I EQ  .EQ.  3)  F 1 ( 1 ) = D EXP  (PI  (1)  ) 

DETERMINE  LEAST-SQUARES  SOLUTIONS  OF  PARAMETERS  P(J) 


20  CALL  ITER(V) 
RETURN 
END 


EXPO00  10 
EXP00020 
EXF0003C 
EXP00040 
EXPO0050 
EXP00060 
EXFOC07C 
EXP00080 
EXP00090 
EXPOOIOC 
EXP001 10 
EXP00120 
EXPOOI 30 
EX  POO  1 4 C 
EX  FOO 1 50 
EX  PO0 1 6 0 
EXPOC17C 
EX  FO0 1 80 
EXPOO 190 
EIP00200 
EXFO021C 
EXPC0220 
EXPO0230 
EXF00240 
EXP00250 
EXPC0260 
EXP00270 
EXP00280 
EX  FOO  290 
EXP00300 
EIPO0310 
EXP00320 
EXPC033C 
EXF00340 
EXPC0350 
EXP00360 
EXPOC370 
EX  FOO  3 80 


SUBROUTINE  ASYM(V) 

IMPLICIT  REAL*8 (A-H.C-Z) 

DIMENSION  V (NMAX.KCLUBN) 

COMMON  /Cl/  N,  NMAX,  KCLUMN,  M 1 < 2)  , NVAR 

COMMON  /C3/  ABORT,  M2,  IA,  M3  (3),  NYC 

COMMON  /C4/  S (8)  , SYX11 

COMMON  /C5/  AN,  A,  B,  C 

COMMON  /C6/  VMEAN(8)  , T(8),  H ( 8 , 8) 

COMMON  /C7/  YCEPT , DELTA 

COMMON  /C8/  Z 1 (6  ) , YDEVSQ,  Z2,  SST 

DIMENSION  SUM (9),  SAVSUM(4) 

LOGICAL*4  ABORT,  IA,  SOLVED 

SUBROUTINE  FOR  DETERMINING  L EAST -SOU AR E S SOLUTIONS  OF  PARAMETERS 
FOR  ASYMFTOT IC-PCWER  FUNCTIONS  CF  FCRM 


WHERE  A MAY  EE  SPECIFIED 

A = YCEPT 
SOLVED  = .FALSE. 

SET  INITIAL  VALUE  CF  Cl 
Cl  = -8.  10 1 DO 

SET  FIRST-ITERATION  LESIG NATCH  TO  1. 
ITERS  = 1 

SET  C INCREMENT  INITIALLY  TO  0.1. 

DC  = 0. IDO 
DO  200  K = 1 , 162 

STEP  INITIAL  C VALUE  BY  DC  INCREMENT. 
Cl  = Cl  ♦ DC 
IF  (K  .NE.  82)  GC  TO  1C 
ITERS  = 1 
Cl  = 0.0C1D0 
10  C = Cl 
20  DO  30  I = 1,  1 
30  SUM ( I ) = 0 . DO 
DO  40  I = 1,  N 
XP  = V (1,2)  **C 
XPSQ  = XP  * XF 
V (I,NVAfl+4)  IS  LN  (XI)  . 


XP 1 = XP 
SUM  (1)  = SUM  (1)  ♦ XF 

SUM  (2)  = SUM  (2)  ♦ XPSC 

SUM  (3)  = SUM  ( 3)  ♦ Xfl 

SUM  (4)  = SUM  (4)  ♦ XP  SO  * V(I,NVAR*4) 

SUM  (5)  = SUM  (5)  ♦ XF  1 * V (I,  1) 

SUM  (6)  = SUM  (t)  ♦ XP  * ( V (1 , 1 ) - VMEAN(I)) 

IF  (IA)  SUM  (7)  = SUM  (7)  ♦ XP  * V(I,1) 

40  CONTINUE 

IF  (.  NOT.  IA)  GO  TO  5C 
B = (SUM  (7)  - A * SU  F ( 1 ) ) / SUH(2) 

G = SUM  (4)  *SUH  (7) -SUM  (2)  *SUM  (5) -A*  (SUM  (1)  *SUN  (4)  -SUM  (2)  *SUH(3)  ) 
YDEVS1  = SYX11  - 2 . DC*  A*S  (1)  ♦ AN*  A*  A - E*B*SUM(2) 

GO  TO  60 

50  B = SUM  (6)  / (SUM  (2)  - (SUH(1)  * SUM(1)/AN)) 

A = (S  (1)  - (E  * SUM  (1)  ) ) /AN 
G = SUM  (5)  - (B  * SU  M (4)  ) - (A  ♦ SUH(3)) 

YDEVS1  = SST  - B ♦ SUM<6) 

60  IP  (K  .EQ.  1)  QSAVE1  = YDEVS1 


V (I,NVAR*4) 


ASYMOO  10 
ASYM002C 
ASYH003C 
ASYM004C 
ASYM0050 
ASYM0060 
ASYMC07C 
AS  YM008  0 
ASYM0090 
ASYM0100 
ASYMO 1 10 
AS  Y M0 1 20 
ASYMO 1 30 
ASYM0140 
ASYMC150 
ASYM0160 
ASYMO 170 
ASYM0180 
ASYM0190 
ASYH0200 
ASYM0210 
ASYN0220 
AS  Y M02  30 
ASYM0240 
ASYM0250 
ASYM0260 
ASYM0270 
ASYM0280 
ASYM0290 
ASYM0300 
ASYM03 10 
ASYM0320 
ASYM0330 
ASYMC340 
ASYM0350 
ASYM0360 
ASYM0370 
ASYM0380 
ASYM0390 
ASYM0400 
ASYM04 10 
ASYH0420 
ASYM0430 
ASYM0440 
ASYM0450 
ASYH046C 
ASYH0470 
ASYM0480 
ASTM0490 
ASYM0500 
ASYB05  10 
ASYM0520 
ASIH0S30 
ASIM0540 
ASYB0550 
ASYB0560 
ASYB0570 
ASYB0580 
ASYBC59C 
ASYB060C 


n n 
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l 

1 

j i 


« 

i 

i 

l 


IF  (K  .Eg.  162)  QSAVE3  = YDEVS1 

ASYM06  10 

IF  (ITEhS  .Eg.  2)  GC  TC  100 

ASYM0620 

IF  (G)  7C.16C.80 

ASYM0630 

70 

H = -1 

ASYM0640 

GO  TO  90 

ASYM0650 

8C 

M = 1 

ASYN066C 

90 

ITERS  = 2 

ASYM0670 

GO  TO  190 

ASYM0680 

100 

IF  (M  .GT.  0)  GO  TC  110 

ASYM0690 

IF  (G)  130, T6C, 120 

ASYM070C 

110 

IF  (G)  120.16C.130 

ASYM07  10 

120 

C = C - {DC  * 0.5DC) 

ASYM0720 

GO  TO  140 

ASYM0730 

130 

IF  (DC  .GT.  C.07I0)  GC  TO  190 

ASYM0740 

C = C ♦ (DC  » C.5DC) 

ASYN0750 

140 

IF  (IA)  GO  TC  142 

ASYM0760 

DDA  = DABS (A/ASTORE  - 1.D0) 

ASYM077C 

IF  (DDA  .GE.  DELTA)  GO  TO  150 

ASYM0780 

142 

DDE  = DABS (S/BSTCRE  - 1.D0) 

ASYN079C 

IF  (DDE  .GE.  DELTA)  GO  TO  150 

ASYM0800 

DDC  = DABS (C/CSTORE  - 1.D0) 

ASYHOP 1 C 

IF  (DDC  .GE.  DELTA)  GO  TO  150 

ASYM0820 

GO  TO  160 

ASYN0830 

150 

ASrORE  = A 

ASYMC84C 

BSTORE  = B 

ASYM0850 

CSTORE  = C 

ASYM086C 

CC  = IC  * 0.5D0 

ASYM0870 

GO  TO  20 

ASYMC88C 

USE  NEC  VARIAELES  FCR  TEMIORARY  SOLUTION. 

ASYM0890 

160 

YDEVSQ  = YDEVS1 

ASYM0900 

IF  (.NOT.  SOLVED  .OR.  YDEVSg  .LE.  CSAVE2) 

GO 

TO  170 

ASYH0910 

A = ASAVE 

ASYM0920 

E = BSAVE 

ASYM09  30 

C = CSAVE 

ASYM0940 

SUM(1)  = SAVSUM(I) 

ASYM095C 

SUM (2)  * SAVSUM (2) 

ASYM0960 

SUM  (4)  = SAVSUM  (4) 

ASYM0970 

YDEVSg  = 0SAVE2 

ASYM0980 

GO  TO  180 

ASYH099C 

STORE  PARAMETER  VALUES  AND  SUM  CF  SQUARES 

CF 

Y RESIDUALS. 

AS  Y Ml  000 

170 

ASAVE  = A 

ASYM10  10 

ESAVE  = E 

ASYM1020 

CSAVE  = C 

ASYB1030 

SAVSUM ( 1 ) = S UM ( 1 ) 

ASYM1040 

SAVSUN(2)  = SUM  (2) 

ASY  B 1050 

SAVSUM(4)  = SUM (4) 

AS  YM 106  0 

QSAVE2  = YDEVSQ 

ASYB107C 

SET  FIRST-ITERATION  DESIGNATOR  TO  1. 

ASYB1G80 

100 

ITERS  = 1 

ASYB 1090 

SET  SOLUTION  DESIGNATOR  TC  TRUE 

ASYB1 100 

SOLVED  = .TRUE. 

ASYB1 1 10 

DC  = 0.  IDO 

ASYB 1 120 

Cl  = Cl  - DC 

ASYB1130 

GO  TO  200 

ASYB  1 140 

190 

ASTORE  = A 

ASYB1150 

BSTORE  = B 

ASYB  1 160 

CSTORE  = C 

ASYB1170 

200 

CONTINUE 

ASYB118C 

IF  A UNIQUE  SOLUTION  FOB  C DOES  NOT  EXIST 

IN 

THE  SPECIFIED  RANGE, 

ASYB1190 

PRINT  A MESSAGE  RELATING  TO  THAT  FACT. 

ASYM1 200 

IF  (SOLVED. ANC.YCEVSC.LT. QSAVE1. AND. YDEVSQ 

• LT 

• USAVE3)  GO  TO  220 

ASYB1210 

1 
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* 


WRITE  (b,  210) 

ASYB1220 

210 

rORNAT  (//'UNO  SOLUTION 

HAS 

PEEK  ECUNC  ECR  THIS  PRCBLEfl  IN 

ASYB  123C 

1 

1 ’THE  PANGE  Cf  -8  TO  «8 

ECR 

c.  •) 

ASYB124C 

ABORT  * .TRUE. 

AS YB 1 2 50 

RETURN 

ASYN1260 

RESTORE  VARIABLES  TC  ORIGINALS. 

ASYB1270 

220 

A = ASAVE 

ASYB128C 

E = BSAVE 

ASYB1290 

C = CSAVE 

ASYB  1 300 

SUM  ( 1 ) = SAVSUB(I) 

ASYB1310 

SUB  ( 2)  = SAVSUB(2) 

ASYB1320 

SUB (4 ) = SAVSUB(4) 

ASYB13 30 

SUB  (8)  = 0.  DC 

ASYB 1 340 

SUB (8)  = 0. DO 

ASYB1350 

T(1)  = 0 . DO 

ASYH136C 

T (2)  = 0 . DO 

ASYB1 370 

T ( 3)  = 0 . DO 

ASYB 1 380 

DO  2 JO  I = 1,  N 

ASYB1390 

XP  = V (1 ,2)  *»C 

ASYB 1400 

VII. NYC)  = A ♦ B » XP 

ASYB1410 

XPL  = V(I,NVA6*4)  ♦ X F 

ASYB1420 

SUB  ( 8 ) = SUB  (8)  ♦ XPL 

ASYB143C 

SUB  (9)  = SUB  (9)  ♦ XFL  * 

XPL 

ASYB144C 

2 JO 

CONTINUE 

ASYB 1450 

IP  (IA)  GO  TC  240 

ASYB 1460 

H (1,  1)  = AN 

ASTH147C 

H (1,2)  = SUB  |1) 

ASYB 1480 

H ( 1 , 3 ) = B * SUB  (8) 

ASYB1490 

GO  TO  25C 

ASYB 1500 

240 

H (1,1)  = 1 . DO 

ASYB1510 

H (1,2)  = O.DC 

ASYB 1520 

H (1 ,3)  = 0. DO 

ASYB1530 

250 

H (2,2)  = SUB  (2) 

ASYB 1540 

H (2,3)  = B • SUB  (4) 

ASYB1550 

H (3,  3)  = 8 * B * SUB  (9) 

ASYB1560 

CALL  SOLVE (3) 

ASYB 1570 

RETURN 

ASYB1580 

END 

ASYB  1590 
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SUBHOUTINE  SOLVE (NSIZE) 

IMPLICIT  REAL»8  (A-H.C-Z) 

common  /c2/  ri  (10) , ind 

COMMON  /C  3/  ABORT 

COMMON  /C6/  Z 1 |S)  , T <S)  , h<8,8) 

LOGICAL*^  ABORT 

EIHENSION  IPIVCT(8),  INEEX(8,2) 

C 

C SUB80UTI NE  FCF  SOLVING  SIMULTANEOUS  EQUATIONS 

C H BECOMES  INVERTED  H;  T BECOMES  SOLUTION  VECTOR. 

C 

IF  (NSIZE  -GT.  1)  GC  TO  10 
IF  ( H ( 1 , 1)  -EC-  C.IC)  GC  10  ISO 
H (1  , 1)  = 1 . D 0 / 8(1,1) 

T(1)  - H ( 1 , 1 ) * 1(1) 

RETURN 

C FILL  OUT  LONER  TRIANGLE  CF  II  MA1HIX. 

10  CO  20  J - 2,  NSIZE 

JM1  = J - 1 
DO  20  K = 1,  JM1 
20  H(J,K)  = H (K  , J) 

C INITIALIZATION 

DC  30  J = 1,  NSIZE 
30  IPIVOT(J)  = C 

C PARTIAL  MATRIX  INVERSION  ROUTINE 

DO  1 20  1=1,  NSIZE 
C SEARCB  FOR  PIVOT  ELEMENT 

HMAX  = 0 . DO 
DO  60  J = 1,  NSIZE 
IF  (IPIVCT(J)  .EQ.  1)  GC  TO  60 
DO  50  K = 1,  NSIZE 
IF  (IPIVOT(K)  - 1)  4 C , 50,  150 
40  IF  (DABS  (HMAX)  .GE.  L A ES  ( H (J , K ) ) ) GC  TC  50 
IRON  = J 
ICOL  = K 
HMAX  = H (J,  K) 

50  CONTINUE 
60  CONTINUE 

IF  (HMAX  .EQ.  0. BO)  GC  TO  150 
I PI  V OT  ( ICOL)  = I PI  VCT  ( ICOL)  ♦ 1 
C INTERCHANGE  RONS  TO  PET  PIVOT  EIEMENT  CN  DIAGONAL 

IF  (IRON  .EQ.  ICCL)  C-C  TO  80 
DO  70  L = 1,  NSIZE 
TEMP  = H(IRON, L) 

H(IROW, L)  = H (ICC  L , L ) 

70  H (ICOL, L)  = TEMP 
TEMP  = T (IRON) 

T (I RON)  = T(ICOL) 

T (ICOL)  = TEMP 
80  PIVOT  = H( ICCL, ICOL) 

INDEX  (I,  1)  = IRON 
INDEX  (I,  2)  = ICOL 

C DIVIDE  PIVOT  RON  EX  FIVOT  ELEMENT 

H (ICOL,  ICOL)  = 1 • DO 
DO  90  L = 1,  NSIZE 
90  H(ICOL,L)  = H (ICOL,L)/PlVCT 
T ( ICOL)  = T (ICCL) /PIVOT 
C REDUCE  NCN-PIVOT  RONS 

DO  110  LI  = 1,  NSIZE 
IF  (Li  .EQ.  ICOL)  GO  TO  110 


SOLVOOIO 
SOLVOO  20 
SO  L V00  30 
S01V004C 
SOLV0050 
SOLVOO 60 
SOLV0070 
SOLVC08C 
SOLV0090 
SOLVO 100 
SOLVO  1 1 0 
SOLV01 20 
SOLVO 1 30 
SOLV0140 
SOLVO 150 
SO l VO  1 60 
SO  LV0 1 70 
SOLVO 1 8C 
SOL  V0 1 90 
SOLV0200 
SOLV02 10 
SOLV0220 
SOLV0230 
SOLV0240 
SOLV0250 
SCLV026C 
SOLVC270 
SOL VO 280 
SOLV0290 
SOLV0300 
SOLV03  10 
SO  L V 0 3 2 0 
SOLV033C 
SOLV0340 
SOL  VO 350 
SOLV0360 
SOL V 0 370 
SOLV0380 
SOL VO 390 
SOL 704 00 
SOLV041C 
SOLV04  20 
SOLV0430 
SOLV0440 
SOLV0450 
SOLV0460 
SOLV0470 
SOLV0480 
SOLV0490 
SOLV0500 
SOLV05  10 
SOL VOS 20 
SOLVQ530 
SOLV0540 
SOLV0550 
SOLV0560 
SOLV0570 
SOLV0580 
SOLV0590 
SOLV0600 


1 
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TEMP  = H(L1,IC0L) 

SOLV06 1 0 

H ( L 1 , ICO  L)  = C . DC 

SOLV062C 

DO  100  L = 1,  NSIZE 

SGLV0630 

100 

H(L1,l)  = H(L1,L)  - H(ICCL,L)  * TfcMF 

SOIV0640 

T(L1)  = T ( L 1 ) - T(ICCL)  * TEMP 

SO LV 06 50 

110 

CONTINUE 

SOLV0660 

120 

CONTINUE 

SOLV0670 

INTERCHANGE  COLUMN’S. 

SOLVO60O 

DO  1 4 C 1 = 1,  NSIZE 

SOLV069C 

L = NSIZE  +1-1 

SOLV0700 

If  (INDEX  (L  , 1)  .EO.  INDEX  (L, 2))  GO  TO  140 

SOLV07  1C 

IRON  = INDEX { L , 1 ) 

SOLV0720 

ICOL  = INDEX (L,2) 

SOLV0730 

DO  1 30  K = 1,  NSIZE 

SOL V 07 40 

TEMP  = H(K.IFCW) 

SOLV075G 

H (K , I SO W)  = H (K  , ICC I ) 

SOLV  076  0 

h (K, ICOL)  = TEMP 

SOLV0770 

130 

CONTINUE 

SOLV078C 

140 

CONTINUE 

SOLV07  90 

RETURN 

3OLV080C 

ERROR  MESSAGE 

SOLV08 10 

150 

WRITE  (6,  IOC)  IND 

SOL  V 08  20 

160 

FORMAT  ('OZERO  DETERMINANT  IN  SUBROUTINE  SOLVE.  THIS  RUN 

SOLV08  30 

1 1 HAS  BEEN  TERMINATED.  '/  1H0  , 'SOLVE  WAS  LAST  CALLED  FROM  A4) 

SOLV0840 

ABORT  = .TRUE. 

SOLV085C 

RETURN 

SOLV0860 

END 

SOLV087C 

. i 


I 
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SUBR0UI1NE  ITEB(V) 

IMPLICIT  REAL*8 (A-H.C-2) 

DIMENSION  V (NMAX.KCLUHN) 

COMMON  /Cl/  N,  NMAX,  KCLUMN,  HI,  NIVP1 
COMMON  /C2/  1EQ,  M 2 ( 8 ) , L1M,  INC 
CuKMON  /C  3/  ABORT,  M3  (3),  N X SET  , NYCCl,  NYC 
COMMON  /C9/  S<8),  SYX  (8,8) 

COMMON  /C5/  AN,  P (P) 

COMMON  /C6/  VMEAN  (8)  , T(8),  H(8,8),  PI  (8)  , SD£V(8) 

COMMON  /C7/  Z1,  LEITA,  22,  KA 

COMMON  /C8/  DPT,  23(7),  3 ST 

LOGICAL*  4 ABORT,  SCLVfD,  IEQ3 

DIMENSION  F P (0  ) , FTEMF  (3,8)  , C(J),  DIEF(R) 

DATA  I NC  3/ * IT  E P * / 

SUBROUTINE  FOR  DETERMINING  L E A ST-SQ  U AR  ES  SOLUTICNS  CF  PARAMETERS 
FOR  NCN-LINEAR  FUNCTIONS  WHERE  AN  ITERATIVE  PROCEDURE  IS 
REQUIRED  (VIZ,  POWER  (IEC  = 3)  AND  EXPONENTIAL  (IEC=^)  FUNCTIONS) 


SET  SUBROUTINE  INDICATOR 

IND  = 1ND3 

IEQ3  = I EQ  .EQ.  3 

SET  INITIAL  GUESSES 

DO  1 0 J = 1 , NIVP 1 

P(J)  = P1(J) 

RECALCULATE  SUM,  SUMSQ,  MEAN  ANE  STD  DEV  FOR  ACTUAL  (NON  LCG)  DATA. 

DO  20  J = 1,  NIVF1 

S(J)  = 0 . DO 

DO  20  K = J,  NIVP1 

SYX  (J,K)  = O.DO 

DO  30  I = 1,  N 

DO  30  J = 1,  NIVP 1 

S(J)  = S(J)  ♦ V(I,J) 

DO  30  K = J,  NIVPt 

SYX  (J  ,K)  = SYX(J,K)  ♦ V(i,J)  * V (I  , K ) 

DO  40  J = 1,  NIVF1 
VMEAN  (J)  = S(J)  / AN 

SDEV(J)  = DSQRT  ( (SYX  (J  ,J)  - S(J)  * S(J)  / AN)  / DET) 

SST  = SYX(1,1)  - S{1)  * S{1)  / AN 
SET  SOLUTION  INDICATOR  TO  FALSE 
SOLVED  = .FALSE. 

DO  31C  1=1,  LIM 

COMPUTED  Y VALUES  AND  PARTIAL  CE  Y WITH  RESPECT  TO  PARAMETER  A. 

IF  (I  EQ  3)  GC  TC  80 
DO  70  I = 1,  N 
TEMP  = P (1) 

DO  60  J = 2,  NIVP 1 

TEMP  = TEMP  ♦ P(J)  * V(I,J) 

V (I  , NYC)  = DEXP(TEMF) 

GO  TO  110 

DO  100  I = 1,  N 
TEMP  = P ( 1 ) 

DO  90  J = 2,  NIVP 1 

TEMP  = TEMP  * V(I,J)**P(J) 

V ( I , N YC)  = TEMP 

IF  A SOLUTION  HAS  BEEN  OETAINEE  (SOLVED=TRUE)  , GO  TO  ENDING. 

IF  (SOLVED)  GC  TO  330 
CLEAR  H AND  T MATRICES. 

DO  130  1=1,  NIVP 1 
T ( I)  = O.DO 


: 
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do  120  J = I,  NIVP1 
120  H(I,J)  = 0 .DO 
130  CONTINUE 

DO  170  I = 1,  N 

C PARTIAL  OF  Y PUNCTION  WITH  RESPECT  TC  PARAHETEPS 

FP(1)  = V(I.NYC) 

IF  (I  EQ 3)  FP  ( 1 ) = V (I,  NYC ) / P ( 1 ) 

DO  140  J = 2,  NIVP1 
140  FP  (J)  = V (I,  NYC)  * V(I,J»NXSET) 

C DIFFERENCE  BETWEEN  OBSERVED  Y AND  CALCULATED  Y (Y  RESIDUAL) 

YDIFF  = V(1, 1)  - V (I, NYC) 

C CALCULATE  H AND  T HATRICES. 

DO  160  II  = 1,  NIVP1 
DO  150  JJ  = II,  NIVP1 

150  H(II,JJ)  = H (II,  JJ)  ♦ (FP  (II)  * FP(JJ)) 

T ( II)  = T (II)  ♦ YDIFF  * FP  (II) 

160  CONTINUE 
170  CONTINUE 

C SOLVE  FOR  CORRECTIONS  TO  PREVIOUS  SOLUTIONS. 

CALL  SOLVE (NIVP1 ) 

IF  (ABORT)  RETURN 

C FIND  WHICH  FRACTIONAL  PART  OF  CORRECTION  TERHS,  WHEN  ADDED  TO 

C PARAMETER  VALUES,  GIVES  LOWEST  SUB  OF  SQUARES  OF  Y RESIDUALS. 

TEBP  = 1.0D0 
180  TEBP  = 0 . 6 DO  * TEMP 
DO  250  J = 1,  3 
FI  = TEBP  * (J  - 1) 

DO  190  !C  = 1 , NIVP1 
190  PTEBP  (J,  K)  = P(K)  ♦ T (K)  * FI 
Q (J)  = 0 • DO 

DO  240  1 = 1,* 

YTEHP  = PTEBP (J, 1) 

IF  (I EQ 3)  GO  TO  210 
DO  20 0 K = 2,  NIVP1 

200  T TEBP  = YTEBP  ♦ PTEBP(J,K)  * V ( I , K) 

YTEHP  = DEXP (YTEBP) 

GO  TO  230 

210  DO  220  K = 2,  NIVP1 

220  YTEHP  = YTEHP  * V ( I,  K)  **  PT  EH  P ( J , K) 

230  YDIFP  = V (I, 1)  - YTEHP 

Q(J)  = Q (J)  ♦ (YDIFF  * YDIFF) 

240  CONTINUE 
250  CONTINUE 
LB  * 1 

DO  260  J = 2,  3 

IF  <Q(LH)  .LT.  Q (J) ) GO  TO  260 

LB  = J 

260  CONTINUE 

IF  (LB  .GT.  1)  GO  TO  280 
DO  270  1 * 1,  NT  VP  1 

IF  (DABS  (T  (I)  * 2.  DO  * TEBP)  .GT.  DELTA)  GO  TO  180 
270  CONTINUE 
280  DO  290  1=1,  NIVP1 

DIFF(I)  * DABS  (PTEBP  (LH,I)/P  (I)  - 1 . DO ) 

P(I)  = PTEBP  (LB,  I) 

290  CONTINUE 

DO  300  1*1,  NIVP1 

C SOLVED  WHEN  RELATIVE  DIFFERENCE  IS  LESS  THAN  DELTA. 

IF  (DIFF(I)  .GE.  DELTA)  GO  TO  310 
300  CONTINUE 

SOLVED  * .TRUE. 


ITER06  10 
ITEP06  20 
ITER0630 
ITER7640 
ITER0650 
ITER0660 
ITER06  70 
ITER06  80 
ITER06  90 
ITER0700 
ITER07 1 0 
ITER072C 
ITER07  30 
ITER0740 
ITER0750 
ITER  07  60 
ITER0770 
ITER07  80 
ITER0790 
TTER0800 
ITERC8  10 
ITER08  20 
ITEP9830 
ITER0840 
ITER0850 
ITER0860 
ITERG870 
ITERC8R0 
ITSR08  90 
ITERC90D 
ITER09 1 0 
ITER0920 
ITER0930 
ITER0940 
ITER0950 
ITER0960 
ITER0970 
ITER0980 
ITER0990 
ITER  1000 
ITER10  10 
ITER  1 020 
ITER  10  30 
ITER  1040 
ITER  10  50 
ITER1 060 
IT BRIO 70 
ITER  1080 
ITER1090 
ITER  1100 
ITER111C 
ITER1120 
ITER  11 30 
ITER  1140 
ITER  11 50 
ITER1160 
ITER  1 1 70 
ITBR1180 
ITER  1190 
ITER1200 
ITER1210 
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c GO  BACK  AND  CONFUTE  NEW  YC  AND  YDEV  FOR  FINAL  CORRECTIONS 
310  CONTINUE 

IF  (SOLVED)  GO  TO  50 
C ERROR  RESSAGE 

WRITE  (6,  320)  LIN,  (J,  PI  (J)  , J,  P(J),  DIFP(J),  J = 1,  NIVP1) 
320  FORMAT  (//'ONO  SOLUTION  HAS  BEEN  OBTAINED  FOR  THIS  RUN 

1 'AFTER  •.  13,  • ITERATIONS.  THIS  RUN  HAS  BEEN  TER HIN AT  ED. • / 

2 1 HO , 71,  'INITIAL  GUESSES',  20X,  'SOLUTION  AT  TERMINATION', 

* 20X,  'RELATIVE  DIFFERENCE'  / (/ ' P1(',  II,  •)  =',  D15.5, 

<*  20X,  'PC,  II,  •)  =',  D17.5,  D40.5)  ) 

ABORT  * .TRUE. 

RETURN 

330  NYCOL  = 1 
NXSBT  = 0 

IF  (.NOT.IEQ3)  EA  = DEXP(P(1)) 

RETURN 

END 
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SUBROUTINE  STAT(V) 
IMPLICIT  REAL*8 (A-H,C-Z) 
DIMENSION 
COMMON  /Cl/ 

COMMON  /C2/ 

COMMON  /C3/ 

COMMON  /C4/ 

COMMON  /C5/ 

COMMON  /C6/ 


COMMON  /C7/ 
COMMON  /CM/ 
LOGIC  AL*  4 


V (NMAX.KCLUMN) 

N,  NMAX,  KOLUrN,  NIV,  N1VP1 
Ml  (1  1)  , I E V 

M2 (2 ) , IA,  M3(2),  NYCOL,  NYC,  NYDIV 
S (8)  , SYX  (8,8) 

AN,  A 

VMEAN(P),  Z 1 (80)  , SDEV(fl),  RM  (R,  8) 
Z2  (2)  , ELEVM,  RK  (4)  , SEYSQ 
DEI,  LEI,  DE2  , CD,  SEY,  CV,  YOEVSQ, 
IA 


FV  ALUE, 


SUBROUTINE  FOE  CALCULATING  STATISTICS 
C 

RDEVM  = O.DO 
YZERO  = O.DO 
YDEVSQ  = O.DO 
DO  2C  I = 1,  N 

V (1,NYPEV)  = V(I,NYCCL)  * V(I,NYC) 

C SUM  OF  SQUARES  OF  Y FESIDUALS 

YDEVSQ  = YUEVSQ  ♦ V(I,NYCEV)  * V(I,NYLEV) 

IF  (V  (I,NYCCL)  .NE.  O.DO)  GO  TO  10 
YZERO  = YZERC  ♦ 1.L0 
GO  TO  20 

10  RDEVM  = RDEVM  ♦ DAES  (V  (I  , NYDEV)  / V (I, NYCOL)) 

20  CONTINUE 

RDEVM  = RDEVM  / (AN  - YZERO) 

SEY SQ  = C. DO 

IF  (DEI  .GT.  O.DO)  SEYSQ  = YDEVSQ  / CF1 
C STANDARD  ERROR  OE  THE  ESTIMATE  CF  Y 

SEY  = DSQRT (SEYSQ) 

C COEFFICIENT  CF  VARIATION  (DECIMAL) 

CV  = SEY  / VMEAN  (1) 

IF  (IA)  SST  = SYX<1,1)  - 2.D0*A*S(  1)  ♦ A * A + A N 

RATIO  = C. DO 

IF  (SST  .GT.  O.DO)  RATIO  = YDEVSQ  / SST 
C COEFFICIENT  CF  DETERMINATION 

CD  = 1 . DO  - RATIO 

IFV  = 0 

DENOM  = DF2  * RATIC 
IF  (DENOM  .GT.  O.DC)  GO  TC  30 
IFV  = 2 
GO  TO  40 
C F VALUE 

30  FV  ALUE  = DF 1 * CT  / EENOH 

IF  (FVALUE  .GT.  1.CD  + 08)  IFV  = 1 
40  DO  50  J = 1,  NIV 
J1  = J ♦ 1 
DO  50  K = J1,  NIVP1 
RM  (J,K)  = 0.  CO 

DENOM  = SDEV  (J)  ♦ SDEV(K)  * OFT 

IF  (DENOM  .NE.  O.DO)  RM(J,K)  = (SY  X ( J , K)  -S  ( J ) *S  ( K ) / A N) /D  E NOM 
50  R M ( K , J ) = BM(J,K) 

RETURN 

ENU 


STATOO  10 
STAT0020 
ST  ATOO  30 
STAT0040 
STAT005C 
STAT006G 
STAT0070 
STAT0080 
ST  AT  009  C 
STAT01 00 
STATO  1 10 
ST  AT  0 1 20 
STATC1 3C 
STAT0140 
STATO 150 
STAT0160 
ST  AT  C 1 7 C 
ST AT  0 1 8 C 
STATO 190 
STAT0200 
STATC21 C 
STAT0220 
ST  A T 02 30 
ST  AT  02  40 
STAT0250 
ST  A T 0260 
ST  AT  02  7C 
STAT0280 
STAT0290 
STAT0300 
ST  AT  C 3 1 C 
STATO  320 
STAT0330 
STAT034C 
STAT0350 
STAT0360 
STAT0370 
STAT038C 
STAT0390 
STAT0400 
STAT04  10 
STAT0420 
STAT043C 
STAT0440 
STAT0450 
STAT0460 
STAT  04  70 
STAT048C 
STAT0490 
STAT0500 
STAT0510 
ST  AT  0520 
STAT0530 
STAT0540 
ST  ATC  550 
STAT0560 


' 
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SUBROUTI NE  TVAL(V) 

IMPLICIT  REAL*8 (A-H,C-Z) 

DIMENSION  V(NMAX.KCLUMN) 

COMMON  /Cl/  N,  NMAX,  KGIUrN,  NIV,  Ml  (2),  NF 
COMMON  /C3/  M2,  LINEAR,  IA,  M3(4),  NVCEV 
COMMON  /C 4/  V V (7 2),  COV  (8,8) 

COMMON  /C5/  AN,  P(8).  SE(8),  TR(8),  SICLEV(e),  EET  A (8)  , VMSQ(P) 
COMMON  /Ch/  V MEA  N ( 6)  , Z1(S),  H(6,R),  72(8),  SDFV<8) 

COMMON  /C7/  73(4),  Ck,  W W ( 2 ) , SEYSQ 
COMMON  /C8/  74,  EF 1 , Z5  (2) , SEY,  Z6,  YCEVSQ 
LOGICAL*  4 EINEAE 

DATA  PI/3.14159  26535  85793/ 

SUBROUTINE  FOR  CALCULATING  T-EEIAFED  STATISTICS 

NSTANT  = 1 * 1 A 
IF  (.NOT.  LINEAR)  GO  TO  40 
DO  10  J = 2,  NP 
BETA(J)  = 0.C0 

IF  1SEEV(1)  .NE.  O.DC)  BETA  ( J)  = P (J)  * SDEV(J)  / SDEV(1) 

10  CONTINUE 

CALCULATE  V A R I ANCE-CC V AE I ANC E MATFIX  FCP  LINEAR  CASE. 

DC  12  I * 2,  NF 
COV  (1,1)  = 0 . DO 

IE  (VMSQ(I)  .GT.  O.DC)  COV(1,I)  = -SEY  SO  * VMEAN(I)  / VMSC(I) 
COV  (1,1)  = CCV  (1  ,1) 

DC  12  J = 2,  NP 

12  COV(I,J)  = SEYSQ  * H(I-I.J-I) 

IF  (IA  .EQ.  1)  GC  TO  6 C 

CALCULATE  THE  STANDARD  ERROR  OF  P(1)  IF  P(1)  IS  NOT  SPECIFIED. 

SUM  = C . DO 

DO  30  I = 1,  NIV 

SUM1  =0.00 

DO  2C  J = 1 , NIV 

20  SUM1  = SUM1  ♦ VMEAN ( J ♦ 1 ) * H(J,I) 

SUM  1 = SUM1  * VMEAN  (1+1) 

SUM  * SUM  ♦ SUM1 
30  CONTINUE 

COV(1,1)  = SEYSg  * (1.D0  / AN  ♦ SUM) 

GO  TO  60 

CALCULATE  V A R I AN CE-CC VAR  I A NC E MATRIX  FCP  NONLINEAR  CASE. 

40  DO  50  I = NST ART , NP 
DO  50  J * NSTART , NP 
50  COV(I,J)  = SEYSQ  * H(I,J) 

60  IDF  1 = DPI 

CALCULATE  STANDARD  ERRORS,  T-RATIOS,  AND  SIGNIFICANCE  LEVELS. 

DO  170  L = NSTABT , NF 
SE(L)  = 0 . DO 

IF  (COV  (L,  L)  .GE.  0.  DO)  SE(L)  = DSQRT  (COV  (L,  L)  ) 

TR  (L)  = 1.D5C 

IF  (SE(L)  .GT.  O.DC)  IR(L)  = P(I)  / SE(L) 

AA  = CAT  AN2  (T  B (L)  , DSQRT  (DF1) ) 

IF  (IDF1  - 2)  70,  EC,  90 
70  SIGLEV(L)  = 2. DO  * AA  / PI 
GO  TO  160 

30  SIGLEV  (L)  = CSIN  (AA) 

GO  TO  160 

90  IF  ( I DF 1 .GT.  3)  GO  TO  100 

SIGLEV  (L)  = 2. CO  * (AA  ♦ CSIN(AA)  * CCOS(AA))  / PI 
GO  TO  16C 


TV  A LOO  1 0 
TVALC02C 
TVAL0030 
TV  ALOO  40 
TV  AL0050 
TV  AL006  0 
TVAL0070 
TV  AL0080 
TV  AL009C 
TV  ALO 1 00 
TVALO 110 
TV  ALO 129 
TV  AL  0 1 30 
TVAL01 40 
TVAL0150 
TVAL0160 
TVAL0170 
TVALO 180 
TVAL0190 
TVALC20C 
TVAL0210 
TV  AL  0 2 20 
TVAL0230 
TVALC240 
TV  AL0250 
TVAL0260 
TVALC27C 
TVAL0280 
TV  AL0290 
TV  AL0300 
TVAL0310 
TVAL0320 
TVAI0330 
TV  ALO  34 0 
TVAL0350 
TVAL0360 
TV  AL037C 
TVAL0380 
TVAL0390 
TV  AL04  00 
TVAL041C 
TVAL0420 
TVAL043C 
TVAL0440 
TVAL0450 
TV  AL0460 
TVAL0470 
TVAL048C 
TV  AL0490 
TVAL0500 
TVAL051C 
TVAL0520 
TV  AL0530 
TV  AL0540 
TV  AL0550 
TV  AL0560 
TVAL0570 
TVAL0580 
TV  AL0590 
TVAL0600 
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100  Z = DCOS  (AA) 

ZM  = Z * Z 

IF  (IDF  1/2*2  .EQ.  IDF1)  GC  TC  130 

II  = (ID F 1 - 3)  / 2 

SUN  = Z 

EROD  = SUN 

DO  110  I = 1,  II 

AJ  = 2 * I - 2 

PROD  = PROD  * (AJ  ♦ 2. DC)  * ZN  / (AJ  ♦ 3. DO) 
IF  (DABS  (PROD)  .LT.  1.D-10)  GO  TO  120 
SUN  = SUN  ♦ EROD 
110  CONTINUE 

120  SIGLEV(L)  = 2. DO  * (AA  ♦ DSIN(AA)  * SUN)  / PI 
GO  TO  160 

130  II  = IDF  1/2  - 1 
SUN  = 1.C0 
FROD  = SUM 
EQ  1 4C  1 = 1,  II 
AJ  = 2 * 1 - 2 

PROD  = PROD  ♦ (AJ  ♦ 1.EC)  * ZN  / (AJ  ♦ 2. DC) 
IF  (DABS  (EROD)  . LT . 1.D-10)  GO  TO  150 

SUN  = SUM  ♦ EROD 
140  CONTINUE 

150  SIGEEV(L)  = ESIN  (AA)  * SUM. 

160  SIGLEV(E)  = DABS ( 1 . D C - DABS  (SIGIEV  |L) ) ) 

170  CONTINUE 

C CALCULATE  THE  DU R B I N - W AT SCN  STATISTIC. 

DW  = 0.D0 

IF  (YDEVSQ  .EQ.  0.C0)  EETORN 
DO  180  I = 2,  N 

DIFF  = V(I.NYLEV)  - V(I-I.NYDEV) 

180  DW  = EH  ♦ DTFF  * DIFF 
DW  = EW  / YDEVSO 
RETURN 
END 


TV AL 06  1C 
TVAE0620 
TVAL063C 
TV  AE0640 
TV  AL065C 
TVAL0660 
TVAL0670 
TVAL0680 
TV  AL069C 
TV  AL0700 
TVAL0710 
TV  AL072C 
TVAI0730 
TV  AL07  40 
TVAL0750 
IVAL0760 
TV  AL077  0 
TVAL078C 
TVAL0790 
TVALC800 
TVAL08 10 
TVAL0820 
TVALC8 JC 
TVAL0840 
TV  AL0850 
TVAL0860 
TVAL087C 
TVAE0880 
TV  AL0890 
TV  AL0900 
TVALC91C 
TVALC920 
TVAL0930 
TV  AL094C 
TV  AL0950 
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/C2/ 
/C  3/ 
/C4/ 
/C5/ 
/Cb/ 
/C7/ 
/Cfl/ 


COMMON 
COBBON 
COBBON 
COBBCN 
COBBON 
COBBON 
COBBCN 
LOGIC  AL*  4 
LOGICAL*  4 
DIMENSION 
DIMENSION 
EQUIVALENCE 
DATA 
DATA 
DATA 
DATA 
DATA 


M2,  NP,  LABEL 


SUBROUTINE  OUT  1 
IMPLICIT  REAL»8 (A-H,C-Z) 

COBBON  /Cl/  N , N 1 (2) , NIV,  NIVP1, 

IEQ.  B3  ( 10)  , I E V 

ABORT,  LINEAR,  IA,  IGUESS,  NXSET,  NYCOL 
Z 1 (7  2)  , COV  (8  ,8)  , Z 2 ( 3 0)  , HEAD  (9)  , L NH  (8 ) 

Z3,  STATS  (8,5) 

V MEA  N ( 8)  , Z4  (72)  , E1(8),  SCEV(8),  PM  (8,8) 

Z5 (2) , RCEVM,  EA,  DR,  XV,  YV 

CET,  LEI,  Cf 2 , Ct,  SEY,  CV , YDEVSC,  RVALUE,  SST 
ABORT,  LINEAR,  IA,  IGUESS,  NIV1,  NIV234,  NIV567 
SMALL,  I E Q 35 

YXLABL  (8)  , BIABEL  (8)  , F(fi),  IBCD(8),  LYX  (8)  , IDG  (4) 
F V (2) 

(P(1),  STATS(1,1)),  (YXIAEL(I)  , HEAD  (2)  ) 
LYX/'Y','X1','X2','X3’,'X4','X5','X6','X7'/ 
IN/'LN'/,  CON.STANT/'  (CCN' , ' STANT) •/.  LELANK/'  •/ 

I BCD/ • A ' , ' E • , 'C • , ' I • , ' E ' , • E ' , ' G ' , *H '/,  DLANK/ ' •/ 
SPE, CIFIED/'  SPE1,  'CIFIED'/,  LA/'A'/,  DVP/O.IDO/ 
LNLPSV/'  (LN)  '/,  FV/'  > 1 C * * 8 1 , • I N F IN  ITE  • / 


SUBROUTINE  FCR  PRINTING  SUMMARY  TAELE 


OUT10010 
OUT  1C02G 
OUT  10030 
OUT  10040 
OUT10050 
OUT  1006C 
OUT  10070 
OUT10080 
OUT  1C09C 
OUT10100 
OUT  10  1 10 
OUT  1 0120 
OUT  1 0 1 3 C 
OUT  10140 
OUT101 50 
OU  T 10 1 60 
OUT  10170 
OUT101 80 
OU71C19C 
OUT10200 
OUT  10210 
OUT  1 0220 


c 

OUT  1023C 

HRITE  (6,  10) 

OUT10240 

10 

FORMAT  ( 1 HO , 49X,  'SUMMARY  TABLE') 

OUT  10250 

IE  (IEQ  .GE.  6)  HRITE  (6,  20) 

OUT  10260 

20 

FORMAT  (1H0,  34X,  ' NCTE  --  STATISTICS  ARE  EASED  ON 

LOGAR ITHMS' ) 

OUT10270 

1 

IEQ35  = IEQ  .EQ.  3 .CF.  IEQ  .EQ.  5 

OUT  10280 

DO  30  J = 1,  NIVE1 

OUT  1 02  90 

I 

ELABEL(J)  = ELANK 

OUT  10300 

LNH(J)  = LELANK 

OUT  103 1 0 

1- 

| 

IF  (LABEL  .NE.  0)  ELABEL  ( J)  = YXIAEI(J) 

OUT  10320 

30 

CONTINUE 

OUT  1 03  30 

i 

IF  (NYCOL  .GT.  1)  LN H ( 1 ) = LN 

OUT  10340 

IF  (NXSET  .EQ.  0)  GC  TO  50 

OUT10350 

DO  40  J = 2,  NP 

OUT  10360 

40 

LNH(J)  = LN 

OUT10370 

50 

IF  (LINEAR)  GC  TO  80 

OUT  10380 

HRITE  (6,  60) 

OUT  10390 

60 

FORMAT  ( 1 HO , 52X , 'STANDARD',  23X,  'SIGNIF'  / 1H  , 

• PARAMETER'  , 

OUT  10400 

1 

1 4 X , 'VALUE',  25X  , 'ERROR',  1CX,  ' T - R AT  10 • , 9X,  ' 

LEVEL'  ) 

OUT  104  10 

IF  (IEQ35)  HRITE  (6,  70) 

OUT  10420 

70 

FORMAT  ( 1H+ , 34X , 'INITIAL  GUESS') 

OUT  10430 

BAX  = 4 

OUT10440 

GO  TO  100 

OUT  10450 

80 

HRITE  (6,  90) 

OUT 1 0460 

90 

FORMAT  ( 1H0 , 52X , 'STANDARD',  23X,  'SIGNIF',  9X,  ' 

BETA'  / 

OUT  10470 

1 

1 H , 'PARAMETER',  14X,  'VALUE',  25X,  'ERROR',  10X 

, 'T-RATIO', 

OUT  10480 

2 

! 9X,  'LEVEL',  9X , 'CCEFF') 

OUT  1 04  90 

MAX  = 5 

OUT  10500 

100 

HRITE  (6,  110) 

OUT10510 

110 

FORMAT  ( 1 H ) 

OUT  10520 

LL  = LBLANK 

OUT10530 

IF  (IEQ  .EQ.  6)  LL  = LN 

OUT  10540 

f 

IF  (.NOT.  IA)  GO  TC  140 

OUT10550 

120 

HRITE  (6,  1 3 C ) LL,  LA,  SPE, CIFIED,  P(1) 

OUT10560 

130 

FORBAT  (1H  , A3,  A2,  A5,  A6,  E14.5,  15X,  3F15.5) 

OUT  10570 

GO  TO  190 

OUT10580 

140 

IF  (IEQ  .NE.  6)  GO  TC  160 

OUT  10590 

DEXPA  = DEX  P ( F ( 1 ) ) 

OUT10600 
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SMALL  UKXPA  . LT  . t)  V t .Oil.  DEXIA  .01.  1 . U «> 

IF  ( SPALL  ) WRITE  (6,  ISO)  PLANK,  LA,  btANK,  f LANK,  LEXPA 
XL'  {.NOT.  SMALL)  WRITE  (b,  1 3U)  PLANK,  LA,  PLANK,  PLANK,  CEXPA 
150  FORMAT  ( 1H  , A3,  At,  A5,  Ab,  E 1 8 . 5 , 1 1 X , Ifl'..  5) 

IF  (XA)  GO  TC  190 

IbO  SMALL  = LA 83  (P  ( 1 ) ) .IT.  DVP  .OF.  PALS(F(1))  .GT.  1 . Lb 

XF  (.NOT.  SMALL)  WKI1F  (b,  130)  LL, LA, CON, START,  (S?A1G(1,R),K-1 

IF  { SPALL  ) WRITE  (6,  160)  LL  , L A , LON , ST  A NT , (STATS  { 1 , K)  , K*  1 

IF  ( . NOT  . I EG  36)  GO  TC  190 

SMALL  = CABS  (PI  (1)  ) .LT.  LVP  .CF.  CABS(P1(1))  .GT.  I.Lb 
IF  (.NOT. SMALL)  WHITE  (fc,  17C)  f 1 ( 1 ) 

IF  < SMALL  ) WHITE  <6,  1 8 C ) ri(1) 

170  FORMAT  ( 1 H ♦ , TJ2,  115.5) 

180  FORMAT  { 1 H ♦ , T32,  LIS. 5) 

190  DO  220  J = 2,  NP 

SMALL  = DABS  < P ( J ) ) .LT.  DVP  .OF.  LABS  (P (J) ) .GT.  I.Db 
IF  ( S MALL  ) WRITE  (6  , 200)  LDCD(J),  YXIABL(J), 

1 (STATS  (J,K)  , K = 1,  MAX) 

IF  (.NOT.  SMALL)  WRITE  <b,  210)  LtsCC(J),  YXLAEL(J), 

1 (STATS  (J,K)  , K = 1,  PAX) 

200  FORMAT  (4X,  A4,  Ab.  T19.5,  11X,  4F15.5) 

210  FORMAT  (4X,  AM,  A3,  X15.6,  15X,  4F15.5) 

If  { . NOT  . I EG  3 5)  GC  TC  220 

SMALL  = CABS  (XI  (J)  ) .LI.  LVP  .CF.  CATS  (Pi  (J)  ) .GT.  1.06 
IF  (.NOT. SMALL)  WRITE  (6,  170)  F1(J) 

IF  ( SMALL  ) WRITE  (6,  180)  I1(J) 

220  CONTINUE 

WRITE  (6,  23  C ) 

230  FORMAT  (1H0) 

N I V 1 = NIV  . EC.  1 

NIV234  = NIV  .LF.  4 .AND.  .NC.T.  NIV1 

NIV5b7  = NIV  .GE.  5 

IF  ( N I V 1 ) WRITE  (6  , 240) 

240  FORMAT  (1H0,  37X,  'STANLARD'  /IP  , 'VARIABLE*,  16X,  'MEAN',  9X , 
1 'DEVIATION') 

IF  ( N I V 2 34)  WRITE  (6,  250)  (LNH(J),  J * 1,  NIVFl) 

IF  ( N I V 2 3 4 ) WRITE  (fc,  260)  (YXLAEL(J),  J = 1,  N I VP  1 ) 

250  FORMAT  ( 1 H , 62X,  'CORRELATION  MATRIX*  /1H  , 37X,  'STANDARD*,  1 

1 b(A2,1JX)) 

260  FORMAT  (1H  , 'VARIABLE',  15X,  'MEAN',  9X,  'DEVIATION*, 

1 5(7X,A6)) 

IF  (NIV5b7)  WRITE  (b , 270)  (LNH(J),  J = 1,  N I VP  1 ) 

IF  (NIV567)  WRITE  (6,  280)  (YXIAEL(J),  J = 1,  NIVP1) 

270  FORMAT  (1H  , 73X,  'CORRELATION  MATRIX'  /1H  , 37X , 'STANDARD',  . 
1 S ( 8 X , A 2 ) ) 

280  FORMAT  (1H  . 'VAPIABIF',  15X,  'MEAN',  9X,  'DEVIATION', 

1 5 X , 8(2X,A8) ) 

WRITE  (b,  11C) 

DO  320  J = 1 , NI VP  1 

WRITE  (6,  290)  LNH(J),  LYX(J),  ELABEL(J),  VMEAN(J),  SDEV  (J) 

IF  (NIV2  34)  WRITE  (6  , 300)  (RM(J,K),  K = 1,  NI  V P 1) 

IF  (NIV567)  WRITE  (6,  310)  (RM(J,K),  K = 1,  N I VP  1 ) 

290  FORMAT  (1H  , A3,  A4,  A8,  2F15.5) 

300  FORMAT  ( 1 H ♦ , 45X,  5F15.5) 

310  FORMAT  ( 1 H ♦ , 50X , 8F10.5) 

320  CONTINUE 

LNLR  = LELANK 

IF  (LNH  ( 1)  .EQ.  LN)  LNLR  = LNLRSV 
WRITE  (6,  230) 

WRITE  (6,  3 J C ) CD,  LNLR,  RDEVM 

330  FORMAT  (IX,  'COEFFICIENT  CF  DETERMINATION  (UNADJ),  R SQ',  F13.1 


n n n 


\ ' 


M3. 5 / 


1 SX.  'MEAN  Of  ABSOLUTE  RELATIVE  EFVIATICNS  ',  A4,  FU(.i| 
WRITE  (6,  34C)  SEY,  C V,  LKLR , YtEVSC,  LNIR,  SST 
340  FORMAT  (IX,  'STANDAhl  ERROR  Cf  1ST  IRATE',  F29.5, 

1 5X,  'COEFF  VARIATION  {STL  FRR  ESI  / KEAN  Y CBS)', 

2 IX,  'SUM  OF  SQUARES  CF  RESICUAIS  *,  A4,  F23.5, 

3 5 X , 'SUM  OF  SQUAHES  TOTAL  ■,  A4,  F3C.5) 

If  (IFV  .EQ.  0)  WRITE  (6,  350)  fVALUf 

IF  (IfV  .GT.  0)  WRIT  f (6,  360)  fV(IFV) 

350  FORMAT  (IX,  *E  VALUE',  34X,  F14.5) 

36C  FORMAT  (IX,  *f  VALLE',  40X,  AH) 

WRITE  (6,  370)  DW 

370  FORMAT  (1H  + , 60X,  ' D L RB I N -W A T SC N STATISTIC',  f32.5) 

IF  (TEC  .EP.  5)  WRITE  (6,  380)  EA 
380  FORMAT  (IX,  'Y  INTERCEPT* , F44.5) 

IF  (IE0  .EQ.  2)  WRITE  (6,  390)  XV,  YV 
390  FORMAT  (IX,  'X  CCORCINATE  OF  VERTEX',  F33.5, 

1 5 X , 'Y  CCORCINATE  CE  VERTEX',  F33.5) 

I DG ( 1 ) * DF1 
IDG  ( 2 ) = DF 2 
IDG  ( 3)  = DFT 
IDG  (4 ) = N 
WRITE  (6,  40C)  IDG 

400  FORMAT  (IX,  ' DEGREES  Of  FREEDOM  FOR  ERROR',  127, 

1 5 X , 'DEGREES  CF  FREEDOM  CUE  TO  REGRESSION',  119  / 

2 IX,  'TOTAL  DEGREES  CF  FREEDOM',  131, 

3 5 X , 'NUMBER  OF  DATA  POINTS',  134) 

IF  (NP  .EQ.  1)  RETURN 


PRINT  VARIANCE-COVARIANCE  MATRIX 


410 

420 


430 


440 

450 


460 

470 


480 

490 


NSTART  = 1 

IF  (IA)  NSTART  = 2 

IF  ( N1V567)  WRITE  (6,  410) 

FORMAT  (1H0/  T42,  'VARIANCE-COVARIANCE  MATRIX*  / ) 

IE  (.NOT.  N I V 567 ) WRITE  (b,  420) 

FORMAT  <180/  T12,  'VARIANCE-COVARIANCE  MATRIX'  / ) 

IF  (IEQ  .EQ.  6)  GO  TO  460 

WRITE  (6,  430)  (LECD(J),  J = NSTART,  NP) 

FORMAT  (T  3,  6 ( 1 3 X , A1)  ) 

DO  450  1 = NSTART,  NP 

WRITE  (6,  44  C)  LECE(I),  (COV(I,J),  J = NSTART,  NP) 

FORMAT  (T6,  Al,  2X,  6(2X,  D12.5)  ) 

CONTINUE 

RETURN 

IF  (IA)  WRITE  (6,  43C)  (LBCD(J),  J = 2,  NP) 

IF  (.NOT.  IA)  WRITE  (6,  470)  (LECD(J),  J = 2,  NP) 

FORMAT  (T 1 3 , 'LN  A',  7{13X,  Al)  ) 

DO  490  I = NSTART,  NP 

IF  (I  .EQ.  1)  WRITE  (6,  480)  (CCV(I,J),  J = 1,  NP) 

FORMAT  (13,  'LN  A',  ?X,  8 (2X , C12.5)  ) 

IF  (I  .GT.  1)  WRITE  (6,  440)  LECD(I),  (CCV(I,J),  J = NSTART,  NP) 

CONTINUE 

RETURN 

END 


OUT  1 1 220 
OUT1 1230 
OUT  1 1240 
OUT1 1250 
CUT  11260 
OUT  11270 
OUT  1 1 2 80 
OUT  11290 
OUT  1 1 300 
OUT  1 1 3 1 C 
OUT  11320 
OUT1 1330 
OUT  1 1 3 4 C 
OUT1 1350 
OUT  11360 
OUT  1 1370 
CUT11380 
OUT  11390 
OUT1 1400 
OUT1  1410 
OUT  1 1420 
OUT  11430 
OUT  1 1440 
OUT  114  50 
OUT  11460 
OUT  1 14  7C 
OUT1 1480 
OUT  1 1490 
OUT1 1500 
OUT  1 151C 
OUT  1 1 520 
OUT  1 1530 
OUT1 1540 
OUT  11550 
OUT  1 1 56  0 
OUT11570 
OUT1 1580 
OUT  1 1590 
OUT1 1600 
OUT  1 1 6 1 0 
OUT  1 1 6 20 
OUT  1 1630 
OUT1 1640 
OUT  11650 
OUT1 1660 
OUT  1 16  70 
OUT  1 1680 
OUT1 1690 
OUT  1 1700 
OUT11710 
OOTI 1720 
OUT  1 1730 
OUT11740 
OUT  1 1750 


COMMON  /Cl/  N,  N MAX,  KOLUMN,  NIV,  NIVPI,  NVAR,  NP,  LABEL,  ID  OUT2004C 
COMMON  /C2/  I EC,  NF1,  NP2 , NP3,  LZZYES,  Ml,  INCUT,  H2(16),  N PAGE  OUT20050 
COMMON  /C3/  A ECN I , LINEAR,  IA,  M3,  NXSET,  NYCCL,  NYC,  NYDEV  OUT20060 
COMMON  /C4/  Z 1 (1 66)  , BEAD  (9)  , 1 KH ( 8)  OUT20070 


COMMON  /C 5/  AN 

COMMON  /C7/  Z 2 (2 ) , REEVE! 

COMMON  /CX/  E XV ( 8)  , NEXV,  NV1 

DIMENSION  ELTER(2),  HEALS  (8)  , TR  A N ( 8)  , VHD(6),  LETR  ( 35)  , LNX(8) 
LOGICAL*  4 ABORT,  NE1,  NF2,  LNCUT,  IICTS,  ID,  SKIP 
ECUIVALENCE  (YIN,  I N E ( 1 ) ) , (YP,  HEAC(2)) 

LATA  IELANK/'  •/,  ELANK/'  •/ 

DATA  CBSERV/' OBSERVED'/,  COM P/ • CC MPUT EE •/ , BBC/'  E'/, 

1 RESIE/ 'RESIDUAL'/,  R E L/ • R F L AT  I V E •/ , BD/'  D'/, 

2 EVIATN/' EVIAT ION •/ 


OUT  20080 
OUT20090 
OUT2010C 
OUT201  10 
OUT20120 
OUT  20 1 30 
OUT  20 1 4 0 
OUT20150 
CUT2016C 
OUT  20 1 7C 


C 

C 

c 

c 


DATA 

LETR/ ' A ' , ' E ' , 

•C'  , 'C  ,'E*  , 

•E','G','H','I','J' 

, *K' , 'L* ,'M' 

1 

•C'.'P*. 

•C', ' R ' , • S ' , 

,'Y','Z','1' 

2 

•3',  '4', 

•5', '6'  ,'7', 

•f,  '9'/ 

DATA 

TRAN/' 

Y : ' , 'XI  : 

',  ' X2:  ',  • x 3 : ', 

' X4:  ',  ' X5: 

1 

•XE:  «,  'X7: 

V 

DATA 

ALTER/ 

• MODIFIED' , 

•NOT  USED'/ 

DATA 

VHC/  • 

<X2)  ' , • 

(X3)  ' , • ( X 4 ) ' , • 

(X5) 

1 

• 

(Xb)  • , • 

<X7)  '/ 

SUBROUTINE  POP  PRINTING  INPUT  DATA,  COMPUTED  Y VALUES, 

Y RESIDUALS,  AND  PERCENT  Y DEVIATIONS 

NIVPS  = NIVP1 
IF  (NEXV  .EQ.  0)  GC  TO  5 
NIV  = NIV  ♦ NEXV 
NIVPP  = NIVP  1 ♦ 1 
NIVPI  = NIV  ♦ 1 
5 ITEMP  * NP  ♦ I A 

IF  (NP  .EC.  1)  ITEMP  = 0 

LINES  * NP  ♦ NIVF1  ♦ N ♦ N/5  ♦ IINEAF  ♦ IA  ♦ IA  ♦ ITEMP 
IF  ( N PAG  E .EC.  2)  LINES  = LINES  ♦ 2 


IF  (IEQ 

• GT.  6)  LINES  = 

LINES  » 

2 

IF  ( I EC 

.EC.  6)  LINES  = 

LINES  ♦ 

3 

IF  ( I EC 

• EQ.  2)  LINES  = 

LINES  * 

1 

IF  (LABEL  .EC.  2)  LINES 

= LINES 

* 

PLOTS  = 

NP 1 .OR.  NP2  .CB 

. NP  3 . 

GT 

SKIP  = 

NYCOL  .EO.  1 .OR. 

LNOUT 

DEVNAX 

= -1.C10 

DEVHIN 

= 1.C10 

RDEVH  = 

0.  DO 

NREPS  - 

(N-1)/35  ♦ 1 

NYPLOT  = NYCCL 
NXPLOT  = NXSET 
NP3PL  = NP3 
NP3E0  = 0 

IF  (IEfi  .EC.  0)  GO  TC  260 
DO  10  I = 1,  NIVPI 
HEADS  (I)  = BLANK 
10  LNX  (I)  = LNH  (I) 

HEADS  (1 ) = CESERV 

IF  ( NXSET*  NYCCL  .EC.  1 -OR.  INCCT)  GC  TC  30 
DO  20  I = 1,  NIVPI 
20  LNX  (I)  = I BLANK 


' N • , OUT  20 1 80 
•2 ' , OUT20190 
OUT2020C 
’,  OOT20210 

OUT20220 
00T2C23C 
OUT20240 
OUT20250 
OUT  20260 
OUT20270 
OUT20280 
OUT20290 
OUT203C0 
OUT203  10 
OUT20320 
OUT2033C 
OUT20340 
OU  T 20  350 
OUT20360 
OUT20370 
OUT20380 
OUT20390 
OUT20400 
OUT2041C 
OUT20420 
OUT204  30 
OUT20440 
OUT20450 
OUT20460 
OUT20470 
OUT20480 
OUT204  90 
OUT20500 
OUT205 10 
ODT20520 
O0T20530 
OUT20540 
OUT20550 
O0T2056C 
ODT20570 
OUT20580 
O0T20590 
OUT20600 
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I 

I 


NYCCL  = 1 
N X5ET  = C 

30  IF  (NVT  .Eg.  C)  GC  1C  Jb 
EO  32  I = 1,  NIVF1 

IF  (EXV(I)  .NE.  0.  DO ) HEADS  ( I)  * AITEF(I) 

32  CONTINUE 

IF  (NEXV  . EQ  . 0)  GC  1C  36 
DO  34  I = NIVFP,  NIVF1 
EXV(I)  = O.DC 
HEAD(I*1>  = VKD (1-2) 

HEADS  (I)  = ALTER  ( 2 ) 

LNX(I)  = IELANK 
34  CONTINUE 

36  IF  (LINES  .LE.  21)  GC  TO  40 
C PRINT  TITLE  CN  NEW  PAGE 

CALL  TITLE2 
40  WHITE  (6,  50) 

50  FORMAT  (/1H0,  49X,  ' TABLE  OF  RESIDUALS') 

LINES  = 0 
NTIMES  = 1 

60  IF  (NIV  .GT.  4)  GC  TC  90 

WRITE  (6,  70)  (HEADS(J),  J = 1,  N1VP1),  COUP,  RESID,  R EL 
WRITE  (6,  80)  HEAD  ( 1 ) , (LNX(J),  HEAC(J*1),  J-1.NIVP1), 

1 LNX  ( 1)  , YP,  LNX(1),  YF,  EED,  EVIATN 
70  FORMAT  (1H0,  10X,  B(7X,  AH)  ) 

80  FORMAT  (1H  , 2X,  AF,  e<4X,A3,A8)  ) 

GO  TO  120 

90  WRITE  (6,  100)  (HEADS  ( J)  , J * 1,  NIVP1),  CCMP,  FESID,  REL 
WRITE  (6,  110)  HEAE(I),  (LNX(J),  HEAt(J*1),  J=1,NIVP1), 

1 LNX(1),  YP,  INX(1),  YF,  EP,  EVIATN 
100  FORMAT  ( 1 H 0 , 10X  , 11  (3X,  AH)  ) 

110  FORMAT  ( 1 H , 2X , AH,  11(1X,A2,AE)  ) 

120  WHITE  (6,  20C) 

IF  (NTIMES  .EL.  2)  GC  TO  140 

YZERO  = C. DO 

DO  210  I = 1 , N 

IF  (LINES  .LT.  40)  GC  TO  140 

CALL  TITLE2 

LINES  = C 

WRITE  (6,  130) 

130  FORMAT  (/1H0,  44X,  'TABLE  OF  RESIDUALS  (CONTINUED)*) 
NTIMES  = 2 
GO  TO  60 
140  LET  = IBLANK 

IF  (PLOTS)  LET  = LET  H (I- 1) /NR  E PS  ♦ 1) 

VID  » BLANK 

IE  (ID)  VID  = V ( I , NV  AN) 

YC  * V ( 1 , N YC  ) 

YDEV  = V (I.NYCEV) 
ir  (SKIP)  GO  TO  150 
YC  = DEI F (V (I, NYC)) 

YDEV  = V(I,1)  - YC 

150  IF  ( V (I , NYCO L)  . NE.  C.DO)  GO  TC  160 
YZERO  - YZERC  ♦ 1 . CO 
R DEV  = O.DO 
GO  TO  170 

160  RDEV  = YDEV  / V(I,NYCOL) 

170  RDEVH  = RDEVM  ♦ DABS  (RDEV) 

IF  (DEVHAX  .LT.  RDEV)  DEV  MAX  = RDEV 

IF  (DEVB1N  .GT.  RDEV)  DEV  MIN  * RDEV 

IF  (NIV  .LE.  4)  WRITE  (6,  180)  l FT,  VID,  V(I,NYCOL), 


OUT206  1C 
OUT206  20 
OUT2063C 
ODT20b40 
OUT206  50 
OUT20660 
OUT2067C 
OOT2068C 
OUT  206  90 
ODT20700 
OUT  207 1 C 
OUT  207  20 
OUT207  30 
OUT  2074C 
OUT2075C 
OUT  2076C 
OUT20770 
OUT2C78C 
OU  T 20790 
OUT20800 
OUT20H  1C 
OUT20820 
OU 1208  30 
OUT20840 
OUT20850 
OUT2C86C 
OUT20870 
OUT2088C 
OUT20890 
OUT20900 
O0T2C910 
OUT20920 
OUT209  30 
OUT20940 
OUT2C950 
OUT20960 
OUT2097O 
OUT20980 
OUT20990 
OUT21000 
OUT21010 
OUT  2 1 020 
OUT21030 
OUT21040 
OUT21050 
OUT21060 
ODT21070 
OUT21080 
OUT21090 
OUT21100 
OUT21 1 10 
OUT21120 
OUT21130 
OUT2  1 140 
OUT21150 
OUT  2 1 1 6 C 
OUT21170 
OUT  2 1 1 80 
OUT21190 
OUT21200 
OUT  2 1 2 10 


| 


-102- 


1 (V  (I.J4NXSLT) , J = 2,  NIVPl),  YC,  YCEV,  RDEV 
IF  (NIV  .GI.  4)  WRITE  (6,  190)  IET,  VIC,  V(I,NYC0L), 
1 (V (I ,J4 NXSET) , J = 2,  NIVPl),  YC,  YCEV,  RDEV 
) FORMAT  ( 1H  , A 2,  At,  8F15.5) 

) FORMAT  ( 1 H , A2,  AP,  11F11.3) 

LISES  = TINES  ♦ 1 

IF  (1/5*5  .EC.  I)  WRITE  (f. , 200) 

} FORMAT  ( 1 H ) 

) CONTINUE 

RDEVM  = RDEV  M / (AN  - YZERO) 

WRITE  (fa,  220)  DEVMIN,  RCEVM,  CEVM.AX 
) FORMAT  ('OMINIMUM  RELATIVE  DEVIATION  = ',  F10.5, 


F10.5,  ' 

I 'MEAN  ABSOLUTE  RELATIVE  DEVIATION  F8.5,  ', 

i 'MAXIMUM  RELATIVE  DEVIATION  =',  F10.5) 

IF  (NVT  .EC.  0)  GO  TO  226 

IF  (N  IV  PS  .LE.  fa)  WRITE  (t,  222)  (T  S AN  (I ) , EX  V (X)  , 1 = 1 , N I VPS) 

IF  (NIVPS  .GE.  7)  WRITE  (fa,  224)  ( T R A N 1 1 ) , EX  V ( I)  , 1=  1 , N I V F S) 

FORMAT  (/•  TRANSFORMATION  FACTORS  --  6(A4,  F10.5,  2X)  ) 

FORMAT  (/•  TRANSFORMATION  FACTORS  --  ',  6(AU,  E10. 5,  2X)  / T28, 
1 2 ( A 4 , 1 1C. 5,  2X)  ) 

IF  (.NOT. HOTS)  RETURN 
IF  ( . NOT . N P 1 ) GO  TO  *40 
CALL  TITLE2 

WRITE  Ife,  23C)  RESIC,  YLN  , YP,  COMP,  Y L t. , YP 

FORMAT ( 1 HO,  43X,  A8,  IX,  A3,  At,'  VERSUS  >,  AS,  IX,  A3,  A 


OUT21220 

COL) , OUT21230 

OUT2  1240 
OUT21250 
OUT2126C 
OUT2127C 
OUT  2 1 280 
OUT21290 
OUT21300 
OUT2  13  10 
OUT21320 
',  OUT2133C 

•,  •,  OUT  21340 

OUT21350 
CUT  2 1 36C 

| ,1  = 1 , NIVPS)  OUT 2 1 37C 

) ,1=1,NIVFS)  OUT 2 1 380 

10.5,  2X)  ) OUT21390 

10. 5,  2 X)  / T28 , OUT21400 

OUT2  14  10 
OUT2142P 
OUT  2 1 430 
ODT2  1440 
OUT21450 

A8 , IX,  A3,  A8/)OUT214feO 


CALL  ELOTYX  ( N , V(I.NYCEV)  , V ( 1 , N YC)  , C , 0 ) OUT21470 

IF  (.NCT.NP2)  GO  TO  25C  OUT21480 

CALL  TITLE2  OUT21490 

WRITE  (6,  2 3 C ) OBSERV,  YLN,  YP,  COME,  YLN,  YE  CUT2150C 

CALL  ELGTYX  ( N , V(I.NYPIOT),  V ( 1 , N YC ) , C , 0 ) OUT21510 

IF  (NP3PL  .EC.  0)  RETURN  OUT21520 

IF  (NP3PL  -LE.  NIV)  GO  TO  310  OUT21530 

r F (NF3PL  .LT.  8)  RETURN  OUT21540 

SPJPL  = 1 OUT21550 

IE  (NIV  .EC.  1)  NP3EC  = 1 OUT2156C 

GO  TO  310  OUT 2 1 570 

FOLLOWING  SECTION  FOR  FLOT-OMY  CPTICN  OUT21580 

IF  (NP3  .Eg.  0 .ca.  NP 3 .GE.  8)  NP3PL  = 1 OUT21590 

IF  (NP3PL  .GT.  NIV)  FETURN  OUT21fa0C 

LINES  = 0 OUT  2 16 10 

WRITE  (6,  27  C)  HEAC(1),  H EAD  ( 2)  , H E A C ( K P 3EL4  2)  OUT21620 

FORMAT  ( 1 H 0 , 17X,  'CESEHVED’/  1H  , 2X  , A8,  7X,  A8,  7X,  A8  / ) OUT2163C 

LNH(NE3PL4l)  = IELANK  OUT21640 

L N H ( 1 ) = IBLANK  OUT21650 

DO  300  I = 1,  N OUT21660 

IF  (LINES  .LT.  4 C)  GO  TO  280  OUT2167C 

CALL  TITLE2  OUT21680 

WRITE  (6,  27C)  HEAC(I),  HEAD (2)  , H I A C ( N F 3 F L 4 2)  OUT21fa90 

LINES  = C OUT  2 1 70 C 

LET  = LETR  ( (1-1) /NREES  * 1)  OUT21710 

VID  * BLANK  OUT21720 

IF  (It)  VID  = V ( I , NV  AR)  OUT21730 

WHITE  (6,  29C)  LET,  VID,  V(I,1),  V (I , NF3FL* 1)  OUT21740 

FORMAT  ( 1 H , A2 , A6,  2F15.5)  OUT21750 

LINES  = LINES  4 1 OUT21760 

IF  (1/5*5  .EC-  I)  WRITE  (6,  200)  OUT21770 

CONTINUE  O0T21780 

CALL  TITLE2  OUT21790 

WRITE  (6,  230)  0 BS  E R V , LN  H ( 1)  , Y P , 0 E S E P V , LN  H ( NF  3PL*  1 ) ,H  EA  D (NP  3P  L*2 ) OUT2  1 800 
LZERO  = LZZYES  OUT21810 

CALL  ELOTYX  ( N,  V(I.NYPLOT),  V ( 1 , N P 3PL * 1 *NX PICT) , NF3EQ,  LZERO)  OUT21820 
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subroutine  PLCTYX  ( N,  Y,  X,  NEC,  LZ  EEC  ) 

IMPLICIT  HEAl*8 (A-H,C-Z) 

E IMF  NS  ION  Y (N)  , X (N) 

COMMON  /C2/  I EC 
COMMON  /C3/  Ml,  LINEAR 
COMMON  /C5/  Z 1 , A,  B,  C 

COMMON  /C9/  CUMtlY  (1535)  , OUTPUT  j 15 , 45)  , LI,  LX(3) 

LOGICAL*  4 LINEAR,  LZEFO 

LOGICAL*  1 CCLRCU  (120,45),  LI,  LX,  NUMBER,  ISYMBL,  LETR  (35) 
EQUIVALENCE  (CCLEON(I),  OUTPUT(I)),  (LI,  IN) 

V 

K'.'L'.’M' 
Y *,  *Z  * , ’ 1* 


1 *0*, 

2 ’3', 

NREPS  = (N-1)  / 35 

YMAX  = Y (1) 

YMIN  = Y (1) 

XMAX  = X (1) 

XMIN  = X (1) 

DC  10  I = 2, 


/. 

ISYMBL/’ . 

V 

V, 

, IELANK/’ 

V, 

ELA 

NK/’ 

’/,  ZFROL/ 

• E ■ 

’.’C’.’C’, 

•E’, 

*F’, 

>C’  , 

•H’.’I’.'J' 

,p. 

, *C',  ’R* , 

•S’, 

' T ’ , 

’ U ' , 1 

'V', ’W* ,’X* 

•4  1 

1 » ' - 1 . '6 1 » 

•7’, 

’6’, 

• 9 '/ 

1 


N 


IF 

(X 

(I) 

. LT. 

XMIN) 

XMIN 

= X 

(I) 

IF 

(Y 

(I) 

.LI. 

YMIN) 

YMIN 

= Y 

(I) 

IF 

(X 

(I) 

• GT. 

XMAX) 

XMAX 

= X 

ID 

IF 

(Y 

(I) 

.GT. 

YMAX) 

YMAX 

= Y 

(I) 

CONTI 

NUE 

PLOT  FROM  ORIGIN  IF  SELECTED 
IF  (LZERC)  XMIN  = 0.  10 
IF  (LZERC)  YMIN  = C.CO 

FIND  SCALE  FACTORS  FCk  45  LINES  AND  120  SPACES 
XSCAL  = (XMAX  - XMIN)  / 120. CO 
YSCAL  = (YMAX  - YMIN)  / 45. DO 
WRITE  (6,  20)  YMAX 
20  FORMAT  (’ONAX  VERT  = • ,F 14. 5/) 

WRITE  (6,  30) 

30  FORMAT  (1H«,10X,1H  ,12(' .’)) 

FIND  RECIPROCAL  CF  SCALf'fOR  MULTIPLICATION 
RXSCAL  = 1 . DC  / XSCAI 
RYSCAL  = 1.DC  / YSCAI 

ELANK  OUT  PLCT  PAGE  WITH  EQUIVALENT  e-EYTE  BLOCKS 
DO  40  I = 1,  15 
DO  40  L = 1,  45 
40  OUTPUT (K , L)  = BLANK 
14  = IBLANK 

DETERMINE  WHERE  ZERO  LINE  IS,  IF  ANY 
NZERO  * 0 

IP  (YMIN  * YMAX  .LT.  C)  NZERO  = YMAX  * RYSCAL  * 1 
IF  (NZERO  .EC.  0)  GO  TO  6C 
DO  50  K * 1,  15 
50  OUTPUT(K, NZERO)  * ZFRCL 
60  DO  70  I * 1,  N 

DETERMINE  HORIZONTAL  AND  VERTICAL  CHARACTER  POSITION 
K = (X(I)-XMIN)  » RXSCAL  ♦ 1 
L = 45  - (Y  (I) - YMIN)  * RYSCAI 

ACCEPT  RIGHT-HAND  PLCT  BOUNDRY  AS  WITHIN  LAST  CHARACTER 
POSITION  AND  UPPER  ECUNCRY  AS  WITHIN  FIRST  LINE 
IF  (K  .Gl.  120)  K * 120 
IF  (L  .LT.  1)  L * 1 

MOVE  PLOT  CHARACTER  TO  BEGINNING  OF  INTEGER  TEST  WORD 


LI  * COLROH  ( K , L) 


PLOTOOIO 
PLCT0020 
PLCTOO  30 
PLOTOO  40 
PLCTC050 
PLOTC06C 
PLOTOO  7 0 
PLCT0080 
PLCT009C 
PLCT010C 
PLOTO  110 
PLCT012C 
, 1 N ■ , PLOT01 30 
, 1 2 ' , PLOTO 1 4C 
PLOTO 1 50 
PLOTO 1 60 
PLCTC1 70 
PLCT  0 1 8C 
PLCT0190 
PLCT020C 
PLCT  02  10 
PLOT0220 
PLOT  0 2 3C 
PLCT02  4 0 
PLCT0250 
PLOT0260 
PLCT027C 
PLOT0280 
PLOT0290 
PLCT0300 
PLCT03  1 0 
PLOTO  320 
PLOT03  3 0 
PLOT034C 
PLOT0350 
PLO70360 
PLCT037C 
PLOTC3  80 
PLOTG390 
PLOT0400 
PLOT04 10 
PLOT0420 
PLOT04  30 
PLCT0440 
PLCT0450 
PLOT0460 
PLOT0470 
PLOTC480 
PLCT0490 
PLOT0500 
PLOT05 1 0 
PLCT0520 
PLOT0530 
PLOT054  0 
PLOT055C 
PLCT0560 
PLOT0570 
PLOT0580 
PLCT0590 
PLCT0600 


‘ 

1 
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FOR  OVERLAPPING  POINTS,  FLOT  » 

IF  (14  .NE.  IELANK  .AND.  14  .NE.  ILINF)  COLROW(K.L)  = NUHEER 
FLOT  SINGLE  FICT  POINT  AS  ALPHANUMERIC 

IF  (14  .EO.  IBLANK  .CR.  14  .EO.  ILINE)  CCLRCN|K,L)  = LETR((I-1)/ 


1 NREPS  ♦ 1) 
70  CONTINUE 


80 


I') 


90  CONTINUE 

100  WRITE  (6,  TIC)  OUTPUT 
110  FORMAT  ( 1 H , 10X,  15A8 

WRITE  (6,  30) 

WRITE  (b,  1 2 C)  YHIN,  XMIN,  XMAX,  YSCAL,  XSCAl 
120  FORMAT  (*0MIN  V ERT=  • , F 1 4.  5/ • MIN  IIOHZ  = ’ , F 1 4.  5 , 8bX  , • M AX  HORZ=' 
1 F14.5/'CVERT  INCREMFNT= • ,F1 2. 5/'  HCRZ  INCFEMFNT-' , F 1 2 . 8) 
RETURN 
END 


PLOT06  10 
PLOT  C6  20 
PLCT06  3 0 
PLOT064C 
PLCT0650 
PLCT066  C 


IF 

(NEC  .NE. 

1)  GO  TO  100 

PLCT0670 

XX 

= XMIN  - 

C.5DC  * 

XSCAL 

PLOT  06  80 

DO 

9 0 K = 1 , 

120 

PLOT  06  9C 

XX 

= XX  ♦ XSCAL 

PLCT0700 

YY 

= A ♦ B * 

XX 

PLCT07  10 

IF 

(LINEAR) 

GC  TO  80 

PLOT  0720 

IF 

( I EO  .EC. 

2)  YY  = 

YY  ♦ C ♦ XI 

* 

XX 

PLCT073C 

IF 

(IEC  .EC. 

3)  YY  = 

A ♦ XX*  * E 

PLOT0740 

IF 

(IEC  -EC. 

4)  YY  = 

A ♦ E * XX**C 

PLCT0750 

IF 

(IEC  .FO. 

5)  YY  = 

DEXP  (YY) 

PLOT  u76C 

L = 

45  - (YY 

- YMIN) 

* PYSCAI 

PLCT0770 

IF 

(L  .LI.  1 

• OR  . I 

.GT.  45)  GC 

TO 

90 

PLOT  C78C 

LI 

= COLROW  ( K , L) 

PLOT  07  90 

IF 

(14  .EC. 

IfcLANK  . 

CH.  14  .EC. 

ILINt)  CCL  ROW  (K , L)  = ISYMEL 

PLCT0800 

PLCT08  10 
PLCT0820 
PLOT 08  3 C 
PLCT0840 
PLOT085C 
PLOT  086  0 
PLCT0870 
PLOT  C 8 8 C 
PLOT  0890 


1 
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ASSEMBLER  ROUTINE  ****************  MEMR0010 


README RO 

START 

0 

MEMR002C 

ENTRY 

MEflRE 

MEMH0030 

ENTRY 

DATEP 

MEMH0040 

EXTRA 

IBCCM* 

MEHP0050 

EXT3N 

F IOC  S * 

NEMROCbO 

♦ 

♦ 

KEMRQC7C 

* 

FIRST  ENTRY 

POINT.  THIS  ROUTINE  PICKS  UP  OPERANDS 

ONE 

♦ 

MEMR008C 

♦ 

AND  TWO  AND 

STORES  THEM  AT  THE  TWO  FULL  WCP13  AT  EUFLOC* 

M E MPOO  90 

* 

♦ 

MEMR01CC 

USING 

».  15 

MEMRC1 10 

HEflH  E 

B 

*♦12 

EFANCH  AROUND  NAME 

MEMP0120 

DC 

X L4  * C 70  00  CO  C 

* ROUTINE  NAME  FOR  CALI  TRACE 

MEMRO 1 30 

DC 

C LA  'CORE ’ 

ROUTINE  NAME  ECR  CALL  TRACE 

MEMR01 40 

SIM 

14.3,12  (13) 

SAVE  REGISTERS 

MEMR0150 

LN 

2,  J,C(1) 

FETCH  CIERANL  ACDEESSES 

MEBRO 160 

L 

3,0  (3) 

FETCH  CFERANL  2 (LENGTH) 

MEHR0170 

STM 

2, J,EUFACR 

STORE  EUFFBF  ICC  AND  LENGTH 

MEMP018C 

LA 

1 , CCRE2 

R1=A (SECOND  ENTRY  POINT) 

KEMH0140 

LA 

3 , C LGAD 

SIT  BASE  REGISTER  FOE  CLOAD 

MEMR02CC 

BALR 

2.3 

LINK  TC  NOELEY  IECCM  ACCON 

MEMR021C 

Ltl 

14,3,12  (13) 

RESTORE  REGISTERS 

MFMRC220 

SR 

15,15 

SUPPRESS  VARIAELE  FETIJFN 

MEMB0230 

BR 

14 

RETURN 

MEMR0240 

DROP 

15 

MEMR0250 

• 

SECOND  ENTRY 

POINT.  IECCM  ENTERS  AT  CCFE2  THINKING 

IT 

* 

MEMR0260 

♦ 

WENT  TO  F1CCS.  THIS  ROUTINE  SIMULATES  FIOCS  BY  POINTING* 

MEMR0270 

CORE  2 


OUTPUT 


RETURN 


CHOVE 

* 

* 

* 


TO  BUFFER  ADDRESS  AND  LENGTH  STORED  EY  FIRST  ROUTINE. 
IBCCM  IS  RESTORE!  TC  NORMAL,  FOLLOWED  EY  RETURN  TO 
IBCCN.  A WRITE  BUFFER  IS  IM1ALIZED  TC  BLANKS  BEFORE 
IBCCM  FILLS  IT  TC  ALLOW  T FCRFAT  TC  WORK  CORRECTLY. 


USING 

*,1 

ST 

4, SAVE4 

SAVE  R4. 

LR 

4,1 

R4  = A (CCFE2) 

USING 

COR  E 2 , 4 

DROP 

1 

LR 

1,0 

R 1 POINTS  TC  FICCS  CAIL  PARAMETERS 

TM 

1 (1)  ,X'  OF  • 

TEST  ECR  OUTPUT,  FIRST  TIME 

BO 

OUT  F UT 

EFANCH  TC  FIRST  OUTPUT  ROUTINE 

L 

1, VFIOCS 

H1=A(E1CCS)  TO  RESTORE  IECOM 

LA 

3, CLOAD 

SET  BASE  REGISTER  FOR  CLOAD 

BALR 

2,  3 

LINK  TC  MODIFY  IECOM  ADCON 

LM 

2, 3 , EUFADfi 

ICAD  ARRAY  ACF  AND  LFNGTH 

B 

RETURN 

EFANCH  TO  RETUFN  TO  IECOM 

LM 

2, 3.EUFACR 

L CAD  ARRAY  ADR  AND  LENGTH 

M V I 

0(2)  , X 1 4 0 1 

BLANK  FIRST  EUFFER  LOCATION 
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